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ABSTRACT 


Studies of the adsorption of the secondary alcohols, 2-propanol, 
2-butanol and 3-pentanol, on y-alumina and y-alumina doped with 2% 
sodium hydroxide were carried out from room temperature to 300°C using 
infrared spectroscopy. On y-alumina four types of surface species 
were observed: (i) alcohol hydrogen-bonded to surface hydroxy] groups; 
(ii) an alkoxide-type structure chemisorbed on Lewis acid sites, alts; 
(iii) a partial double bond type structure which on desorption gives 
an olefin; and (iv) a carboxylate species which is stable even at 
high temperatures. The physically adsorbed species are mainly dehydrated, 
forming a mixture of olefinic isomers. On 2% sodium hydroxide-doped 
y-alumina, in addition to the species i, ii, and iv, a ketone-like 
species was detected which on desorption gave a ketone. Two types of 
carboxylate species may exist on sodium hydroxide-doped y-alumina. A 
mechanism for the formation of carboxylate species, identical for both 
the catalysts studied, is proposed. The three alcohols follow similar 
mechanisms for dehydration and dehydrogenation, with 3-pentanol re- 
acting at a lower temperature than encountered with 2-butanol and 2- 
propanol. 

The catalytic activity of the y-alumina surface has been in- 
vestigated, using the IR absorption technique, in terms of strong and 


weak Bronsted-acid sites and Lewis-acid sites. Doping of the y-alumina 
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with NaOH was used to eliminate Bronsted-acid sites. The influence of 
the two types of sites upon catalytic activity and selectivity was 
shown by the simultaneous dehydration and dehydrogenation of 3-pentanol 
on y-alumina or on several doped aluminas of varying NaOH content. The 
selectivity exhibited towards the two reactions was shifted from de- 
hydration to dehydrogenation by eliminating Bronsted-acid sites. 
Finally, simultaneous kinetic and mechanistic studies have 
been performed by IR spectroscopy to provide both quantitative con- 
version measurements and spectra for the species adsorbed/reacting on 
the surface of the solid catalyst. A recycle reactor was designed 
so that steady-state measurements could be taken of either the recir- 
culating gas concentration or complete spectra scans through the 
catalyst pellet. On the basis of simultaneous spectral and kinetic 
data, the kinetics of 2-propanol dehydration on y-alumina and of 2- 
propanol dehydrogenation on 8% NaQH-doped y-alumina were modelled 
adequately with the appropriate Langmuir-Hinshelwood equations, in- 
ferred from the surface mechanism, with a minimum of statistical 


correlation effort. 
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CHAPTER I 
INTRODUCTION 


Classically, the measurement of the kinetics of a homo- 
geneous chemical reaction has been used to deduce the mechanism of the 
reaction. In heterogeneous catalysis, it is well-established that the 
solid-catalyzed reaction proceeds by the adsorption of reactants on 
"active sites" on the surface, reaction of the adsorbed species with 
either another adsorbed species or a molecule from the surrounding 
phase, and finally desorption of the products to regenerate the surface 
active site. Accordingly, from observations of the product distribu- 
tion and the reaction kinetics, reaction mechanisms have been pro- 
posed. One major difficulty that arises in such a procedure is that 
more than one set of mechanisms can be proposed that fits all of the 
observed data quite well. The differences in fit may be so slight as 
to be explainable entirely in terms of experimental error. Therefore, 
if a number of alternative mechanisms fit the data equally well, the 
selected rate equation is one of good fit but not necessarily one that 
represents the true mechanism. 

The reason for searching for the "actual" mechanism is that, 
in reactor design, extrapolation to new or more favorable operating 
conditions is much more safely done. Since a study of kinetics alone 
can not always provide a conclusive reaction mechanisms, the mechanism 


problem remains to be solved. 
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Infrared spectroscopy has proven to be of great value in 
determining the structure of molecules. In recent years, its applica- 
tion to the study of surface chemistry has provided one of the most 
direct means of deducing the active sites on the catalyst, the inter- 
actions and perturbations that occur on the surface during adsorption, 
and the structure of the adsorbed species. Such results, would be 
anticipated to be very useful in inferring mechanisms of solid- 
catalyzed reactions. 

Although it was reported!) that the adsorbed species ob- 
served by infrared analysis may be good candidates for consideration 
as intermediates of the reaction, it is sometimes difficult to make 
direct observations on reaction intermediates in the infrared spectrum 
because they may only exist in small quantities. For a reaction inter- 
mediate to be detectable in an infrared spectrum, it must have a re- 
latively high surface concentration. This would be expected to be 
true only in cases where the slowest reaction step follows the forma- 
tion of the intermediate. Therefore, the possibility exists that less 
strongly chemisorbed species, subsequently existing in low surface 
concentrations, may be the reactive intermediate which may not be de- 
tectable in the infrared analysis. A combination of both kinetics and 
infrared spectral analyses would be thus expected to provide more in- 
formation for obtaining reaction mechanisms than by either method alone. 

In a few cases (2°3) , infrared study of adsorption and surface 
reactions were used in conjunction with independently obtained kinetic 


data to suggest kinetic models which were mechanistically sound. How- 























-. - a, 
ai a ae 
F LP in be 

a 

a 7 7 , 

nt sulsv Jsexp to od ot nevorg 26ri vaoszortaade bevsrtal wr 

a ms 
~s9f faqs 237 ,2169y Ins087 nl .astusefom to esd ound si pntntn v2 eb 

tzom 91 to sno bebfvotq esd yatetmads sostywe To vbute: ‘ont! ot or: 
“iint sty .Jeyfstso ont no 2adfe sviio& sAd pntoubeb to ensser Fo: 
‘ DA + 

Hottqro2zbs pniwb sostiwe sit no Was0 t5Ad3 ONS BeNORG brs enotios 

sd bDiluow ,esfueot owe .2efasqe bedvoebs shy to. eiutounde os 7 

-bifoz to emetnsizen BArry stint nt futsev yrev ed: ot bategt » $n, 


—_, 


, znotsaeen best a | 


‘ : - - ; ( [ }, ; : ar . 
~do asineqe bedigebs aft edt.‘ ‘betyoasy 25w If iguons TA . oo 
Fi : : iy ny aw -* 
Norteysbfenos yot 2otsbtbns2 boop 3d vem si bovs7tat. yd be ve 9 
- we — 


Ae 


16m of $fuoetttb eat: amoe ef tt Lnotdaset Sit 70 eatotbematat it 2 
muiso9qe bariettnt edd nt 2s tisarred it fortosey no eno tevrsedo dost 
age oe: 

-Y¥sint nottonst § 107 .2stst 'SRGUP Iism2 nt t2etxs vito yan vents mn 
~ : seal 
" > 
-91 6 sv6ed Jeum St ,muydo5qe ber ett a fs nf aidetoaseb: ad of on rb 
sd oF bstosqxs sd bluow etdTt .notdeevinssnoo sostwe eit wlovt 
sarvot sit ewol{ot qes2 nofsosat teswole edt evanw 29269 vt Uno uy: 
a " + ? a, ‘, Ms - 
ees! add eterxs yitltdte2oq shit ..stotarshtT ots thematar ond hy of; 


sastiue wot nr prtiel [x xisnaupaedue  estoaqe bedvoztmsds vipne 






-9b sd Jon ysm datdw sis theanatat. s aviszn6ey sit ad oon : ‘snot eam ~ 
r iy aftontA djod to pottsntdmos A 2hey Fens bovsttnt ony ake ‘ds at a4 
es ag, : ae 

“ oak som sbiverg bede bases aurlt sd bf fins mtdictsig Atos 3s banet 


non 
-aneinstioem noriossy entatstdo x0t not 
eee an a (63) gus 


ever, this approach has limitations in that different forms of the 
catalyst, different pretreatments of the catalyst, and different 
temperatures and reactant pressures may be necessitated by the two 
different sets of experimental conditions for the kinetic and infra- 
red spectral studies. 

To date, spectral analysis has been combined with simultaneous 
kinetic studies where both infrared spectra and kinetic data are ob- 


(4-8) | However, it was shown?? that ob- 


tained under batch reactions 
taining both infrared spectra of adsorbed species and kinetic data 
simultaneously under steady-flow conditions facilitates the determin- 
ation of diffusion limitations and eliminates complications of changing 
mechanisms which may occur under batch reaction conditions. Goldsmith! !°) 
and Baddour et ai'9) have applied simultaneous kinetic and spectroscopic 
studies in a flow system to the study of carbon monoxide oxidations on 
silica-supported palladium. However, their samples of the exit stream 
were analyzed on a gas chromatograph and steady-state operation was 
assumed when there was less than 2% variation in the outlet composition 
over a 30-minute period. 

The purpose of this thesis has been to develop and to apply 
a technique for the simultaneous measurement of the reaction kinetics 
and the spectra of the adsorbed surface species by using infrared 
analysis for a solid-catalyzed chemical reaction. This information 
was then to be used to identify the kinetically important surface 


species, to determine a realistic mechanism for the reaction, and then 
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to elucidate the kinetic models for various catalyst modifications. 
Steady-state operation can be observed directly by setting the infrared 
spectrophotometer at a single frequency of interest and by recording 

the changes in the sample beam intensity as a function of time! 1), 

This technique appears to be simple, versatile and accurate in obtaining 
the kinetic data. 

The heterogeneous alumina-catalyzed dehydration of alcohols 
has been known since the end of the eighteenth century. Although the 
dehydration reaction is generally agreed to involve surface acidic 
sites, widely different views are still held concerning the mechanism 
of this reaction, The use of y-alumina as a catalyst also has been 
studied intensively for many decades yet, the nature of the sites re- 
sponsible for their activity remains controversial. Earlier, in a 
kinetic study of dehydrogenation of 1-propanol over solid catalysts, 


Masudaya\ 2 


reported that the role of the alumina was altered from 
that of a dehydrating catalyst to one of the dehydrogenating catalyst 
by doping with sodium hydroxide. The relation of the catalytic activity 
and selectivity to the supposedly responsible acidic sites on the y- 
alumina has not been clearly revealed. For these reasons, the decom- 
position of secondary alcohols, C, to Ces on the dual function catalysts 
made of sodium hydroxide-doped y-alumina was chosen to study the re- 
action mechanism using the proposed infrared technique. 

Some of the spectral analysis in this study was obtained with 


a Perkin-Elmer Model 621 infrared spectrophotometer connected on-line 


to an IBM 1800 computer. In the study of adsorbed species, elimination 
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of the background spectra of the catalyst by obtaining infrared data 
on-line saves a great amount of work and of time. The system components, 
the operating procedures and data processing methods are described in 


Appendix A. 
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CHAPTER II 
SURVEY OF RELATED WORK 


2.1 Dehydration and Dehydrogenation of Alcohols 

The dehydration of simple alcohols on alumina catalysts leads 
mainly to the formation of ethers, olefins andwater. At lower temper- 
atures, the dehydration of unbranched alcohols (up to 1-hexanol) yields 


ethers, olefins and water’ 2! 


whereas for secondary alcohols ether has 

not been detected as a dehydration product. Product distributions for 

the dehydration of 2-propano] ‘!4-17) | 2-butano] 14215518) and 3-pentanol ‘2>!8) 
have also been investigated. The dehydration of 2-propanol leads to 

the formation of propylene, with trace quantities of isopropyl ether 

at low reaction temperatures. The dehydration of 2-butanol yields 

]-butene and, cis- and trans- 2-butenes at al] the reaction temperatures 
studied. The dehydration of 3-pentanol gives predominantly 2-pentenes 


(18) | The cis 


and traces of 1-pentene according to Pines and Manassen 
to trans ratio of the 2-pentene found is 2.5 and 1-pentene is thought 
to be a secondary product by isomerization of 2-pentenes at reaction 

temperatures above 300°C. Thus the dehydration of secondary alcohols, 


C4 to Ce can be represented as follows: 


Al503 


Ree Cio Cua 
\ + Ry aC Shen Ro rr HA0 : (221) 


2 
OH 
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The dehydrogenation of alcohols on sodium hydroxide-treated 
alumina catalyst has also been studied (2>19) | In this case the de- 


hydrogenation reaction proceeds as follows. 


Ry = or - R 


OH 


2 NaOH+AI,0, 
~ OR, - f sae (2.2) 


The catalyst used in their studies was y-alumina treated with 2N 
sodium hydroxide aqueous solution. Besides ketone and hydrogen, some 
dehydration products were found. 


(19,20) 


In some studies » the dehydration and dehydrogenation 


of alcohols accompanied by thermal decomposition were also observed. 
However, the thermal decomposition can be avoided by studying the re- 
actions at lower temperatures (below 450°C) in a glass reactor. 

A comprehensive survey of the alcohol dehydration was pub- 
lished by Winfield (2!) for work prior to 1960. The two major issues 
seem to be the mechanism and the kinetics of the reaction, and the 
type of active sites of the catalysts. To date, these questions are 
still the subjects for several investigations. 

The catalytic activity of a surface is strongly dependent on 
its structure and also its properties in addition to the surface com- 
pounds of the substrates and reaction products. Henceforth, the 
eeruetine and properties of the y-alumina surface, the adsorption 
structures of both alcohols and their dehydration as well as dehydro- 
genation products will be reviewed before the kinetic and mechanistic 


Studies are described. 
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2.2 The Structure of y-Alumina 


The y-alumina is one of the transition aluminas most commonly 
used for catalytic purposes. It is formed by controlled thermal de- 
composition of the hydrate, bayrite, or of the oxyhydrate, boehmite, 
at 450°C. Above 600°C, however, transformation to the delta phase 
occurs. The thermal process can be represented by the following 
sequence: 


450°C 600°C 1050°C 1200°C 
Boehmite > y-Al,0, > 5-Al,0, > 6- and a-A1,0. * a-A150. 


The crystalline structure of y-alumina is a defect spinel 


lattice (22), 


The unit cell of the crystal consists of 32 oxide ions 
with 21 1/3 aluminum ions arranged at random in the 16 octahedral and 

8 tetrahedral positions of the spinel structure. Ina single-crystal 
electron diffraction study, Lippens and deBoer (23) found that the oxygen 
lattice of y-alumina should be fairly well ordered. Consequently, the 
disorder in the y-alumina is determined mainly by disorder of the 
aluminum ions. Saalfeld and Mehrotra ‘24) determined the cation dis- 
tribution by Fourier synthesis of the electron diffraction patterns 

and found that the octahedral aluminum sublattice was fully occupied 
and that the necessary vacant sites were randomly distributed over the 
tetrahedral interstices. The details of the crystallographic structure 
determination have been discussed by Lippens (25) . 


The necessity for maintaining the structural integrity of the 
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catalyst throughout the investigation is likewise apparent. Often the 
change in water content of the catalyst is accomplished by treatment 
at various higher temperatures; care must be taken during this treatment 


(26) studied 


to eliminate possible phase transitions. Pines and Haag 
the effect of both the treatment temperature (up to 1000°C) and the 
water content upon alumina activity towards isomerization of cyclohexene. 
They found a pronounced decrease in activity between 800 to 900°C, which 
could be the result of a phase transition of the catalyst (27) since 
delta alumina forms at 600°C. Care must be taken not to alter either 
the phase or modification of an alumina since by doing so one con- 


currently changes its surface site configuration and thereby its in- 


trinsic activity. 


2.3 The Surface of y-Alumina 
Though many authors (28,29) believe the [111] face of the 


oxygen sublattice to be energetically favored, Lippens (25) , based on 
x-ray and electron diffraction studies, concludes that the [100] plane 


(30) that 


of y-alumina is preferentially exposed. Since it is believed 
the surface is terminated in anions, which are more readily polarizable 
than cations, the outermost aluminum ions thus Jie below the surface 
plane. These aluminum ions provide the electrophilic centers by which 
water is dissociatively held on the surface; the electron pair of the 
OH eraup being positioned over the aluminum ion while the dissociating 


proton forms a hydroxyl via binding to lattice oxygen (31) | This re- 


action proceeds until the surface is completely filled with OH groups. 
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At higher water vapor pressures more water is bonded in a multilayer 

physical adsorption process but this water can be removed easily by a 
drying procedure at about 120°¢'32) | Upon heating at higher temper- 

atures, the opposite reaction occurs, elimination of water by re- 


combination of the two OH groups, which leaves a lattice oxide jon and 


an aluminum ion site at the surface in the following manner'33»34) , 
H H 
/ ‘p 0 
0 0 ae US a (2-3) 
| | Al Al 
Al Al 


The aluminum ions thus formed may be regarded as acidic centers of the 


(35,36) (37,38) | 


Lewis type which can be selectively poisoned by organic bases 


However, the OH groups were reported to be very weak acidic Bronsted- 
(37,39,40) 


(35) 


acid centers which can react with ammonia (pK, = 5) to form 


ammonium ions but not a strong enough acid to react with pyridine 
(pK, = 9) to form pyridium ions (37) , Tamele(27) depicted the attach- 
ment of water to the alumina as the completing of the octet of electrons 
of aluminum and the conversion of a Lewis-acid site to a Bronsted-acid 


site as shown in reaction (2.4) 


H.-H 
07 TIAL pi ite that 
ie oH 0 es alee A ee (2.4) 
0 7° 
| 508 


An adjacent lattice oxygen provides an electron pair for sharing with 
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a proton. In this way two hydroxyl groups are created. 


H..H 
:0: H 
Chad Ali tue: Qh te Gdenvblaann 0 (2.5) 
0 a 


The lattice OH groups thus formed are therefore positionally below 
the plane of the parent OH groups. 
In recent years, the most extensive study of the alumina 


surface was made by Peyj (3994142) 


using gravimetric and infrared 
Spectroscopic methods. For an undried y-alumina, infrared bands were 
observed near 3300 and 1650 cm”! and these bands were assigned to 

the stretching and bending frequencies, respectively, of molecular water 
adsorbed on the alumina surface. After evacuation at 720°C, the spectrum 
has shown three well-defined peaks at 3698, 3737 and 3795 em C41). 
These bands are essentially symmetric and are attributed to the stretch- 
ing vibrations of isolated surface hydroxy] groups. Deuteration causes 
these bands to disappear and to be replaced by three equivalent bands 

at 2733, 2759 and 2803 cm), All bands are shifted in frequency by a 
constant factor of 0.738 which is close to the theoretical value (0.73) 
for the isotope shift. This provides good evidence that the original 
absorption bands were due to the stretching vibrations of surface 


(41) ] 


hydroxyl groups. Also, it was observed that the 3698 cm hydroxy] 


groups were able to exchange their H for D faster than the other two 


groups and therefore the groups corresponding to the band at 3698 cm” | 
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were proposed at be the most acidic of the three hydroxyl groups. 

The aluminum ions exposed by vacancies in the oxide lattice 
were also considered by Peri (39) and others (43) as Lewis-acid centers 
which provide the principal sites for the adsorption of ammonia mole- 
cules (a Lewis base). Thus, ammonia should react with the alumina 


Surface to give a chemisorbed species as follows: 


—9 0 H 
“x ~ | 
—O— Al + NH, ==-0—Al + N—H (2.6) 
y | 
—0 -0 H 


However, ammonia reacts with the alumina via the chemisorption process, 
reaction (2.6), and other processes as well. The adsorption is compli- 
cated by interactions of the ammonia with surface hydroxyl groups and 
surface oxide ions, and thus does not provide immediate evidence of 
these Lewis-acid sites. Pyridine was found not to react with hydroxy] 
groups and oxide jons on the alumina surface. Therefore, the spectrum 
of adsorbed pyridine can be interpreted as showing the presence of 
these Lewis-acid centers. The details of the analysis have been given 


(37) and no further discussion is attempted here. 


by Parry 
During the infrared study of the adsorption of ammonia on 
the alumina surface, Peri (39) also reported that a new hydroxyl band 


appeared at 3725 cm! 


This was considered to indicate the presence 
of the oxide ions on the alumina surface. The chemisorption reaction 


can proceed as shown in reaction (2.7). 
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ff once. (207) 


Peri (35) 


proposed that the above reaction as a completely ionic pro- 
cess. Other workers (44) have often represented the bonding between 
atoms in the surface in covalent terms. In this way reaction (2.7) 


could be written as follows: 


0 H 
Re Sep TES heat / 
A] BY oe NHS On eH (2.8) 
| ) 3 | 2 
Al Al 
Le ES 


The true nature of the bonding in oxides probably lies between the 
extreme ionic and covalent descriptions ‘44) , 

Several interpretations of the nature of active sites on the 
y-alumina surface are available. One view proposed(4°) that the aluminum 
or oxide ions exposed at the surface act as active sites and that the 
surface was free of hydroxyl] groups. This view is contradictory to 


(41) 


that presented by Peri and Hannan who observed that alumina surface 


is completely hydrated at lower temperatures by investigation of water 
adsorption and analysis of the infrared spectra of water adsorbed on 
y-alumina. 

The other view proposed that the alumina surface contains 
cation defects and strained sites, which result from the dehydration 


46,47) 


of the original hydrated surface! A more detailed model was 


42) 


proposed by Peri who obtained the model by a computer simulation 


from the dehydration of a fully hydrated surface on (100) crystal planes. 
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The resultant computer model predicted a surface containing five dif- 
ferent types of hydroxy] groups. The various types of hydroxyl site 
are shown in Figure 2-1 and the corresponding infrared vibrational 
frequencies are shown in Table 2-1. This model may be considered more 
or less representative of the dehydration mechanism and the structure 
of alumina surface. However, the validity of the assumptions 
proposed by Peri (42) have been subjected to criticism by Lippens (33) , 
1.€. for alumina samples heated at temperatures above 600°C, the assump- 
tion that a spinel type alumina exposes only (100) planes on the surface 
is very unlikely energetically. 

Although Peri's model could be used to explain some of the 
experimental observations on the alumina surface, conflicting results 
were obtained in this study. These wil] be discussed later in this 


thesis. 


2.4 Effect of Sodium Hydroxide on the Structure of Alumina Surface 


In a study of reactions of 1-propanol on solid catalysts, 

Vasudeva‘ !2) Showed that the activity of alumina towards the dehydra- 
tion of 1-propanol was completely destroyed when the alumina catalyst 
was treated with 2N sodium hydroxide solution. It was attributed (> 19) 
to the destruction of the "acid centers" on the alumina upon treatment 
with sodium hydroxide. 

Recently, in an infrared study of l-propanol adsorption on 
sodium hydroxide doped y-alumina, Deo and Dalla Lana‘ 48) reported that 
a 2% of sodium hydroxide content in y-alumina was enough to react with 


] 


the surface hydroxyl groups of 3785 cm to form = Al - 0 - Na. Thus, 


bt 
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Figure 2-1. Computer Model of a Dehydroxylated 
Alumina Surface (42). The five types 
of OH are assigned the frequencies 
in Table 2-1. 
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Table 2-1 


IR Spectrum of Hydroxy] Groups on Alumina* 


—_—_—_—_—_—————————— 


Band Wave number, cm”! No. of nearest 
Oxide neighbors 


————_—_—_[—"_$—"—_——__ 


A 3800 4 
B 3780 3 
C 3744 2 
D 3733 
E 3700 0 


*With reference to Figure 2-1 
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17 
the reaction was considered to be with the Bronsted-acid sites. Un- 
fortunately, no data were presented for the effect of sodium hydroxide 
on Lewis-acid sites, and the complete reaction mechanism can not be 
obtained at this stage. However, Pines and Haag (2°) suggested the 
possibility that doping with sodium hydroxide could give rise to an 
acid/base surface complex of the type shown in reaction (2.9) 


| 
— Al +NaQH + [-Al— OH] Nat. (2.9) 
| 


This reaction is not in agreement with the data reported by Ross 


(49) 


and Bennett who found no evidence for the formation of such a 


complex in any detectable amount. 


2.5 Infrared Studies of the Adsorption of Alcohols on y-Alumina Catalysts 


Infrared spectral studies of the adsorption of alcohols on 
alumina surface have been previously reported and four different types 
of surface species were observed: 

1. physically adsorbed alcohols via hydrogen bonds §45»48,50) 

2. an aluminum alkoxide type of structure formed by dissociative 
chemi sorption (45»48,51-56) | 

3. a carboxylate structure formed by oxidation of the 
alcoho] (45»4851-53,57 58) 

4. adsorbed water, one of the dehydration products. 
It was reported ‘®9) that the adsorption of 3-pentanol on the sodium 


hydroxide doped y-alumina at higher temperatures yields the corres- 


ponding ketone and hydrogen. Therefore, the role of adsorbed ketone 
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and hydrogen will also be discussed. 
2.5.1 Physically Adsorbed Alcohols 

The infrared absorption of free hydroxyl groups of alcohols 
should be expected to appear in the range of 3600 to 3650 cm”! , e.g. 
spectra of vapor phase alcohols. The OH stretching vibrations shift 


] 


from ~ 3600 to 3350 cm for all the liquid phase alcohols (C, to C,). 


This frequency shift can be explained by the formation of hydrogen 
bonds and the extent of the shift is a measure of the strength of the 
hydrogen bonds. All alcohols adsorbed on alumina give a strong OH 
stretching bands centered around 3500 cm”! This can be interpreted 
to be the formation of weak hydrogen bonds as compared to the 


3350 cm! for liquid alcohols. 


(54) (56) 


Whereas Babushkin et al and Boreskov et al thought 


that they could still detect physically adsorbed ethanol even at high 


(45,51) ( 


desorption temperatures, Greenler et al methanol and ethanol), 


53) (ethanol, 1- and 2-propanol), kage) '45) 


Treibmann and Simon! 
(1-propanol and 1-butanol), and Deo and Dalla Lana (48) (1-propanol ) 
found that alcohols are readily desorbed on evacuation even at room 
temperatures. NMR studies showed that the hydroxyl group of the 
alcohols interacts directly with the alumina surface 00) it must 
therefore be assumed that hydrogen bonds are formed. Kage (45) pro- 
posed, without direct spectroscopic evidence, that the structures, 


[I] and [II], are present on y-alumina surface. 
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Although Knozinger et al reported the presence of structures [I] 
and [II], Deo and Dalla Lana‘ 48) observed very little hydrogen-bonded 
species when hydrogen atoms of the surface hydroxyl] groups were re- 
placed by sodium atoms. This, therefore, indicates that structure 
[II] is a more probable representation of the hydrogen-bonded alcohols. 
This is also in agreement with the results obtained by Kryloy'61) | and 
Jain and Piliai (62) | 
2.5.2 Aluminum Alkoxide Structure 

Ethanol, on adsorption at lower temperatures (20 to 100°C) 
gives rise to a surface compound whose infrared spectrum is almost 
identical with that of the corresponding aluminum alcoholate >!) , For 
comparison, both spectra are shown in Figure 2-2. Similar results have 
been obtained with 1-propano] (45253) | 2-propanol (93) , 1-butanoy (45) | 
iso-butano1 (94) , tert-butanol° and benzy] alcohol '°2). The dif- 
ferences in the infrared spectra between surface alkoxides and aluminum 
alcoholates are attributed to the different ratios of - OR radicals to 
Al atoms. 


There are two possible ways for the formation of alkoxide 


species by dissociative adsorption: 
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R H 
ae ee 
0 0 0 
it a ] | 
Al Al + ROH > Al + Al (2.10) 
Ne Ge aN Bees 
H 
Os 
0 
| { 
Pie + ROH > Al - OR + HA0 (2. 11) 


In reaction (2.10) a hydroxy] proton is split off from the 


alcohol to form a alkoxide structure. The proton then adds to a 
Surface oxide ion to form a hydroxyl group. Reaction (2.11) involves 
a hydroxy] group from the alcohol and a proton from the surface hydroxy | 
group. The C - 0 bond of the alcohol is broken and the resulting 
hydroxyl group reacted with a surface proton to form water. 

Although the somewhat lower energy of rupture of the C - 0 
bond in the alcohol favors the second possibility, reaction (Zak) 
the lattice energy liberated on incorporation of the alkoxide or hydroxy] 
group favors the rupture of the O-H bond !©3) | Treibmann and Simon (23) 
concluded from displacement reactions that the alkoxide formation in- 
volves incompletely coordinated aluminum ions. Kagel ‘45) reached the 
Same conclusion on the basis of mass spectroscopic residual gas analyses 
after adsorption. 
2.5.3 Carboxylate Structures 

Oxidation of the alcohols during adsorption on alumina above 


170°C to give carboxylate structures has been observed for methanol {51 »58) , 
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(51,53 ,56 ,58) 45 48,64) 


ethanol ‘ 1-propano1‘ 1-butanot (45 +64) and 


2-propanol (93958) | 


Also corresponding surface compounds have recently 
been found on the adsorption of benzyl] alcohol and isobuty] alcoho1 (65) 
on alumina above 120°C. 

In the study of adsorption of methanol, Greenler (>!) observed 
that on heating the methoxide to higher temperatures the spectrum 
shows three new infrared absorption bands at 1597, 1394 and 1377 cm, 
By comparison with the published spectra of various metal formates , 
this indicates that the surface methoxy group has decomposed and pro- 
duced a surface formate [H - oa ] . The formate ion may be treated 
as a simple XZY., molecule that gives rise to a symmetric and an asym- 
metric C - 0 stretching vibration. For aluminum formate these are 
found at 1370 and 1600 cm! respectively. Consequently, the bands 
found at 1377 and 1597 cm! in the spectrum of the surface species 
are attributed to these modes. The third absorption band at 1394 cm”! 


| for the aluminum formate are 


for the surface species and at 1390 cm 
attributed to a C - H bending vibration. Further conformation of the 
nature of this species was obtained by isotopic substitution of cls 


for c!@ 


in the reactant alcohol. The resultant spectrum showed that 
the isotopic shifts were in good agreement with calculated values. 
Thus, interpretation of the nature of the surface species is considered 
to be satisfactory. 

From an extensive infrared study, Fink (98) reported mechan- 


isms for the formation of the carboxylate compounds as follows: 
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OH 
C,H-OH + CH.CHO (adsorbed) + surface acetate + H 


ole 3 teale) 


2 


OH 
CH,CHOHCH., > CH,COCH., (adsorbed) + surface acetate + CH, (S13) 


where OH denotes surface hydroxy] group. 

All the carboxylate structures are observed at temperatures 
up to 600°C. In view of its high stability, this structure could 
hardly be important as an intermediate in surface reactions which 


occur below 300°c (57 66) | 


However, it was reported (67) that, above 
300°C, the dehydrogenation of ethanol on alumina could proceed via an 
acetate structure. This is not in agreement with the experimental 
results reported by Deo and Dalla Lana‘48) | 
2.5.4 Adsorbed Water 

The adsorption of alcohols below the reaction temperature is 
always accompanied by the appearance of a band at 1600 cm! when the 
infrared spectrum is scanned at room temperature. This band increases 
in its intensity with increases in the adsorption temperature and it 
is not believed to be due to the alcohol or any of its surface com- 
pounds but due to adsorbed molecules of water (68), 

It is generally assumed that molecules containing oxygen 
and having two pairs of free electrons are preferentially adsorbed 


via two passive hydrogen bonds as in structure crm '44) 
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This type of adsorption bonding via two hydrogen bonds requires two 
adjacent hydroxyl groups on the surface. It is believed(©9) that at 
200°C, there are still hydroxy] groups covering 50 to 80% of the 
alumina surface and thus structure LIII] represents the adsorption 
of water from room temperature to 170°C. 

However, adsorption on isolated hydroxyl groups such as 
structure [IV] should be possible, particularly when two adjacent 


hydroxyl sites are not available. 


[IV] 


TALAS LE 


< 


It was also proposed ‘©3) that structure [IV] is more probable when 
water is adsorbed in the presence of excess alcohol, e.g. adsorption 
of alcohols at higher temperatures. 
2.5.5 Adsorbed Ketone 

It is well established in the literature that ketones wil] 
form stable complexes with inorganic Lewis-acid halides such as AICI. 
and BCI. The complex formation takes place via the oxygen atom of 
the carbonyl group, and the donation of electrons from the oxygen atom 
into the vacant orbital of the Lewis acid leads to a delocalization of 
the electrons in the double bond'70) | In the case of alumina, Lewis- 
acid sites are believed to be present on the surface and thus a similar 


reaction should occur: 
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I 
0 (2.14) 
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When acetone is adsorbed on alumina surface, a strong shift in infrared 
frequency of the free carbony] group is observed 98) , The type of the 
chemical bond of the chemisorbed complex is similar to that of electron 
donor-acceptor complexes. For comparison, infrared absorption peaks 


of acetone in different states are listed in Table 2-2. 


Table 2-2 


Infrared Frequencies (cm!) of Acetone in Different States 


eee 


: Acetone Acetone on Acetone on 
Assignment Sire 
liquid BCI, Solid(7!) y-aq,0, (58) 
| Zo 
eee ee ee eS 
v0 = 0 1710 1625 1620 
6 CH, 1433 1400 1460 
545CH3 1358 1362 1380 
vo =C 1225 1262 1250 


—_—_—— OC Cr 


2.5.6 Adsorbed Hydrogen 
So far, very little attention has been paid to the study of 


adsorption of hydrogen on y-alumina. Russell and Stokes (72) reported 


that hydrogen is adsorbed on alumina surface at temperatures above 
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500°C. According to these authors the amount of hydrogen adsorbed is 
small, comprising less than 0.2 cm?/g (S.T.P. Ai4No adsorption mechanism 
was given. It is expected that molecular hydrogen can be dissociatively 
adsorbed on y-alumina only with difficulty because the free energy change 
for the dissociation reaction is about 200 K cal/mole (72), However, 


Peri (35) 


and Duken et a1'73) Proposed that ammonia can be dissociatively 
adsorbed on y-alumina surface. The proton from ammonia reacts with a 
surface oxide ion to form a hydroxy] group. More recently, in an infra- 
red study of chemisorption of acetylene on alumina, Bhasin et ay (34) 
proposed that protons are accepted in the vicinity of OH groups form- 

ing hydrogen bonded groups (Reaction (2.15)), and to some extent by 


neighboring oxide ions to form isolated hydroxyl groups (Reaction 


(2.16)). 


\ 
H 
H oe 
pas 0 
iba: | x Per (2.15) 
Al Ae 
aN ee man eee 
0 Sa nee 
H 
0 vA 
whe 0 
H+ Al ees ean (2.16) 
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ES PAK 
They claimed that this is evidenced by the rapid appearance of the 


broad band at 3600 cm™! due to hydrogen-bonded species and Significant 


increase in the intensity of the infrared bands due to the surface 
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27 
hydroxyl groups. 


Several investigations (©»45+48.51) 


elucidating the mechanism 
of the formation of alkoxide from alcohols have been reported. It is 
generally agreed that the proton, which splits off from the alcohol, 


is adsorbed on the surface oxide ions. 


2.6 Mechanism of Alcohol Dehydration on y-Alumina 

Many reaction mechanisms have been postulated for the de- 
hydration reaction which differ primarily in the intermediate com- 
pounds. They can be classified as follows: 

(a) Pines and Haag’’*), Kryloy(®!) and Jain and Pillai (©) 

prefer an adsorption of the alcohol through the oxygen of the alcoholic 
hydroxyl group on incompletely coordinated aluminum ion. The C - OH 
bond of the alcohol can therefore be polarized, making the OH group 


easier to remove. A proton from the g-carbon is thus abstracted by 


an oxide ion to form water, i.e. 


H oH 
| | 
Rj— C—C—R, , 
eo : Ric ee Serer Papen eat) = 
10 HN Ry - gu glio) et 
My |’ 
ce (aula) 
ety er 


Although, the participation of basic sites (oxide ions) in the de- 
hydration of alcohols on alumina catalyst has recently been proved by 


poisoning experiments with tetrocyanoethylene‘/°), the participation 
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28 
of the aluminum ions has not been proved. 

(b) In an infrared study of the adsorption of ethanol on y- 
alumina, Bremer et a1 (57) observed an infrared absorption band at 
1080°cm™! and attributed this band to the C-0 vibration which is due 
to the alcohol coordinately bound to the Lewis center (aluminum ion). 
Using the principle of the least structural change, which has been 
applied in the past to the proton exchange, they proposed a reaction 


mechanism for the dehydration of 2-propanol on y-alumina as follows: 


H H H H 
| | ' | | 
eit Gicohes . 12) telg-fetano ‘ 
Gold 0" 4H 
of H As H 
H 0” Ly 0~ 
| 
—0— Al-0 — AI— 0 — —0—Al — 0 —AI— 0— 
a 
Habactton Oral H3C — CH - CH, 
fast H - H H H ee 
y~! na (0 as kell 


| H 
y — —- 


Hy 


0 + CH, - CH = CH 


I + H,0 (2.18) 
—0— Al— 0— Al— 0— 


2 2 
The mechanisms (2.17) and (2.18) are similar in that both reactions 
involve Lewis-acid sites (aluminum ions). However, by poisoning these 


sites on alumina surface selectively with pyridine, many authors (63+76 377) 
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29 
concluded that Lewis-acid sites were not involved in the dehydration 
reactions. 

(c) The formation of olefins by dissociation of surface alkoxide 
groups has been discussed in literature (6931.78.79) | In a study of 
alcohol dehydration on silica-alumina, DeBoer and Visseren (2) re- 
ported that the reaction was not governed by the acid character of 
the catalyst but by the total number of aluminum atoms Present on the 
catalyst surface. They concluded that the alcoholate formation may 
be the intermediate process in the dehydration reaction. Surface 
alkoxide species are indeed formed during the adsorption of certain, 
but not all, alcohols. For example, tert-butanol which solely forms 
olefins on dehydration at lower temperatures (85 to 160°C) does not 


52). 


form surface alkoxide structure | The participation of the surface 


alkoxide species in the dehydration of secondary alcohols was re- 
jected recently(2!) , 

(d) A carbonium ion mechanism similar to that assumed for liquid 
phase reaction in acidic media is also proposed by several authors (8 »83) 
for the dehydration of alcohols on alumina catalyst. A proton from 
the surface hydroxyl group is assumed to transfer to the adsorbed 
alcohol molecule. The participation of surface proton (from surface 
OH groups) is supported by the need of water to be present on the 
catalyst (63) | However, Knozinger ©) reported that (CD) COD gives a 
high primary kinetic isotope effect, contrary to the expectations 
of a carbonium ion mechanism, in which the rate of C - 0 heterolitic 


fission should be rate-determining and should not therefore be in- 
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fluenced by substitution of hydrogen with deuterium. Moreover, for 
a carbonium ion to be formed by addition of a proton, the alumina 
surface must have contained some fairly strong acid centers of 
Bronsted type. The hydroxyl groups present on alumina surface were 
found to be only weakly acidic as discussed in Section 2.3, so that 
ionic addition of protons to alcohols is very unlikely. 

(e) From the above discussions, it seems that an alcohol mole- 
cule held to the surface by hydrogen bond is the most probable re- 
action intermediate. Thus the dehydration reaction, as proposed by 


48) 


Deo and Dalla Lana! » can proceed as follows: 
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| - H 
Lnyarracsite 
See + R, - CH= CH- R, + H,0 (2.19) 
a. 
H 0 
| 
pee Ah 


During a study of the adsorption of isopropyl alcohol on alumina, 


Yakerson et ay (84) 


pylene while alcohol bonded to a lattice aluminum atom forms acetone, 


with formation of a complex with acetate character as a side reaction. 


The rate-determining step can be established by measurement of the 


kinetic isotope effect. Investigations of this type have so far been 


published (85) for the dehydration of tert-butanol on y-alumina. No 
primary isotope effect was observed in the dehydration of (CH3) COD, 


suggested that hydrogen-bonded alcohol yields pro- 
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3] 
in which heavy water is formed on deuterated surface; the rate can not, 
therefore, be determined by the desorption of water or the cleavage 
of the R - H bond. On the other hand, a large primary isotope effect 
is found in the dehydration of (CD,),COD, thus the removal of the 
B-hydrogen in the Reaction (2.19) must be the slowest step in tert- 


butanol dehydration. 


2.7 Kinetics of Alcohol Dehydration on y-Alumina 


In order to understand the mechanism of a heterogeneous re- 
action, a knowledge of the reaction kinetics is necessary. However, 
only a few studies have been reported on the kinetics of the dehydra- 
tion of alcohols on y-alumina. In the following sections, in light 
of earlier investigations related to this study, the kinetics of the 
alcohol dehydration are discussed. 

The experimental results of Brey and Krieger °°) for the 
ethanol dehydration on y-alumina at 350 to 400°C, represent the 


following rate equation best: 


kL Salen 


22400 . rate = >> 
KPa ts KwPw 


(2.20) 
This equation is derived on the assumptions that surface reaction is 

the rate determining-step, and that ethanol and water are strongly 
adsorbed while ethylene is not. This model is essentially a Hougen-Watson 
type for a single-site mechanism with the assumption that the reaction 


is irreversible. 
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More recently, de Mourgues et ai (16) have studied the de- 
hydration of 2-propanol on alumina in a flow reactor at low temperatures 
(100 to 160°C) and low alcohol partial pressures (8 to 23 mmHg). In 
order to ascertain the inhibiting action of the reaction products, 
they introduced propylene as well as water into the propanol feed. 
Their results show that the reaction is zero order with respect to 
the partial pressure of 2-propanol and that the reaction is inhibited 
by water, but not by propylene. On the basis of these observations, 
they claimed that a mechanism with a Hougen-Watson equation and com- 
petitive adsorption on the same site does not fit the experimental 


data well. The following rate equation was proposed 





ae) 


where By and BY are the total number of sites which can be occupied 
by alcohol and water, respectively; On and Op, are the fractional 
coverages of By and Be respectively, and Kn and K are some constants. 

In a study of dehydration of isobutanol, Knozinger ‘®3) formu- 
lated an empirical kinetic equation in which the reaction rate was 
Proportional to the square root of the partial Pressure isobutanol, 
which is indicative of the dissociative adsorption of the alcohol. 
The equation is given as follows: 

ro 7Pp 


‘{.—:. (Puce) 
(YP, + bP.) 
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where Yo and r are reaction rates at time zero and time t. Pa and Be 
are partial pressures of alcohol and water, and b is a constant. This 
is in conflict with the reaction mechanism proposed by Knozinger in 
his earlier investigations. 

In later works on the dehydration of alcohols on alumina 
catalyst, zero and first order kinetic equations were hypothesized in 
some cases, or else no kinetic equation was given at all. 

Several investigations have used the rate equation of Hougen- 
Watson type to represent the reaction mechanism. From the kinetic 
data for the dehydration of 2-butanol, Kittrell and Mezaki (86) pro- 
posed a dual-site mechanism with surface reaction as the rate-con- 
trolling step. The same conclusion was reached by Carra et a1 (87) 
for the dehydration of cyclohexanol and by Bremer et ay (57) for the 
dehydration of 2-propanol. On the other hand, Butt et a1 (8889) | 


(61) and Kasaoka et ay (90) chose a single-site mechanism with 


Kry lov 
Surface reaction as the rate-controlling step for the formation of 


olefins from alcohols. 


2.8 Activation Energy for Alcohol Dehydration on y-Alumina 


In the literature, different results on the activation energy 
of dehydration of alcohols do not agree. The apparent activation 
energies for the dehydration of alcohols on y-alumina catalyst are 
summarized in Table 2-3. 

The difference in the two values of Reference (91) was ex- 
plained to be due to the difference in the heats of desorption of 


water at the two different temperatures, which is 13 K cal/g-mole. 
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Table 2-3 


Apparent Activation Energies for the Dehydration of Alcohols* 


Alcohol Temperature Activation energy References 
Range (°C) Kcal/g-mole 
2-propanol 350 16 (21) 
2-propanol 105-145 26 (91) 
2-propanol 170-200 39 (91) 
2-propanol 100-160 32 (16) 
2-butanol 350 14.9 (21) 
isobutanol 350 16.4 (21) 
isobutanol 160-240 30.0 (52) 


——————— eee 


*Activation energies are found to be independent of system pressure 
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The value of 16 K cal/g-mole obtained for the dehydration of 2-propanol 
at 350°C was explained in a similar way: at 350°C the rate of reaction 
can be approximated by a first order rate equation with respect to 
2-propanol. The apparent activation energy would thus be the dif- 
ference between the true activation energy (estimated to be 25 


K cal/g-mole 16), 


, and the heat of adsorption of 2-propanol (approxi- 
mated to be equal to the heat of liquification of 2-propanol, which 


is 10 K cal/g-mole). 


2.9 Effect of Doping Sodium Ion on y-Alumina on the Dehydration of 
Alcohols 


(78) 


Topchieva et al showed that the dehydration of ethanol 


was not affected by the presence of sodium ions whereas other workers (2649557) 
Claimed that the addition of sodium hydroxide suppressed the tendency 
of y-alumina to dehydrate alcohols. These investigators did not examine 
the influence of the sodium ion impregnation upon the catalytic selectivity 
of y-alumina. However, on impregnating alumina with KOH, Kuriacose et 

a1 (92) reported a decrease in the ratio of the rate of dehydration to 

the rate of dehydrogenation of 2-propanol from 12.20 to 1.75. Pisman 


94) also found evidence of 


et ay (93) and, Maatman and Vande griend 
1-butanol dehydrogenation over KOH- and KC]-treated alumina catalyst. 
Kinetics and mechanism of dehydrogenation of alcohols on the sodium 
hydroxide-treated y-alumina are not available in literature at this 


time. 
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2.10 Kinetic Model Discrimination Techniques 
The model discrimination procedure is often regarded as the 

Process of selecting the best model from the various eligible models, 
assuming an adequate model is included among the models proposed for 
the system studied. This problem has been studied extensively and 
useful techniques have been reviewed by several authors (95+96) | After 
the preliminary screening of the data to eliminate poor models, model 
discrimination consists of the following three steps: 

1. parameter estimation, 

2. computing the model discrimination index - likelihoods 
of the data, and 


3. sequential design of experiments. 


The techniques, which cover the above three areas, developed by Singh (98) , 
were adopted in this work. Details of the methods are Summarized as 
follows: 
2.10.1 Model Discrimination Procedure 

Several criteria have been used as indices of mode] perfor- 


mance ( 98) 


» the earliest of them being the minimum sum of errors 
Squared. However, one of the most successful criteria is the Bayesian 
expected likelihood procedure based on a Bayesian interpretation of 
probability. The Bayesian approach has been used for model discrimin- 
ation in this thesis and will be described below (28) , 

Suppose that after preliminary screening, it is concluded 


that the true model form is a member of the following set of m models, 
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which comprises the metamodel: 
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Y= O(% a) +e (2.23) 


where metamodel is defined as the mode] Space containing all the 
potential mechanisms which are considered Simultaneously, and 
X = vector of independent variables 


= vector of unknown parameters corresponding to the 
,th 


ay 
mode] 
y = dependent variable 


th 


€. = error term in the i mode] 


j 


The models may be in differential, rather than algebraic, form but 
that will not affect the argument to follow. 

Uncertainty about which model form is the true one is ex- 
pressed quantitatively in the form of a discrete probability distri- 
bution, here Pr is the probability that the qin model form is correct. 


Obviously, the set of PL's satisfy the axioms of probability, i.e. 
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esd area Ln P, = 1.0 (2.24) 


Conceptually, the Pris are personal or subjective probabilities in the 
Bayesian sense and not objective or frequency probabilities in the 
classical sense. Bayesian probabilities are a measure of our current 
State of knowledge about an uncertain quantity or hypothesis; they 
are modified with the acquisition of new information in accordance 
with Bayes' theorem. 

Before any data are taken, probabilities are assigned to 
the different model forms. In the absence of strong information, it 
is conventional to make these "prior Probabilities" al] equal. Ina 
Similar way, lacking specific information about the parameters, it is 
usually assumed that their values are uniformly distributed over the 
range - ~ to + ~ (called the improper uniform distribution). 

Given n sets of observations ty; > Reel c ele 2. 1s 9, tie 
posterior probabilities, Pie can be calculated in accordance with 


Bayes' Theorem. 


P. L(y|¢.) 


Pj" EPLUTS ies lis aaa) (2.25) 
J 


where Ps designates the probability prior to observing the data y and 
L(y|¢;) is the expected likelihood which is used to designate the pro- 
bability of observing the data, given that the model form os is correct. 


One disadvantage of the Bayesian discrimination procedure 
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A oe ca 
where ask =k 


element of the parameter vector aj that minimize 
the sum of errors squared for the has mode] 
C; re the diagonal element of the variance covariance 


matrix of the parameters, Cc; for the jth model 


In order to compute ci, the following equations are required: 


39 (x5 a.) ies 
Mea = ; 
=jik 0 as Te 
jk a a 
xX, = Xsiketl = Ly ¢#-u-f'h Ram al 2. .Y} (2.29) 
=n x r design matrix 
A. = xX." X, (2.30) 
eo ean, 


The c. can now be obtained by equation (2.31) 


Comoe Ad | (2.31) 


2.10.2 Design of Experiments 


Only in recent years has there been due attention paid to 
the planning of experiments that would be expected to result in maxi- 
mum discrimination among the rival models. The idea is to explore the 


entire operating space to locate a region in which predicted responses 
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4] 
of the various competing models, based on existing information, will 
provide maximum discrimination. The expected entropy change criterion 
used for design of experiments in this thesis is due to singh (28) and it 
is both exact and analytically simple unlike the earlier criteria that 


were only approximately correct (2°) , 


The underlying basis is the 
concept of entropy as a measure of the uncertainty associated with a 
Situation as outlined by Kulback (97), For a set of m models after ob- 


serving n data points, the entropy may be defined as 


m 
EB emernenh 13 Grp (2.32) 


th mode] being true after n 


where P. r is the probability of the i 
observations. A hypothetical (nt1)°¢ experiment is carried out and 
the expected posterior probabilities Ps nt] are now used to compute 


E It is the difference - R=E -E that is a measure of the 


n+] ° n n+] 
improvement in our knowledge. The following expression for - R gives 


the negative of expected entropy change 
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where of = error variance of a single measurement 
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t X. _) t 2 


PF nt = X0'Xi on Ayn) Xp o 


Xn = 1x r design vector for (nt1)S* row of the design 
matrix X 
n+] 
Yin = prediction from the jth gota correponding to Xp 
2 ae at 24 
Cnt] = ba Pian Yi yne] 


“A 


m 
Ynt1 ~ = dy Pian Yi nt] 


The hypothetical (nt1)S* experiment is carried out at various experi- 
mental conditions and the negative value of the expected entropy change 
is examined to obtain a region where - R is maximum. 

Computer programs used for the above calculations are given 


in Appendix B. 
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CHAPTER III 
REAGENTS AND CATALYSTS 


3.1 Feed Materials 

The nitrogen used in the kinetic study was of research grade 
and was supplied in cylinders by the Union Carbide Chemical Company. 
The purity specified by the manufacturer was 99.97% nitrogen. Mass 
spectral analysis on a typical sample showed that this gas contained 
no water or carbon oxides. 

The alcohols, 2-propanol and 2-butanol, were obtained from 
Fisher Scientific Co. Ltd. and were of spectroscopic grade. The 3- 
pentanol was obtained from Aldrich Chemical Co. Inc. and is also of re- 
search grade. The purity of each alcohol was assessed by gas chromato- 
graphy using a di-n-decyl phthalate-celite column. No contaminants were 
detected. Infrared spectra of these alcohols in the vapor phase showed 
only the characteristic peaks of the alcohols. No further purifications 
were attempted, except in adsorption studies where these alcohols were 
purified by repeated freezing and thawing under vacuum. 

Acetone, mixed with 2-propanol as feed was used in some of 
the experimental runs. The acetone, obtained from Fisher Scientific Co. 
(Catalogue No. A-20) was found to be pure by infrared and gas chromatographic 


analyses. It was used as received. 
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3.2 Alumina Catalysts 
5 yar | Preparation of Catalysts 


Only one type of alumina, known as Alon C, was used in this 
work and it was supplied by the Cabot Corporation, Boston, Massachusetts. 
The typical properties of this alumina, as specified by the manufacturer 
are listed in Table 3-1. Alon C, a trade-name of Cabot Corporation, 
is a fumed alumina consisting predominantly of the gamma form. It is 
made by the hydrolysis of aluminum chloride in a flame process. This 
Process produces a material of extremely smal] particle diameter, high 
Surface area and high purity. The y-alumina powder was used to prepare 
thin infrared-transparent catalyst wafers by pressing in a die under 
high pressure (2 ton/sq.cm.). 

The y-alumina doped with sodium hydroxide (Fisher Scientific 
Co. Catalogue No. 5-312) was prepared by adding the required amount of 
sodium hydroxide to the Alon C with a minimum amount of water to make 
a thick slurry. The water was removed by heating the slurry at 120°C 
in air and the resulting cake was crushed into fine powder for preparing 
catalyst wafers. In this way, catalysts with 1, 2, 3, 5, 8 and 15% by 
weight of sodium hydroxide in Alon C were obtained. Only one catalyst 
batch for each level of NaOH content was used in preparation of all 


catalyst samples examined in this thesis. 


Since the catalyst with 5% sodium hydroxide content had approxi- 
fe) 
mately 10 Na* ions per 100 (A), and since the maximum possible 
hydroxyl concentration on the alumina surface was estimated to be 13 OH 


0 . 
per 100 (A)? (99) » Slightly higher than 5% of sodium hydroxide 
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Table 3-1] 


Typical Properties of Alon C 


Note: 


Color and Form: 

X-ray Structure: 
Alumina Content*: 
Ignition Loss: 

Metallic Oxides**: 

Avg. Particle Diameter: 
Surface Area: 


pH (10% Aqueous Sus- 
pension): 


Specific Gravity: 
Loose Density: 
Bag Bulk Density: 


Refractive Index: 


White Powder 
90% Gamma form 
99% minimum 
4.5% maximum 
0.2% maximum 
0.03 micron 


100 m?/gram 


4.4 

3.6 

1.8 =e2a bs Curt t. 
3 - 4 lbs./cu. ft. 
1.70 


*Excludes physically and chemically combined 


water 
**Other than Al,0, 
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46 
content would correspond roughly to the fully doped level. 
3.2.2 Surface Area of Alumina Catalysts 
Surface areas of the alumina samples, after pressing into 
the wafers, were determined by the standard BET method (100) | Before 
each determination, the samples were evacuated for 8 hours at 350°C. 
Nitrogen was used as the adsorbate at the temperature of liquid nitrogen 


(- 196°C). Table 3-2 represents these results. 


Table 3-2 
Sodium Hydroxide Content and 


Surface Areas of Doped y-Alumina 


Catalyst Sample NaOH Content Surface Area 
gm/100 gm A1,0. (sq.m)/gm A10. 
A 0 9] 
B ] 75 
C 2 74 
D 3 55 
E 5 45 
F 8 39 
G 12 26 
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3.2.3 Pore Size Distribution for the Alumina Catalyst 
The pore size distribution was calculated according to the 
Procedure reported by Gregg and Sing! 101), The adsorption measure- 
ments were made with the catalyst wafer of pure y-alumina as adsorbent 
and nitrogen gas as adsorbate. The adsorption isotherm is plotted as 
shown in Figure 3-1, and the pore size distribution is shown in Table 


3-3 and is plotted in Figure 3-2. 


Table 3-3 


Pore Size Distribution for the Alumina Catalyst 


ry) ; 
Vv, cm Relative 

P/Py (Salers) rh Distribution 

ee ee ee ee OB ay se 
0.9751 86.50 398 0.00 
0.9727 79.62 364 0222 
0.9177 77.46 124 0.00 
0.8979 TiS ats: 100 0.22 
0.8794 ara tig! 85 0.90 
0.8754 51.92 83 5.39 
0.8428 26.42 66 2:03 
0.8181 18.68 58 leg 
0.7871 18203 49 0.82 
0.7499 9.52 42 0.48 
0.7106 7.40 36 0.26 
0.6079 5.66 27 O213 


ans 2 ee Eee eae ee 
Notes: p/p. = adsorbate gas phase pressure/satu- 


° ration vapor pressure 
Vv = volume of the adsorbate adsorbed 
r. = pore radius of the adsorbent 
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Figure 3-1. Adsorption Isotherm of Y-Alumina 
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CHAPTER IV 


APPARATUS AND EXPERIMENTAL PROCEDURES 


4.1 Study of Adsorption of Alcohols 


4.1.1 Infrared Cell 

A detailed diagram of the quartz reactor-infrared cell for 
the study of adsorbed alcohols is shown in Figure 4-1. The catalyst 
wafer was placed between two quartz rings (2.4 cm i.d.) and was moved 
in and out of the infrared beam by means of an external magnet. Two 
sodium chloride windows were sealed to both ends of the tube with a 
Silicone rubber potting compound (Fisher Scientific Co., Catalogue 
No. 4-769-5), which is 100% solid after curing and is stable up to 
250°C. The cell can hold a vacuum of 107° mm Hg and the potting com- 
pound is flexible enough to compensate for the different thermal ex- 
pansion coefficients of quartz and sodium chloride. To remove the 
NaCl windows for cleaning, if necessary, the cells were heated in an 
oven at 350°C for 4 hours and allowed to cool to room temperature slowly. 

This type of cell enables the catalyst sample to be heated 
in the heating section and moved to another section of the cell where 
infrared spectrum can be recorded. The heating section was made of a 
25 ohm heating coil whose outside surface was covered with a one inch 
thick asbestos insulation, so that the cell can be heated up to 650°C. 


The cell design did not facilitate a direct measurement of 
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Figure 4-1. Infrared Cell for Adsorption Studies. 
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52 
the interior temperature. However, a thermocouple was placed between 
the heater and the exterior quartz surface and the heating temperatures 
were recorded. A check of these temperatures against other values ob- 
tained with a thermocouple placed inside the cell under 10 torr of air 


indicated that the interior temperature was always substantially lower: 


Exterior temperature, °C Interior temperature, °C 


650 465 
500 370 
400 300 
300 220 
200 130 
150 100 
100 80 


Thus the catalyst temperature was calculated by the thermocouple 
placed outside the quartz surface and was then calibrated to the in- 
terior temperature by the above table. 
4.1.2 Vacuum System 

The vacuum system consisted of a mechanical roughing pump 
(the Welch Scientific Company, Skokie, T]llinois, Model No. 1405), a 
two-stage mercury diffusion pump, and a liquid nitrogen cold trap. The 
liquid nitrogen trap was installed to prevent oil and mercury particles 
from diffusing into the infrared cell. The vacuum in the cell was 
measured with a McLeod vacuum gauge and was about 5 x hee mm Hg under 


operation conditions. 
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4.1.3 Preparation and Pretreatment of Catalyst Wafer 

For the preparation of a catalyst wafer, about 100 mg. of the 
catalyst powder was pressed into a thin infrared-transparent wafer. 
Commercial dies are available and are equipped with facilities for 
pressing under vacuum. This process does not always lead to best 
samples. The author has achieved better results by not evacuating 
during the pressing process and by pressing the powder between two 
pieces of recording-chart paper (Perkin-Elmer Co., Part No. 221-1613). 

The die used in this study was fabricated in the machine 
shop, Department of Chemical and Petroleum Engineering, University of 
Alberta and was similar to the commercial dies (2.54 cm. dia.) except 
no evacuation line was provided. 

To prepare the catalyst wafer, the finely divided powder 
was pressed at two tons/sq. cm. between the two pieces of paper for 
10 seconds. On release of the pressure, the sample was easily re- 
moved fom the paper by a razor blade. The diameter of the wafer was 
2.54 cm and the thickness of the wafer was about 0.1 mm (roughly 
15 mg/sq. cm.). The catalyst wafer was then placed in situ in the 
infrared cell. 

The pretreatment of the catalyst involved degassing at 400°C 
for two hours, heating in oxygen at 10 cm Hg pressure for two hours at 
400°C, and again degassing at this temperature overnight. After cooling 


to room temperature, the catalyst was ready for the adsorption study. 
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4.1.4 Adsorption, Desorption and Recording of Spectra 

The catalyst wafer, after pretreatment, was placed in the 
sample beam of a Perkin-Elmer Model 621 Infrared spectrophotometer 
installed with a Perkin-Elmer Model 21 presample chopper. The base- 
line spectrum of the catalyst was recorded using an identical evacuated 
cell without a catalyst wafer in the reference beam. Because of 
scattering of the infrared beam by the wafer, especially at the high 
frequency end of the infrared spectrum (3000 to 4000 cm !), the trans- 
mission may be considerably reduced to anywhere between 5 to 40 percent. 
In order to obtain a good spectrum, the reference beam had to be at- 
tenuated and the slit opened to compensate for the loss in energy. 

The slit opening was thus set directly proportional to the square root 
of the reciprocal of the original transmittance. For the required 
accuracy of the spectrum, slower scanning speeds were also necessary. 
Detailed operating instructions are available in the Perkin-Elmer 
Model 621 Manuals (Perkin-Elmer Corporation, Conneticut, U.S.A. 
Catalogue No. 990-9479). 

After purification by repeated freezing and thawing on the 
vacuum rack, the alcohol was introduced at about one cm Hg absolute 
pressure into the infrared (sample) cell. To obtain a spectrum due 
solely to the combination of catalyst and added adsorbate, it is necessary 
to correct for any spectral adsorption due to the gas phase present 
in the cell. Often the cell pressures are low enough (0.5 to 2.0 cm 
Hg) that the gas spectrum is negligible, but, if it is appreciable, 


its absorption effects must be eliminated. This is most conveniently 
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done by using an identical cell with no catalyst wafer but containing 
the alcohol at the same pressure as in the sample cell in the reference 
beam of the spectrophotometer. The y-alumina shows strong infrared 
absorption at frequencies below 1000 cm”! and this forms a lower limit 
to the useful frequency range for such adsorption studies. 

After the spectrum of the adsorbed molecules at room temperature 
was recorded, the cell was then either pumped off or heated to various 
temperature levels. Depending on the sequence of experiments being 
performed, various combinations of heating and/or degassing cycles 
were used. The infrared spectra were always measured at room tempera- 
ture during the experimental cycles. To evaluate the extent of re- 
action, the gas phase was analyzed by infrared and gas chromatographic 


Or mass spectral analyses. 


4.2 Simultaneous Studies of Kinetics and Mechanism 
4.2.1 Equipment 
Simultaneous studies of the reaction rate and the infrared 
spectra of adsorbed species for the decomposition of an alcohol were 
performed in an infrared cell which was designed to be operated as a 
recycle reactor. The main advantages of using a recycle reactor are, 
1. near-isothermal conditions exist throughout the catalyst 
bed and, 
2. high recirculation rates eliminate diffusional effects. 
The equipment used in the study of dehydration and dehydro- 


genation of 2-propanol consisted of an alcohol feed system, a recircu- 
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56 
lation pump, a vaporization system, a reactor (IR cell), the temperature 
control and recording system, a product collection system and a gas 
spectral analysis system. A schematic diagram of the apparatus is 
shown in Figure 4-2. The recirculation line was heated to prevent the 
condensation of the vapor products in the line. Heating tapes and 
voltage controllers (Variacs) were used to maintain the line temperature 
at 170°C. In the following sections, each system will be discussed in 
detail. 

(1) Feed System 

Both gas and liquid feeds were introduced into the vaporiza- 
tion system simultaneously (see Figure 4-2). The gas feed used in this 
work was nitrogen and it was supplied from a pressurized nitrogen 
cylinder. The pressure of the nitrogen supply at the cylinder was 
determined by a gas pressure regulator, and the flow-rate of nitrogen 
was controlled with a needle valve. The flow-rate of nitrogen was 
measured with a soap bubble meter. A Matheson 601 stainless steel 
rotameter was connected to ensure a steady nitrogen supply during an 
experimental run. 

A Sage Model 355 syringe pump (Sage Instruments, Inc., White 
Plains, New York) was employed to inject liquid feed, 2-propanol, into 
the vaporization system at preset speeds. 

‘The syringe pump was designed to infuse fluids at continuous ly 
variable flow rates with a ratio 25,000 to 1 at the maximum and minimum 
settings with any single syringe. The pump motor drives the gears on 


top of the pump which engage mating racks on the drive carriage moving 
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58 
it uniformly forward to drive the syringe piston smoothly. An electronic 
feedback circuit regulates the permanent magnetic DC motor to a repro- 
ducibility of + 0.3% full scale, independent of changes in back pressure 
or + 10% variations in line voltage. 

Liquid 2-propanol was fed into the system by a Hamilton 
90 cc gas-tight syringe (Catalogue No. 1050) which features a teflon- 
coated piston in the glass barrel. Very precise and constant feed 
rates with average variations of less than 0.4% were obtained with 
this syringe. The calibration of the feeder is given in Appendix C. 

(2) Vaporization System 

The vaporizer was made from 2 cm 0.D. 22 cm. long Pyrex glass 
tubing. At the inlet, a 1/8 in. Kovar joint was connected to the 
syringe by 3/32 in. 0.D. stainless stee] tubing. Another inlet tubing, 
1/8 in. 0.D., was used to introduce nitrogen gas. The outlet of the 
vaporizer was a 12/5 female ball joint socket which was coupled to the 
inlet of the infrared cell (reactor). A 3/32 in. iron-constantan 
thermocouple was placed inside a 2 mm I.D. glass thermocouple wel] 
near the outlet of the vaporizer. This thermocouple was connected 
to the Honeywell temperature controller (SCR trigger controller, Model 
R7161J) which controlled the vapor feed temperature at 170°C. The out- 
Side of the vaporizer was covered with 1-1/2 cm thickness of asbestos 
insulation. 

(3) Recirculation Pump 
A magnetic solenoid pump was incorporated into the recircu- 


lation loop. Efficient circulation of the product stream was important 
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59 
Since it would ensure that no concentration gradient existed within 
the gas phase and that no mass transfer effects influenced the measured 
reaction rate. 

The double acting solenoid pump was of stainless steel] con- 
struction with teflon-coated piston, valves and O-ring seals. A 
sectioned sketch of one-half of the pump is shown in Figure 4-3. The 
two halves of the pump are identical. The O-ring connector on the 
pump discharge was omitted from the drawing. 

Both suction and discharge valves were formed from teflon 
disks which were back-seated with stainless steel Springs. The four 
external connectors were "Swagelok" 1/4 inch straight-thread connectors. 
The suction valve disks were seated against the threaded end of the 
connectors. These connectors had been previously machined to provide 
the valve seat. The discharge valve seats were machined in the pump 
heads. The check valves appeared to be the most critical components 
in the design and operation of the pump. 

The pump cylinder was made of one inch I.D. stainless steel 
tubing. Its total length was 13-3/4 inches. 

The pump piston was constructed of mild steel coated with 
teflon. The clearance between the piston and cylinder wall could be 
varied by adjusting the thickness of the teflon coating. The overall 
piston length was four inches. 

Three coils, 3-2/3 inches long and separated by 1/4 inch fiber 
disks, were wrapped on the pump cylinder. Each coil had approximately 


2300 turns of No. 20 enameled copper wire. A variable voltage (9 to 
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6] 
18 volts D.C.) power Supply with a maximum current of three amperes was 
used to excite the solenoids. At the design condition of 12 volts 
D.C., each coil carried a current of 1.1 amperes resulting in a rating 
of 2530 ampere-turns. 

The external electronics consisted of the D.C. power supply, 

a repeat cycle multi-gang timer and three Switching relays. The timer 
was a three gang, single cam assembly driven by a 60 r.p.m. synchronous 
motor installed with 8-variable speed drives. The three gangs were 

timed so that they were sequentially offset by 90 degrees. The pump 
piston would therefore travel one complete cycle for each revolution 

of the cam. Figure 4-4 shows a wiring diagram of the timer and switching 
relays. 

The designed capacity of the magnetic solenoid pump was ranged 
from 15.6 to 65.5 liters per minute, but the actual capacity was approxi- 
mately three quarters of that expected. This reduction was mainly 
attributed to the efficiency of the check valves. 

The pump was insulated with 1/2 inch of asbesto tapes so that 
the steady-state operation temperature was maintained at 170°C. The 
pump was mounted horizontally between the two infrared cells. This 
pump could be operated at a temperature of 170 + 35°C without readjusting 
the thickness of teflon coating on the piston. 

(4) Infrared Cells and Sample Holder 

The pyrex reactor-infrared cell is sketched in Figure 4-5, 

The cell length was 11 cm. A 40 mm diameter NaCl window was sealed 


onto one side of the cell with the silicon rubber potting compound. 
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Another side of the cell was sealed by a 50 mm diameter NaCl window. 
Pyrex vacuum stopcocks were used to isolate the reactor in evacuative 
treatments of the catalyst. The gas inlet and outlet to the cell are 
indicated in Figure 4-5. An O-ring (Viton, 40 mm diameter) joint 
allowed the removal of the front end of the cell for the replacement 
of the catalyst samples. The catalyst wafer, prepared according to 
the procedure described in Section 4.1.3, was supported in the cell by 
the sample holder (see Figure 4-5) and was facing the NaCl windows. 

The central part of the cell was heated by a 20-ohm heating 
coil placed outside the cell while the windows were kept at temperature 
about 170°C by blowing cold air around the ends of the cell. The out- 
Side of the heater was insulated by one inch thickness of asbesto tapes. 
The temperature in the reactor cell was controlled by a variable volt- 
age transformer (Variac, Fisher Scientific Company). The voltages applied 
to the heating coil were measured by a voltmeter and were controlled 
within + 0.1 volt of the set point. In this way, the reactor temperature 
was controlled with a variation of less than 0.3°C. 

A measure of the sample temperature was provided by a 1/8 
inch thermocouple introduced into the cel] via a silicon rubber septum; 
the configuration was such that the junction formed a pressure contact 
with the catalyst wafer. It can be said that the temperature thus 
measured was very close to the true wafer temperature based on the 
following observations: 

(a) The measured wafer temperature was found to decrease by less 


than 0.2°C when the cell was evacuated at 400°C. 
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(b) Under the evacuated conditions and at a temperature of 400°C, 
a 1/50 inch thermocouple gave temperature readings which were identical 
to those measured by the 1/8 inch thermocouple. 

One side of the cell was connected to the magnetic solenoid 
pump using a Kovar joint (1/4 inch) and the other side to the reference 
cell using glass tubing. The reactor cell can be used up to a temperature 
of 450°C and a vacuum of 5 x 107? mm Hg. 

The reference cell was made of pyrex glass tubing, 27 mm 1.D. 
and 11 cm long. Both ends of the cell were sealed with NaCl windows 
(40 mm diameter). The cell was wrapped with a heating tape and thus 
the temperature in the cell was maintained at 170°C by applying 34 
volts on the heating tape. The cell was insulated with asbesto tapes 
to prevent excessive heat loss to the Surroundings. 

The reference cell was used to measure the partial pressures 
of both the reactant and the products. This provided the information 
for calculating the material balance as well as the conversion. 

(5) Temperature Recording System 

The thermocouple (1/8 inch iron-constantan) located in the 
reactor was connected to a Speedomax Type G temperature recorder 
(Leeds and Northrup Co.). The thermocouple was calibrated before in- 
stallation. The temperatures recorded in the calibrated isothermal 
bath between 273 and 673°K versus the voltages measured from the thermo- 
couple can be represented after applying a least squares fit (see 
Appendix C for computer programs), by y = 5.5112 + 18.0013x 
where y = temperature of the bath, °K, 


and x = signal measured, mv. 
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The temperature recorder was checked occasionally by feeding 
a known voltage from a potentiometer into the recorder. Thus the 
temperature on the catalyst wafer could be recorded Precisely during 
experimental runs. 

(6) Product Collection System 

The vapor products leaving the recirculation loop were in- 
troduced through a fluid control valve (Hoke, Inc., Catalogue No. 
4172G2Y) into the condenser-sampler immersed in an ice-water bath. 
For better liquid sample collection, a pyrex glass integral condenser- 
sampler was used (Figure 4-2). The noncondensable gases from the 
liquid sampler were then passed to a 75 cc gas sampling burette and 
then to a soap-bubble meter. The samples thus collected were sent 
for gas-phase chromatographic analysis. The results were used as an 
additional check of the number of compounds appearing in the reaction 
product mixture. 
4.2.2 Operation of Equipment 

The procedure used in preparing the Catalyst wafer was identical 
to that described in Section 4.1.3. The catalyst wafer was placed in 
the sample holder which was then inserted into the center part of the 
infrared cell in such a way that the catalyst wafer formed a pressure 
contact with the thermocuuple. This insured that the catalyst wafer 
remained at the same position for all runs. After the O-ring connection 
was tightened (see Figure 4-5) the infrared cell was fixed in position 
to the recirculation loop through the two glass ball joints (12/5). 


After all joints were sealed with silicone grease and secured, the 
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equipment was ready for an experimental run. 

The vacuum system was similar to that described in Section 
4.1.2 except that the McLeod vacuum guage was replaced by a Pirani 
vacuum gage (Edwards High Vacuum Limited, Model G9). The pretreat- 
ment of the catalyst was performed by evacuating the sample cell at 
a catalyst temperature of 350°C for 12 hours. The base line of the 
catalyst wafer showed an identical spectrum to that obtained in the 
adsorption study (Section 4.1). However, it was learned later in the 
research program that heating the catalyst wafer by continuously flush- 
ing nitrogen gas through the sample cell for 12 hours gave the same 
results as pretreating the catalyst in vacuo. Therefore, except for 
the first catalyst (catalyst (a) in Figure 5-10), all other catalysts 
were pretreated according to the latter method. It should be mentioned 
that during the catalyst pretreatment, the magnetic solenoid pump, the 
recirculation loop heaters and the vaporizer heater were turned on 
since they took four hours to reach a steady-state operation condition. 
Thus, except the alcohol feeder, theentire system was under working 
conditions when the pretreatment of the catalyst was finished. 

At this point the nitrogen gas flow rate was measured by the 
bubble-meter. It was found that the nitrogen flow rate could be main- 
tained constant for all the experimental runs provided the pressure in 
the nitrogen supply cylinder exceeded 300 psig. Thus no adjustment 
of the metering valve was necessary. 

The base line spectrum of the catalyst was taken by an on-line 


computer data acquisition system involving coupling of the spectro- 
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photometer to an IBM 1800 digital computer (see Appendix A). 

After obtaining the base line spectrum of the catalyst, the 
alcohol feeder was turned on and the reactant was passed through the 
vaporizer into the reaction system. The reactant feed rate was ad- 
justed by the syringe pump control dials (see instruction manual of 
Syringe pump, Sage Instrument, Inc.). The partial pressure of the 
reactant in the recirculation loop was controlled by the flow rate of 
the reactant. The reference cell was then inserted into the sample 
beam of the infrared spectrophotometer. The frequency, at which the 
desired 2-propanol absorption band occurred, was set and peak intensity 
was recorded as a function of time. This can be done simply by turning 
the scan clutch control of the spectrophotometer to the “out" position 
at the frequency of interest. After the wafer temperature and product 
mixture composition in the recirculation loop became constant, the on- 
line computer data acquisition system was activated and steady-state 
was assumed when the peak intensity became constant for 20 minutes. A 
data point for the kinetic measurement was then taken. This consisted 
of 

(a) recording the absorption intensity of 2-propanol and propylene 
in case of dehydration studies, or of 2-propanol and acetone in case of 
dehydrogenation. 

(b) inserting the sample cell in the sample beam and the reference 
cell in the reference beam and recording the spectrum of the adsorbed 
species (plus base line), 

(c) analyzing the liquid and gas samples collected in the samplers 
by gas chromatograph. (This provided a check for the presence of any 


side reactions). 
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(d) recording the catalyst temperature and alcohol feed rate. 

(e) plotting the spectrum, after eliminating base line spectrum 
with the computer, on the IBM 1627 plotter. (see Appendix A for the 
plotting program). 

About two hours were required per data point. When the 
temperature was varied, about 30 minutes elapsed before a new s teady 
State condition was obtained. For feed rate changes about 20 minutes 
were required to reach a new steady state. 

4.2.3 Data Evaluation 
(1) Quantitative Analysis of Product Mixture 
Quantitative infrared analysis is similar to routine colori- 


metric analysis and makes use of the Beer-Lambert law: 


I 
: OA) = £2P 
where A, is the absorbance at wavelength A, Ip, and I , the 


intensities of the incident and transmitted beams 
respectively, | 

C the concentration 

& the length of the cell 

€ the extinction coefficient 

P partial pressure of the gas sample 

T absolute temperature and 

R the gas constant. 


It is important to test the validity of this law in the range of con- 
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70 
centrations likely to be employed. 
In the case of gaseous mixtures, if two or more components 
of the mixture absorb at the same wavelength, the total absorbance 
is assumed to be equal to the sum of the absorbances due to each one 


of the components in the mixture, i.e. 


i ms 
A, ; baa C. pF | (Ro (4.2) 


where the subscript, i, refers to the particular species. However, 
Overlapping absorption bands tend to decrease the accuracy of the 
analysis. Several methods such as the ratio method, etc., have been 
developed to overcome the difficulties from the overlapping of bands (102) | 
Fortunately, in this study, it was possible to choose an intense iso- 
lated absorption band for each component in the product mixture. Thus, 


Equation (4.2) reduces to 


Aj) = et i C. Lo = (RT) (4.3) 


Thus frequencies chosen for each component in the reaction mixtures 
are listed in Table 4-1. 

Because the vapor phase spectrum of water was found to be 
inaccurate for the quantitative analysis, and hydrogen gas was found 
to be infrared-inactive, the partial pressures of water from dehydra- 
tion and of hydrogen from dehydrogenation were taken to be stoichio- 


metrically equal to that of propylene and acetone, respectively. 
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Table 4-1 


Frequencies Employed for Quantitative Analysis 


a eee ee 
Applicable partial Frequency, 


Component pressure range, ie =i 
cm Hg 
(a) Dehydration 2-propano] ee 10 1250 
2-propanol nig Ke) 3660 
propylene <a 3106 
(b) Dehydrogenation —  2-propanol < 6 1148 
2-propanol 26 3660 
acetone <Ono 1740 


Since no evidence for the occurrence of reactions other than dehydra- 
tion or dehydrogenation was encountered, this approach was not really 
questionable. 

In infrared quantitative analysis the accuracy in recording 
the peak intensity is very important. In practice, the peak height 
or the area of the band is taken as the criterion of band intensity. 
The former method was chosen in the present study since it was simpler 
and it gave a good straight line correlation according to Beer's law. 
The accuracy of the evaluation of I/T was also found to be important 
in preparing a good calibration curve therefore, the error of the I/Ip 


obtained was estimated by the following procedure. 
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The fractional infrared transmission, I/T), is the ratio of 
the incident infrared radiation transmitted by the gas to that trans- 
mitted in the absence of the absorbing gas. In the present study Ih 
was always greater than 80% of the total transmission. The value of 
"I" decreased, depending on the concentration of the gas in the reference 
cell, sometimes to values as low as 5% of the total transmission. The 


deviation in reading the recorder pen was 0.3%. The error in T/T), is 


given by, 
ADs (F true % eee 
eb HE 
E = oo si ps) t Soa) (4.4) 


Since Ty was always greater than 80%, a value of 80 will be used to 


determine the maximum error. 


Pa 
Emax > (go) > (gpsacg) + (ag S a3) 


m 
i] 
+ 


+ .000047.I + 0.00375 


m™m 
" 
+ 


fa max ~ 7 *9065%, for maximum I = 80% . 


Thus, the error introduced due to misreading the recording chart is 
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negligible. This accuracy was improved even further by reading the 
results with an on-line computer which read the transmission data 
accurate to 0.1% (see Appendix A). 

The calibration curves were prepared by introducing different 
amounts of vapor into the reference cell. The cell was connected to 
a quartz precision pressure gage (Texas Instrument Model 141A) and 
Was evacuated to obtain the zero pressure reading of the gage. After 
introducing the vapor into the cell, the vapor was found to be adsorbed 


(103) that this difficulty 


partly on the glass wall. It was shown 
could be minimized by prior conditioning of the system with the vapor 
sample. However, the author found that allowing a thirty-minute 
stabilization period after introducing the vapor was a more convenient 
way to eliminate the above difficulty. It was found that for all the 
absorption frequencies used in determining the concentration no mutual 
interference of the absorption band was detected, e.g. the band in- 


tensity for 2-propanol at 1250 cm”! 


was not affected by introducing 
different quantities of propylene or water. The absorbance is plotted 
against partial pressure to obtain the working curve at 29°C. (After 
exposing the cell to the infrared beam, the temperature of the cell was 
recorded to be 29°C which was 5°C above the room temperature). For 
higher pressures of 2-propanol (10 < P < 60 cm Hg), the calibration 

was performed by connecting the reference cell to a mercury manometer 
and the entire system was heated at 100°C to ensure no condensation 


Of 2-propanol vapor. The vapor pressure of 2-propanol as indicated 


by the manometer was read with a cathetometer. The absorption band at 
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3660 cm! was recorded by compressing the band intensity (% transmittance) 


by a factor of ya), 


This provides greater accuracy for the band 
intensity measurements. Figures 4-6 to 4-9 and Tables 4-2 to 4-6 show 
the calibration results after fitting by least squares. 
The pressure-broadening effect was found to be unimportant 
Since the absorbance of the gas samples (2-propanol, acetone, and 
propylene) was not changed by increasing the cell pressure to one 
atmosphere by injecting nitrogen gas. 
(2) Evaluation of Kinetic Data 
In the recirculation loop, the recycle rate was found to be 
32.8 liters per minute. The volume of the entire loop was estimated 
to be 0.72 liter. Consequently, the conversion per pass was low, and 
the loop was treated mathematically as being equivalent to a well- 
stirred flow reactor. The reaction rate was thus calculated from 
the following material balance on component qill2), 
(Sapa 10 “A (4.5) 
In order to obtain (-ry)> the fractional conversion of the reactant, 
Xp» should be determined. This was done by determining the partial 


pressures of the reactant in the feed, Po» and in the product stream 


Ps? and then calculating Xn from 


yee PB. (4.6) 






eae sz 
v 4D 







to oteraviIoo [snotsant? 


[stt1sq nig pitinetatsb yd snob eaw etril -saaeieetata 24 blue, ofte.. 


awente er bie 
re 


" 
ne 


ofts(vortog1eds al. © yiy 
ne sit to. smuToy odT csdunte 19g 21stth.B 
219Vi it ,¥itasuesenod...vettl St D ad 
. | 

ontod. 26 {697 I sasntem. bedsexs,26ew gool 


ie 



















Re i 
\ aft piv 2291QMNOS vd bebyoss? Reslancait iE 
ee ea ee nosast.B.yd 


U-h OF +f ini h 


. et nomayu2som etaoegel 
(a Dp AFIT? sedts. atfueer nofisydiisa ond 
4a ontaebsord-syvaze yg ofT 5) 


1-3) 2a hqmna 26p it. dn sonsdioeds., odd ¢ 


Onteasionr vd bapa aro Jon 260. (onet ‘a o* 
26g nspottin pritostnt yd avenge 


bid ottantr. to ngtssulev3. a 


26w S367 notiossy sil. .»velnssy wort beta 


ABHognos -110 


sons isd [etr9sem  prtwot Totem 


¥ a ' 


ag oly ~) nistdo of 9b at 







AB 994 , bast oH ig vi anstano% —~ were 1c 
ait “nor * ae ; afug 6) ns nt brs | 


a 
- 
™ 


< i 


ABSORBANCE 


1.6 


1.0 


9 
co 


4 
a 


0.4 


0.2 


Figure 4-6. 


» at 1148 cm | 


e at 1250 cm™ 
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PARTIAL PRESSURE OF 2-PROPANOL, cmHg 


Calibration of 2-Propanol Absorbance Using 
Beer's Law. 
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Figure 4-7. 


@ AT 3660 cm ~ 


10 20 30 40 
PARTIAL PRESSURE OF 2-PROPANOL, cm Hg. 


Calibration of 2-Propanol Absorbance Using 
Beer's Law. 


76 










tm OBOE TA @ 


Pale em, ha: p a i eee en are 
: i, . . : y Sa " | 2 AY mm 
a a ts Tak Bi a | a ik cok: de , 
. i x 
: 4 | : 
» SY RTE Gh fa Ave, ake 


ABSORBANCE 


Figure 4-8. 


aK 


@® AT 3102 cm 


4 6 8 10 v2 


/ 


PARTIAL PRESSURE OF PROPYLENE, cm Hg. 


Calibration of Propylene Absorbance Using 
Beer's Law. 
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PARTIAL PRESSURE OF ACETONE, cm Hg. 


Figure 4-9. Calibration of Acetone Absorbance Using 
Beer's Law. 
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TABLE 4=2¢ CALIBRATION OF 2=PROPANOL ABSORBANCE 


USING BEER'S LAW (1148 CM=1) 


X = PARTIAL PRESSURE OF 2=PROPANOL»s CM HG 


Y = ABSORBANCE 


THE COEFFICIENTS OF THE POLYNOMIAL» Y = AO + Al * X» 


AO = #0200154 


Al 


0229960 


REGENERATED DATA 


X MEASURED Y OBSERVED Y CALCULATED PCT ERROR 


10287 0385 00384 00292 

20285 0680 02683 0s401 

30053 020915 0«913 00206 

4502 1346 16347 02045 

50443 1629 1629 02019 
VARIANCE = 02000003 


STANDARD DEVIATION = 04001787 


MAXIMUM PCT ERROR = 00401405 
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TABLE 4=3¢ CALIBRATION OF 2=PROPANOL ABSORBANCE 


USING BEER'S LAW (1250 CM=1) 


PARTIAL PRESSURE OF 2=PROPANOLs CM HG 
ABSORBANCE 


x< 
Hou 


THE COEFFICIENTS OF THE POLYNOMIAL» Y = AO + Al * Xo 
AO = 0200270 


Al = 0015966 


REGENERATED DATA 


X MEASURED Y OBSERVED Y CALCULATED PCT ERROR 


02824 02131 00134 12793 
12304 0209 O«210 00526 
12658 00265 O«267 00876 
26300 00373 02369 02930 
30146 00510 02505 00999 
4«550 0730 0.729 00113 
50474 00876 0«876 02068 
60756 12077 12081 00379 
8«155 16305 16304 02019 
10.031 16605 16604 02051 
VARIANCE = 02000007 


STANDARD DEVIATION = 06002763 
MAXIMUM PCT ERROR = 14793472 
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(RT E60 £8061 


21000 
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TABLE 4=4¢ CALIBRATION OF 2=PROPANOL ABSORBANCE 


USING BEER'S LAW (3660 CM=1) 


x 
w 


PARTIAL PRESSURE OF 2=PROPANOL»s CM HG 


=¢ 
u 


ABSORBANCE 


THE COEFFICIENTS OF THE POLYNOMIAL» Y = AO + Al * Xo 


AO = =0200749 


Al 0204846 


REGENERATED DATA 


X MEASURED Y OBSERVED Y CALCULATED PCT ERROR 


50307 00252 02249 02905 
82894 00421 00423 00514 
166682 0800 O+«80l O«ol27 
276986 1350 12348 020104 
386«823 1.873 16874 02048 
49e211 20378 20377 00017 
VARTANCE = 02000002 


STANDARD DEVIATION 


02001670 


MAXIMUM PCT ERROR = 02905624 
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TABLE 450 CALIBRATION OF PROPYLENE ABSORBANCE USING 


BEER'S LAW (3102 CMe1) 


~< 
u 


PARTIAL PRESSURE OF PROPYLENE CM HG 


-~< 
" 


ABSORBANCE 


THE COEFFICIENTS OF THE POLYNOMIAL» Y = AO + Al * Xo 


AQ 


0002012 


Al 0095658 


REGENERATED DATA 


X MEASURED Y OBSERVED Y CALCULATED PCT ERROR 


20521 20483 20431 20067 

46989 40835 4.792 02878 

76514 70011 70207 22808 

102009 90641 90594 00481 

126152 11707 112644 00533 

146725 142100 142105 02041 
VARTANCE = 02009859 


STANDARD DEVIATION = 02099296 
MAXIMUM PCT ERROR = 26808381 
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TABLE 4"6¢ CALIBRATION OF ACETONE ABSORBANCE USING BEER'S 
LAW (1740 CM=1) 


PARTIAL PRESSURE OF ACETONE CM HG 
ABSORBANCE 


< 
nh ft 


THE COEFFICIENTS OF THE POLYNOMIAL» Y = AO + Al * Xo 


AO = =0«00032 
= 0+66321 


REGENERATED DATA 


X MEASURED Y OBSERVED Y CALCULATED PCT ERROR 


O«151 02099 02099 02e052 
06283 0«187 O«187 00140 
00479 0317 00317 Oelll 
0«682 0453 00451 02268 
12001 0663 02663 0«067 
12196 00793 0«792 0.091 
1.357 0+898 02899 00139 
1.568 12039 1.039 02008 
16808 1.6198 1.198 O«011 
22100 12392 16392 02020 

VARIANCE = 02000000 

STANDARD DEVIATION = 02000672 

MAXIMUM PCT ERROR = 0268247 
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Before starting an experimental run, a sodium chloride wafer, 
of the same size as the catalyst wafer, was placed in the sample holder 
with the equipment operating under steady-state conditions. Since 
NaCl was found not to catalyze reactions of alcohol, it was used to 
determine the partial pressure of the reactant in the reference cell 
when no reaction occurs. The partial pressure thus obtained was as- 
sumed to be identical to Po: 

The partial pressure of 2-propanol was measured according to 
the procedure outlined in the previous section [4.2.3 - (1)] and was 


then calculated using equation (4.7): 


: 
é one 973%. 170 
Po ~ Pcalibrated * ee Pcalibrated * Bri (4.7) 


where Peal ihrated is the partial pressure of 2-propanol obtained 
from the calibration curves (Figures 4-6 and 4-7). 
To is the temperature of the reference cell, 170°C, and 
Ty is the calibration temperature, 29°C. 
The partial pressure of 2-propanol in the reference cell was found to 
be proportional to the feed rate by a least Squares fit (see Appendix C 


for computer programs), which can be given by 
Py = 0.09518 + 92.9448 (Fao) (4.8) 


The results are plotted in Figure 4-10 and are presented in Table 4-7. 


The validity of the calibration was checked before and after the experi- 
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Figure 4-10. Relation of 2-Propanol Feed Rate to Its 
Partial Pressure Under Non-reaction 
Conditions. 
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TABLE 4=7e RELATION OF 2=PROPANOL FEED TO ITS PARTIAL 
PRESSURE UNDER NON=REACTION CONDITIONS* 


2=PROPANOL FEED RATEs GRAM=MOLE PER HOUR 
PARTIAL PRESSURE OF 2=PROPANOL IN REACTORs CM HG 


<< a 
oon 


THE COEFFICIENTS OF THE POLYNOMIAL® Y = AO + Al * Xs 


AQ = 0209518 
Al = 92094480 


REGENERATED DATA 


X MEASURED Y OBSERVED Y CALCULATED PCT ERROR 
02021 2e051 22051 02032 
02052 42940 40938 02029 
02052 92591 92594 02032 
00170 150927 150932 00037 
00245 222903 222894 0+036 
00329 30733 300748 02050 
00426 390723 392698 02060 
00496 46290 466260 06062 
02584 540343 540412 Ool27 
02640 590611 59579 02052 

VARIANCE = 02000836 

STANDARD DEVIATION = 020028925 

MAXIMUM PCT ERROR = 006127182 


*NITROGEN FLOW RATE = 00378 G=MOLE/HRe» 
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mental runs and was found to hold quite well. 

After obtaining Pos the sodium chloride wafer was replaced 
by a catalyst wafer and the partial pressure of each component in the 
reaction mixture was collected by measuring the absorbance at the 
characteristic frequencies for each component after steady-state was 
reached. Though the reaction mixture was always more than 50% nitrogen 
and the percentage conversion of the reactant was low, the change in 
molal flow rate was calculated using the procedure described below. 

The partial pressure of each component in the product mixture 


and fractional conversion of 2-propanol were calculated as follows: 


(a Fag + Xn) + b) 


Bae Das a bantce a (4.9) 


p (aiFen( lee x, )o+ ib) 

and ‘ee Xn = 2 = Roy ° —A0 _ A (4.10) 
0 (a Fag + b) 

where P, and Pp, are partial pressures of the jth component and of 


alcohol, respectively, corrected for the volume expansion 
due to reaction and 


th component and of 


Pip and Poy are partial pressures of i 
alcohol, respectively measured in the reference cell. 
The term, (a Fag (] + Xn) + b)/(a Fag + b), is the correction factor for 
the change in molal flow rate due to chemical reaction. Thus by com- 
bining equations (4.10) and (4.5), the reaction rate can be calculated. 


By changing the feed rate, i.e. the partial pressure of the 
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reactant, a set of data in the form of (- ra) = F(p.) could be obtained 


for modeling the kinetics. 
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CHAPTER V 
RESULTS AND DISCUSSION 


5.1 Adsorption and Surface Reactions of Alcohols 
ee ea 2-Propanol on Pure y-alumina 


Figure 5-1 shows the spectra of 2-propanol adsorbed on pure 
y-alumina from room temperature to 300°C. Up to 80°C (Curve 2), one 
observes that the surface hydroxy] groups are displaced considerably 
to the lower frequency with the formation of a large broad hydrogen- 
bonded hydroxyl region. On adsorption, the three bands observed in 
the gas phase at 2981, 2972 and 2888 cm”! have shifted to the lower 
frequencies at 2960, 2930 and 2870 cm), The C-H bending region shows 
almost no change in frequencies except that the intensity of the band 
at 1465 cm”! has increased two-fold over that in the gas phase in con- 
trast to the band at 1373 cm |, These changes are similar to that 
observed during the transition from gas to liquid phases (94) | The 
three distinct bands at 1245, 1145 and 1065 cm! in the 0-H bending 
and C-0 stretching regions, observed in the gas phase spectrum, on 
adsorption merge into a single broad band at 1170 cm, These ob- 
servations suggest that 2-propanol is physically adsorbed through the 
hydrogen-bonding between the -C-OH of the alcohol and the surface 
hydroxyl groups (as shown tn Structure V). 


The close similarity between the spectrum of the adsorbed 


89 
















v AAT BAND 


wolz2ua2i0. duis 2T.ue3A 
2forgolA 10 2tobtoseh sosty 9 | 
Sf ny | ud f tha 


e7uq mo bedvo2zbs Tonsqo1g-$ ¥0 sitoage oat ewode f-2 aS : 
0 «(S$ avid) 3°08 of qU .9°OOE ot Swwtsreqgnet moor moi? sntmuls 
\idersbrenos beoslqeth sys equedp tyxovbyrl soetiwe srt tans sovised 
-nsporbyd bsovd ‘opt. B td nottemot eft dttw yponaupart ewol a8 9 
nt bayisedo. ebned serit ont ,noriqroebs 0 ‘notes Tyxorbyl b 
newol sit ot battrda2 vei” mo 8889S bas Stes [82S 36 sesdq asp § 
ewone nofpsy entbned H-3 aA ao 0888 brig oees , 08S 36 29 


7 has 
4 


brisd oft to ySienstnt edd teddy tqe0x9 estonsupert mt spned> on f if 





-n02 nf seatq 260 shi nf sends 19v0 blot-owd beesorant esd ar) daar ts 
Jerid of vs/imi2 s1s 2egmmdts seedT -./"mo CNET 28 basd Sit of tenn 
AS*) sa zea brupil of esp mort not tenes eft patwh . 
pnibasd H-0 ait nt Pino ea0t bns aert (BhST 36 2bned ens er 
f10 msa0q8 ae ‘e6p” ont nt bev1sedo , notees pete 


-do aeodT ino ba ii an _ dani oa 6 linia ne 


any 








TRANSMITTANCE (%) 





3800 


Figure 5-1. 
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Spectra on the Adsorption of 2-Propanol on 


y-Alumina. 1. Baseline. 2. Alumina and 
2-Propanol at 10 torr and room temperature. 
3. Alumina and 2-propanol at 300°C for 1 hr. 
4. Alumina and water vapor at 5 torr and 
room temperature. 5. Degassing of 3. at 
room temperature for l hr. 6. Degassing of 
5. at 300°C for one hour. 
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species and that of aluminum isopropoxide suggests, in addition to the 
hydrogen-bonded species, the presence of an alkoxide type of structure 
yr (45551) Pumping at room temperature and up to 80°C removed most 

of the physically adsorbed hydrogen-bonded species, leaving the alkoxide 
type structure on the surface. The formation of such an alkoxide species 
on the electron-abstracting Lewis-acid site, mic, is in agreement 

with the work reported(45 48.51) 


The two species, V and VI, can be represented as follows, 


H 4 
| | 
HO CNH, Hc — rain 
0 0 
/ ~ i sig 
fry e 
val Al 
0 V Goal Daa 
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On heating above 100°C and up to 300°C (Curve 3), all of 
the band structure attributable to these two species disappears, 
giving propylene in the gas phase and leaving adsorbed water on the 
alumina surface. This is indicated by the strong band at 1635 cm! 
the 0-H bending vibration. For comparison, Curve 4 shows this band 
when water vapour is adsorbed on y-alumina. This adsorbed water can 


be removed by pumping, even at room temperature. The dehydration 
(18) 
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reaction can proceed as shown in reaction (5.1). 
| 
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Above 130°C, two new bands appear at 1575 and 1465 cm’! 


and 
they increase in intensity with increasing temperature. The observed 
bands could be assigned to the symmetrical and asymmetrical stretching 
vibrations of the COO groups of an acetate (5!) , The formation of a 


(51) 


carboxylate from the alkoxide species as explained by Greenler and 


Kage1 (4°) (58) 


has been subject to criticism by Fink » i.e. it seems 
unlikely that three hydrogen atoms will split from the alkoxide and 
an adjacent hydroxyl group. However, the adsorption of acetone on 
y-alumina at 220°C produced the same infrared characteristics that 
were observed for the proposed carboxylate species VIII obtained with 
2-propanol. Although acetone was not detected as an intermediate 
product from 2-propanol, the adsorbed acetone may have been present 
in the form of species VII preceding the formation of the carboxylate 
structure. Additional evidence supporting this explanation arises 
from studies with the sodium hydroxide-doped y-alumina. The doping 
eliminates or modifies active surface hydroxyl groups. Under these 
conditions, large amounts of the adsorbed acetone-like species VII 
can be found by adsorbing 2-propanol at 130°C. Upon further heating 
above 220°C, large amounts of the carboxylate species were detected. 


The following reaction scheme seems plausible for the formation of 


carboxylate: 
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The presence of methane in the gas products was confirmed 
by high resolution mass spectral analysis. In addition, the acetate 
Structure VIII could be partially hydrolyzed by water above 220°C to 
form acetic acid. These observations further support the proposed 
reaction steps (5.2). 

The formation of ether, particularly at low temperatures, 
as already reported (© | was not observed in this study. It is possible 
that the ether is present in undetectable amounts or that it rapidly 
forms the alkoxide species (79) , Moreover, the formation of ether 
from secondary alcohols is thermodynamically less favorable than the 
formation of olefins (14), 
5.1.2 2-Butanol on y-alumina 

Figure 5-2 shows the results from the adsorption of 2-butano] 
On y-alumina up to 300°C. The spectrum of physically adsorbed 2-butanol 
(Curve 2) shows changes similar to those observed with 2-propanol. 
The products from dehydration were observed in the gas phase as a 
mixture of 1l-butene, trans- and cis-2-butene. Consistent with thermo- 


dynamics, their distribution appeared in the order, 
trans-2- >) cis-2- > 1-butene 


Little of the carboxylate species was noted even at a temperature of 
300°C. Since the carboxylate species is formed from the alkoxide 
species, the reduced amount of alkoxide Present could account for this. 


Two infrared bands detected at 1561 and 1479 cm”! Suggest the presence 
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Figure 5-2. 
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Spectra on the Adsorption of 2-Butanol on 
Y-Alumina. 1. Baseline. 2. Alumina and 
2-butanol at 10 torr and room temperature. 
3. Alumina and 2-butanol at 130°C for i hr. 
4. Degassing of 3. at room temperature for 
one hour. 5. Degassing of 4. at 300°C for 
one hour. 
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of a propionate (48) . 


Analysis of the gas phase confirmed the presence 
of a propionate structure. Methane accompanied the dehydration pro- 
ducts and propionic acid was found upon hydrolysis at 220°C. These 
observations are consistent with the explanation of the acetate 
Structure formation shown in reaction (5.1). It would thus be antici- 
pated that a secondary alcohol may form some carboxylate and an alkane 
fragment upon heating in the presence of y-alumina. 
5.1.3 3-Pentanol on y-alumina 

Figure 5.3 shows results of 3-pentanol adsorbed on y-alumina 
from room temperature up to 300°C. The spectrum of 3-pentanol ad- 
sorbed on y-alumina at room temperature is similar to those for 
2-propanol and 2-butanol, showing the presence of both hydrogen-bonded 
and alkoxide species. 

On heating at 80°C, a new band is observed at 1685 cm’! 
accompanied by the formation of a trace amount of mixture of cis- 
and trans-2-pentene. This 1685 cm”! band occurs at a frequency 
which was considered to be too high to be attributable to the O-H 
bending band of water formed on dehydration. An alternative explanation 
Suggests that this 1685 cm! band is due to some partial double bond 
character, formed as an intermediate just before dehydration occurs 
at higher temperatures. This species may be represented by structure 
IX. 


A similar band at 1675 em! 


was observed on adsorbing trans- 
2-pentene on y-alumina, stable even on prolonged pumping at room 


temperature. This additional evidence supports the explanation that 
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Figure 5-3. 
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FREQUENCY (CM!) 


Spectra on the Adsorption of 3-Pentanol on 
y-Alumina. 1. Baseline. 2. Alumina and 
3-pentanol at 10 torr and room temperature. 
3. Alumina and 3-pentanol at 80°C for 1 hr. 
4., 5. AS per 3. but at 130 and 220°C, resp. 
6. Alumina and 3-pentanol at 300°C for 1 hr. 
7. Degassing of 6. at 300°C for one hour. 
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H,0 — CH = CH — CH, — CH, 


H 
H 
> IX 
0 


a partially double bond form of intermediate is formed before dehydra- 
tion occurs. Above 80°C, the dehydration proceeds at a faster rate 
making it more difficult to observe the intermediate, IX. The 1640 
cm! band due to 0-H bending of adsorbed water increases with temperature 
and above 220°C, l-pentene is formed. 

The carboxylate species is formed at 130°C along with parallel 
formation of the alkoxide species, shown by an increase in the in- 


tensity of the broad band at 1150 cm”! 


The carboxylate structure 
(in this case, propionate) remains stable even when pumping at 300°C. 
The formation of the carboxylate would, in this case, be expected to 


produce ethane analogous to the methane formed in reaction (5.2). 
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5.1.4 2-Propanol on 2% NaQH-doped y-alumina 

Figure 5.4 shows the spectra obtained from adsorption of 
2-propanol on 2% sodium hydroxide-doped y-alumina from room temperature 
to 300°C. At temperatures up to 130°C (Curve 2), the hydrogen-bonded 
species V is formed to a much lesser degree than was observed on the 
y-alumina. The elimination of many surface hydroxyl groups by re- 
action to form =Al-0-Na would account for this reduced physical ad- 
sorption. The presence of the alkoxide species VI was shown by the 
deep band at 1165 cm”! 

Above 130°C, the dehydrogenation reaction occurs forming acetone 
as the product along with minor dehydration indicated by a trace of 
propylene. The bands at 1705, 1640, 1440, 1360, 1330 and 1227 cm”! 
(Curve 3) originate with the acetone adsorbed on ats sites. The band 
at 1640 cm”! results from the adsorbed acetone as well] as from some 
adsorbed water formed during dehydration. In addition to this band 
structure, four more bands are observed at 1572, 1456, 1330 and 1170 


] 


cm”! (Curve 3). The pair of bands at 1572 and 1456 cm are attributed 


to the carboxylate species, similar to those found on y-alumina. The 


1 which also always appear 


remaining two bands at 1330 and 1170 cm 
as a pair and occur in the carboxylate frequency range, were tentatively 
assigned to a second type of carboxylate structure. On pumping from 
room temperature (Curve 5) to 300°C (Curve 6), the band structure due 

to acetone completely disappears, leaving only the four carboxylate 
bands. Adsorption of acetone on the doped alumina at 150°C showed 


the presence of both the adsorbed acetone and carboxylate species. 
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Figure 5-4. 
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3000 2800 1800 1600 1400 1200 1000 
FREQUENCY (CM-:') 


Spectra on the Adsorption of 2-Propanol on 
2% NaOH-doped Y-Alumina. 1. Baseline. 

2. Doped alumina and propanol at 10 torr and 
room temperature. 3. Doped alumina and 
propanol at 300°C for one hour. 4. Doped 
alumina and acetone at 10 torr and room temp. 
5. Degassing of 3. at room temperature for 
one hour. 6. Degassing of 5. at 300°C for 
one hour. 7. Degassing of 4. at 300°C for 
one hour. 
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Upon heating the ketone bands diminished in intensity along 
with an increase in intensity of both pairs of carboxylate bonds. 
Further, mass spectral analysis of the degassed phase showed only 
methane to be formed, in substantial agreement with the proposed 
mechanism shown in reaction (5.2). On this basis, the assignment 
of the additional band pair at 1330 and 1170 cm | to a second carboxy- 
late structure seemed logical. 

The presence of these two pairs of bands may be the result 
of differences in charge distributions and in Sy inetry considerations 
for oxygen atoms and their energy of bonding to the surface (104) | 
By doping y-alumina with sodium hydroxide, it was observed that the 
high frequency hydroxy] groups disappeared. This may change the 
charge distribution of oxygen atoms which are responsible for the 
carboxylate formation. 

Acetone was adsorbed on fresh sodium hydroxide-doped y- 
alumina catalyst and the system was then heated to 130°C (Curve 4) 
to compare adsorbed acetone and its behavior, with the intermediate 
which led to the formation of acetone. This spectrum is very similar 
to that observed on heating the catalyst in the presence of 2-propano] 
(Curve 3). Upon degassing the catalyst exposed to pure acetone at 
300°C, only the carboxylate structure remained (Curve 7). This con- 
firms the formation of an adsorbed acetone-like species, VII, like 


that obtained from 2-propanol which on desorption gives acetone. 
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5.1.5 2-Butanol on 2% NaQH-doped Y-alumina 

Figure 5.5 shows the spectra of 2-butanol adsorbed on 2% 
sodium hydroxide-doped y-alumina from room temperature to 300°C. At 
temperatures up to 80°C, the alkoxide species is formed with very little 
hydrogen-bonded species appearing (Curve 2). Near 130°C, the dehydro- 
genation reaction starts to form species VII as shown by the increased 
bond intensity at 1730 cm! attributed to the C=0 stretching fre- 
quency (Curve 3). After heating at 300°C (Curve 4), the carboxylate 
Structure spectrum became dominant. Heating of adsorbed 2-butanone 
on a newly prepared catalyst wafer gave spectra similar to those from 
2-butanol (Curve 5). Both 2-butanol and 2-butanone produced carboxy- 
lates with bands identical to those observed with 1-propano) (48) on 
the doped y-alumina. 
5.1.6 3-Pentanol on 2% NaQH-doped y-alumina 

Figure 5-6 shows the various spectra of 3-pentanol adsorbed 
on 2% sodium hydroxide-doped y-alumina from room temperature to 300°C. 
At room temperature (Curve 2), the main Species is again the alkoxide 
type. The dehydrogenation reaction occurs at a still lower temperature, 
80°C, and continues up to 300°C with the formation of the stable car- 
boxylate species (Curves 3,4,5,7). Adsorption of 3-pentanone (Curve 
6) formed by heating shows spectral characteristics comparable to 


those expected from reaction (5.2). 
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Figure 5-5. 











3000 2800 1800 1600 1400 1200 1000 
FREQUENCY (CM"!) 


Spectra on the Adsorption of 2-Butanol on 
2% NaOH-doped y-Alumina. 1. Baseline. 

2. Doped alumina and butanol at 10 torr and 
room temperature. 3. Heating 2. at 130°C 
for one hour. 4. Heating 3. at 300°C for 
one hour. 5. Doped alumina and 2-butanone 
at 10 torr and room temperature. 6. Degas- 
sing of 4. at room temp. for one hour. 

7. Degassing of 6. at 300°C for one hour. 
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Figure 5-6. 


FREQUENCY (CM°') 


Spectra on the Adsorption of 3-Pentanol on 
2% NaOH-doped y-Alumina. 1. Baseline. 

2. Doped alumina and pentanol at 10 torr and 
room temperature. 3. Heating 2. at 80°C for 
one hour. 4.,5. As per 3. but at 130 and 
300°C,respectively. 6. Doped alumina and 3- 
pentanone at 10 torr and room temperature. 
7. Degassing of 5. at room temp. for 1 hour. 
8. Degassing of 7. at 300°C for one hour. 
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5.2 Catalytic Properties of y-alumina Using IR Adsorption Technique 
5.2.1 Effect of NaOH Concentration on Bronsted-type Acidity on the 


Alumina Surface 

Figure 5-7 shows the spectra of sodium hydroxide-doped y- 
aluminas in the OH Stretching region (3500-3850 cm |), Figure 5-7A 
shows the spectrum of the pure y-alumina. The related spectra for 
varying levels of NaOH doping are shown in Figures 5-7B to 5-7G. Three 
characteristic bands are observed for pure y-alumina at 3785, 3720 
and 3680 cm! Increased doping results in the disappearance first 
of the 3785 and then the 3680 cm”! bands, the 3720 cm”! band remaining 
largely unaffected. Doping at the 1% NaOH level, eliminated the 
3785 cm”! band, but the intensities of the 3732 and 3680 om™! corres. 
ponded almost identically to those for their counterpart frequencies, 
3720 and 3680 cm | On pure y-alumina. Elimination of the 3785 cm”! 
band also resulted in a slight shift from 3720 to 3732 cm! for the 
middle hydroxy] band, implying increased freedom for its stretching 
frequency. The presence of excess NaOH could be observed in the spectrum 
of the 15% NaOH sample (Figure 5-7G) in the form of bands characteristic 
of NaHCO... The high frequency 3785 cm”! band decreased in intensity 
very rapidly with doping. Since the doping was carried out in the 
"liquid-phase", the acidity of the hydroxy] groups, commonly associated 
with these three bands can be inferred to be Stable at normal conditions, 
and may be considered in terms of the groups' relative reactivity with 
the NaOH. The evidence Suggests that the 3785 cm! hydroxyl group is 


the most acidic, followed by the 3680 cm”! band. The continued presence 
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of the 3720 cm! band irrespective of the level of doping implies that 
it exhibits little, if any, acidity or else, that its acidic behavior, 
thus defined, is suppressed by steric factors or low reactivity. 

Additional insight into the acidic behavior may be inferred 

by studying the catalytic selectivity of y-alumina towards dehydration 
or isomerization reactions. Pines and Haag26) suggested that alumina 
exhibited either strong or weak types of acidic character. Their 
work suggested that the total number of hydroxyl-type acid sites in 
y-alumina can be measured by dehydrating 1-butanol, and that the 
fraction of these which are strong acid sites can be determined by 
cyclohexene isomerization. Heating in 5 torr of cyclohexene at 200°C 
for one hour, sample A (pure A1,03) was found to catalyze the cyclohexene 
isomerization but sample B (1% NaOH) was found to be catalytically 
inactive for this reaction; thus the strong acid sites at 3785 cm”! 
were eliminated by doping. Proceeding further, both samples A and B 
catalyzed the dehydration of 1-butano] under identical reaction conditions, 
B to a lesser extent, hence the weak acid sites alone would be expected 
to remain on the series of doped catalyst, and probably decrease 
with increased doping. This point-of-view is consistent with the 


observation that the 3785 and 3680 cm! 


hydroxyl bands correspond, 
respectively, to strong and weak Bronsted-acid behavior of the surface 
hydroxy] groups. 

Additional infrared spectra were obtained for the hydrated 
surface of y-alumina by slowly drying a wafer of pure y-alumina at room 


temperature and under vacuum over a 48-hour period as shown in Figure 


5-8. Starting with a hydrated surface whose IR spectrum exhibited no 
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Figure 5-8. Spectra on the Degassing of Y-Alumina at 
Room Temperature over Varying Time Periods: 
A, 0.5 hr.; B, 4 hr.; pee Te. 
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bands characteristic of free hydroxy] groups, pumping led first to the 
appearance of the 3720 cm”! band and then to the 3680 cm! band; 
additional pumping having no further visible effect. Since the degree 
of hydrogen-bonding between surface hydroxyl groups and water mole- 
cules on the surface is somewhat proportional to the "acidity" of the 
hydroxyl groups, this also Supports the view that the 3720 cm”! band 
exhibited weaker acidic character than the 3680 cm”! band. The high 
frequency 3785 cm”! band was not observed at room temperature pumping 
indicating that its strongly acidic character resulted in stable 
hydrogen-bonding of surface water molecules. 

As a result of the NaOQH-doping, the hydroxy] groups likely 
formed structures of the type, = Al-0-Na, since the various IR spectra 
did not reveal bands attributable to other types of bonding. 

5.2.2 Effect of NaOH Concentration on Lewis-type Acidity on the 
Alumina Surface . 
Lewis-acid sites are believed to be abundant on the dehydrated 


37) | probably as incompletely coordinated aluminum 


alumina surface | 
ions. On this basis, Pines and Haag'26) also suggested the possibility 
that doping with NaOH could give rise to an acid/base surface comp lex 
of the type shown in reaction (2.9). 
I I 
Al + NaQH > [-A1-OH]" Na* (2.9) 


To determine whether reaction (2.9) had occurred upon doping, it is 


necessary to measure whether the number of aluminum ion sites decreases 
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11] 
with increased NaOH doping. The aluminum ion sites are generally re- 
garded as being Lewis-acid sites. In the presence of pyridine, 


(37) reported that the infrared spectra of the coordinately 


Parry 
bonded species (on Lewis-acid Sites) should differ from those of the 
pyridinium (on Bronsted-acid sites). In this work, neither pure 
alumina nor the doped aluminas exhibited IR spectra of pyridinium 
ions on their surfaces hence their Bronsted-acid sites are too weakly 
acidic to react with pyridine. On the other hand, the spectral in- 
tensities suggested that the amount of Pyridine bonded on the Lewis- 
acid sites remained constant with increased NaOH doping. This evidence 
verifies that only the Bronsted-acid sites (two types of hydroxyl group 
sites) reacted with the NaOH to form = Al-0-Na and thus reaction (2.9) 
does not occur. 
9.2.3 Effect of NaOH Doping upon Dehydration Catalytic Activity 

Since the mechanisms for dehydration or dehydrogenating of 
different lower secondary alcohols were found to be virtually 
identical for each alcohol (see Section 5.1), 3-pentanol was chosen 
as the model reactant to study the catalytic activity and selectivity 
of NaOH-doped y-aluminas. Figure 5.9 shows the results obtained from 
adsorbing 3-pentanol on the various NaOH-doped y-aluminas after heating 
at 200°C for one hour. The catalytic reactions which occur are either 


the dehydration, 


(C,H.) CHOH + Pentene + H,0 (5.3) 
a5 9 2 
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Figure 5-9. Spectra on the adsorption of 3-pentanol on NaOH doped 
y-aluminas at 200°C for one hour. 
Sample numbers correspond to those given in Figure 5-7. 
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or the dehydrogenation, 


(Cols), CHOH + (CAH) CO + Ho (5.4) 


2 


and the relative amounts of reactions (5.3) and (5.4) which have pro- 


| band, from the 


] 


ceeded may be estimated by comparing the 1715 cm 
carbonyl group of the adsorbed 3-pentanone, or by the 1635 cm ~~ band 
for adsorbed water. The spectra for pure y-alumina (Figure 5.9A) 
Shows that only dehydration is catalyzed. The spectra for 8% and higher 
doping levels (Figures 5.9F and 5.9G) show that only dehydrogenation 
has occurred. Increasing doping levels of NaOH content, to between 
5 and 8%,catalyzed varying amounts of both reactions (Figures 5.9B 
to 5.9E). The selectivity of the catalyst towards dehydrogenation 
increased with increasing NaOH content. The influence of excessive 
NaOH doping was also checked by contacting 3-pentanol vapor with a 
wafer of pure NaOH at 200°C. The results showed that neither NaOH 
nor NaHCO. exhibited catalytic activity and hence the role of NaOH in 
reactions (5.3) and (5.4) is essentially that of a "doping" agent. 
5.2.4 Effect of NaOH Doping on Catalytic Selectivity 

The distribution of chemical products 3-pentanone and water 
for series of tests to evaluate the influence of NaOH doping on catalytic 
selectivity in reactions (5.3) or (5.4) was determined by gas chromato- 
graphic analysis, using an eight-foot "Porapak T" (80-100 mesh) column 
at 160°C and a thermal conductivity detector. The relative activity 
of each catalyst sample (100 mg A1,0,) was studied by introducing 5 mm 


Hg of 3-pentanol vapor into the infrared cell at 200°C with a contact 
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time of one hour. Consistent and reproducible spectra and GC analyses 
were obtained in this manner. The experimental results are shown in 
Table 5-1, expressing activities in terms of mole % of 3-pentanol re- 


? of the catalyst surface (see Table 3-2 for surface 


acted per 10 m 
area of each catalyst sample). This activity criterion provided a 
means of correcting for loss in catalyst specific surface area attri- 


butable to increases in NaOH content. 


Table 5-] 
Activity and Selectivity of Y-Al 403 or NaOH-doped Al,033 
Dehydration or Dehydrogenation of 3-pentanol at 200°C 


Catalyst NaOH % Of 3-pentanol *% Of 3-pentanone 
Sample Content converted per 10 m formed per unit 
% catalyst surface total conversion 

A 0 87 0 

B 1 78 0 

C 2 6] 38 

D 3 47 5] 

E 5 16 96 

8 13 100 

G TS 1] 100 


The product distributions show that the % of 3-pentano] which 


has chemically reacted decreases with increasing NaOH doping. In addi- 
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tion, the % of 3-pentanone formed per unit total conversion (via com- 
bined dehydration and dehydrogenation) increased from zero between 
one and two percent NaOH content to 100% between 5 and 8% NaOH content. 
These results quantitatively demonstrate the shift in catalytic selec- 
tivity from dehydration to dehydrogenation, and the decline in catalytic 
activity per unit surface area with level of NaOH doping. 

Assuming that the surface hydroxy] groups act as the cata- 
lytic sites for dehydration (see Section 5.1), comparing the spectra 
Of Figure 5-9 with the product distributions in Table 5-1 suggests 
that the hydroxyl groups with the 3680 cm”! band (samples B to Bard 


possibly, those with the 3785 cm™! band (sample A) catalyze dehydra- 
tion. When these two bands vanished upon doping (Samples F and G), only 
dehydrogenation catalytic activity appeared. Although the 3785 cm”! 
hydroxyl groups could not be isolated from the 3680 cm! groups, their 
greater acidity suggests that they should be more active towards de- 
hydration than the 3680 cm”! groups. 

Figures 5-9C to 5-9F show that the selectivity for dehydration 
decreases with the decreasing intensities of the 3680 cm”! hydroxy] 
band. The 3720 cm! hydroxyl sites are either not accessible to NaOH 
doping or are catalytically inactive towards dehydration since their 
band intensities did not change with increased doping. 

5.2.5 Surface Structure of y-alumina 

Several interpretations of the nature of active sites on the 
alumina surface are available. One view proposed (45) is that the 
aluminum or oxide ions exposed at the surface act as active sites and 


that the surface was completely free of hydroxy] groups. This view 
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is contradictory to a conclusion of this study based upon the infrared 
spectra of alumina being characteristic of surface hydroxy] groups in 


the region, 3600 - 3800 cm |, 


The current widely accepted model, pro- 
posed by Peri (42) | assumes a random configuration of hydroxyl groups 
after dehydration, leaving adjoining residual oxide ions and oxide 
vacancies, and exposed aluminum ions. The variations in the number 
of oxide ions surrounding the hydroxyl groups differentiate them into 
several types and imparts varying Bronsted-type acidity to each of the 
three groups. The three hydroxyl bands at 3680, 3720, and 3785 cm! 
are thus interpreted. 

The general view of the structure of y-alumina surface, 
based upon crystal chemistry proposed by Verwey (22) , suggests that the 
unit cell of y-alumina, a tetragonally formed spinel, consists of 32 
oxide ions with 21-1/3 aluminum ions arranged at random in the 16 octa- 
hedral and 8 tetrahedral positions of the spine] structure. It is _ 


(30) that the surface is terminated in anions and the outer- 


believed 
most aluminum ions lie below the surface plane. These aluminum ions 
provide the electrophilic centers by which water is dissociatively 
held on the surface; the electron pair of the OH being positioned 
over the aluminum ion while the dissociating proton forms a hydroxyl 


31) 


via binding to lattice oxygens | The structure of the hydrated 


aluminum surface can be represented by the following, 
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The lattice OH is positioned on the surface below the plane of the OH 
adsorbed over the aluminum ion. It seems reasonable to assume that 
there are differences in the energies of interaction of surface 

hydroxyl groups with octahedrally and tetrahedrally coordinated 
aluminum ions, and that these differences result in the appearance 

of different isolated infrared bands, characteristic of hydroxyl groups, 
between 3680 and 3785 cm”). 


Peri and Hannan (4! > 105) 


used the exchange reaction between 
deuterium gas and surface aluminol groups to determine the number of 
such groups on their sample. They determined from the IR spectra ob- 
tained during exchange with deuterium that the average infrared ab- 
Sorptivity of the hydroxy] groups was of the order of 8 (10°) cm°/mole.. 
They also suggested that Beer's law could Probably be used to compare 
relative amounts of hydroxy] groups at coverages below 15%, where no 
hydrogen bonding appears to exist. Applying this suggestion to the 
catalysts used in this investigation, the highest frequency OH groups 
were calculated to number about 15% of the lowest frequency groups 
(Figure 5-7A); in reasonable agreement with Pines and Haag( 26) who re- 
ported the fraction of strong acid sites numbered roughly 10% of the 
total acidic sites. Although the absorptivity calculations agree with 
Peri and Hannan's view, their interpretation that the higher acidity 
is associated with the lower frequency hydroxyl groups is not consistent 
with our findings based upon the chemical reaction of Bronsted-acid 
sites with NaOH solutions. 

Lewis-acid sites, such as the incompletely coordinated aluminum 


ion, are believed to be abundant on the dehydrated alumina surface (3/) , 
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Assuming their presence as well as that of the Bronsted-acid sites, 

their relative numbers and corresponding relative catalytic activities 

can be estimated. Using 2 (10/3) site/cm* (106) and 10 (10!3) site/cmé (26) 
for Lewis-acid and Bronsted-acid sites, respectively, a dehydrogenation 

to dehydration selectivity ratio of 1 to 5 would be anticipated for 

both types of site which are equally active. Comparing sample A (total 
Bronsted-acid) and sample B (weak Bronsted-acid) to sample E (total 
Lewis-acid) in Table 5-1,it can be estimated that a value of 15 for the 


% Of 3-pentanol converted per 10 me 


Surface would be approximately 
equal to the value for the complete dehydrogenation of 3-pentano] on 


the fully doped alumina. Thus, 


Fractional conversion of 
Number of total] 
Bronsted-acid sites 3-pentanol on pure Al,04 ou eS a 


= dae 
Number of total “ Fractional conversion of 15 ~ 7 


Lewis-acid sites 3-pentanol on fully doped 


and 


Fractional conversion of 
Number of weak 
Bronsted-acid sites 3-pentanol on 1% NaQH-doped A150. (ome 2 


Number of total ° Fractional conversion of ies ie al 


Lewis-acid sites 3-pentanol on fully doped Al,03 


Since these ratios are of the order of magnitude to be expected, the 
activities per site do indeed appear to be roughly equivalent, and 
differences in catalyst selectivity appear to be the result of relative 
Site accessabilities. Furthermore, the roughly 1% NaOH content of 


sample B eliminated the strong Bronsted-acid type of sites (3785 cm’! 
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hydroxyl band). The relative selectivities can thus be used to esti- 
mate the strong Bronsted-acid sites; they number roughly 10% of the 
total Bronsted-acid sites, (2:83.2 8-9. 3-2) 10.3, the figure reported by 
Pines and Haag. The bees, level for doping of Bronsted-acid 

Sites corresponded to a NaOH content Slightly in excess of 5%. The 
amount of NaOH used in samples B to E did not correlate with the esti- 
mate of numbers of Bronsted-acid sites eliminated by doping and cal- 
culated using Beer's Law. This may be due to the non-uniform dispersion 
Of NaOH on alumina surface during the drying process. 


5.2.6 Mechanisms of Dehydration and Dehydrogenation of 3-pentanol 


on y-Al,0, and NaOH-doped y-A1,03 


Adsorption of 3-pentanol on pure y-Al,03 at temperatures up 
to 150°C resulted in the formation of two types of surface species: 
(X) 3-pentanol Strongly hydrogen-bonded to the surface hydroxy] 
groups and, (XI), an aluminum pentyl oxide type of structure chemi- 
sorbed on Lewis-acid sites. These Structures can be represented by 


the following, 
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The structure (X) was proposed in Section 5.1 based upon its infrared 
spectrum. When the NaOH doping level was increased with the resulting 
elimination of hydroxy] groups, the peak intensity of the 3500 cm | 
band decreased, suggesting a reduction in hydrogen-bonding of type (X). 
These results suggest that the H atom of the hydroxyl group in the 
alcohol is free and that the alcohol molecules are hydrogen-bonded to 
the surface hydroxyl group through the oxygen in the alcohol. 

The interpretation of structure XI has also been discussed 
(Section 5.1) in terms of the alkoxide species, characterized by the 


deep band at 1165 cm), 


Since the intensity of this band and thus 
the amounts of the surface alkoxide species were found to be identical 
for all of the catalyst samples investigated, the view is re-enforced 
that the Lewis-acid sites, which are responsible for the formation of 
species XI, are not affected by NaOH doping of the alumina. Thus, it 
is unlikely that the alkoxide species, XI, appears as an intermediate 
in the dehydration step. The participation of the surface alkoxide 
species in the dehydration of alcohols was also rejected Pocantty ee) 
on independent grounds. Thus, a hydrogen-bonded alcohol molecule seems 
to be the most probable adsorption state from which the elimination of 
water is initiated. 

Several different views of the mechanism of the dehydration 
of alcohols on alumina surfaces exist 89) but this work supports the 


previous reaction schemes (Section 5.1) for dehydration and dehydro- 


genation reaction mechanisms: 
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oe (525) 


XII is a mixture of cis- and trans-2-pentenes 


2""5 
|x = — | 
0 H 0 
| | : 6.6 


Reaction (5.6), involving dehydrogenation requires that an alkoxide 
reaction intermediate be formed. The observed hydrogen-bonded species 
(X) and the alkoxide species, (XI), thus are involved in the dehydration 
and dehydrogenation mechanisms, respectively. They, in turn, involve 
Bronsted-acid or Lewis-acid surface sites, respectively, and their 
activity and selectivity on various NaOH-doped aluminas conform to 


the preceding discussions. 
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5.3 Simultaneous Kinetics and Mechanisms 
In this study 2-propanol was chosen to be the reactant be- 

cause it gives only two reaction products, propylene and water, in 
the dehydration step, and acetone and hydrogen in the dehydrogenation 
step. According to the results in Section 5.2, pure y-alumina, or 8% 
NaQH-doped y-alumina could be used to catalyze either dehydration 
alone or dehydrogenation alone, respectively. 
5.3.1 Blank Run 

To determine the threshold temperature at which the catalytic 
effects of the reactor surfaces became Significant, the temperature 
of the reactor (without a catalyst wafer) was raised in a stepwise 
manner over a period of 10 hours; during this time 2-propanol was 
continuously passed through the reaction system. At a reactor temper- 
ature of 400°C and a reactant pressure of 0.9 cm Hg, conditions at 
which appreciable reaction occurred in the presence of a compressed 
catalyst wafer, no reaction was detected in the empty reactor by IR 
and GC analyses. Since the maximum temperature employed in this work 
was always less than 400°C, the catalytic effects, if any, of the 
reactor surfaces may be neglected. It is concluded that in the kinetic 
studies all of the observed reaction occurred on the catalyst surface. 
9.3.2 Mass Transfer Considerations 

Interpretation of experimental results become difficult when 
mass transfer resistances are significant in comparison with chemical 


(107) 


resistances. Recently, Carberry reported an extensive survey of 
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the effect of mass transfer upon observed reaction kinetics. Since 
catalytic activity is usually exhibited as a function of some key 
experimental variables, the effect of these key variables in each rate- 


controlling domain is summarized in Table 5.2. 


Table 5-2 
Behavior of Various Kinetic Parameters 


in Each Domain of Gontroy 6107) 


ee Se ee St eet ree tircul se 


. Activation sate Fluid 
Domain Energy Order Pellet Size Velocity 

Chemical reaction E n independent independent 
Pore diffusion - E (nt+1) a : 

reaction a 5 (-r))=F(7) independent 

3/2 

Bulk mass ¥ . Peer ya = -4,,\ 1/2 
transfer ~0 first = (-r,)=F(7) x(velocity) 


~<a aes eee ee 


*]/L = external surface area-to-volume ratio for the wafer. 


Prior to considering the reaction kinetics for the dehydration 
of 2-propanol, the influence of mass transfer on the rates of reaction 
must be evaluated. In this regard, two possible contributions to 
diffusional retardation of the reaction rate exist, namely bulk gas- 
film diffusion and pore diffusion: 

(1) Gas-film Diffusion 


To establish the absence of concentration gradients between 
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the bulk gas phase and the catalyst external surface, it was necessary 
to determine the relationship between the bulk flow rate of the re- 
actant in the reactor and the conversion to products. In this study, 
since the rates of reaction of 2-propanol were obtained directly, 
the experimental determination of the mass transfer effect was greatly 
simplified. If the rates of reaction were unaffected by film diffu- 
sion, an increase in bulk flow rate in the recirculation loop (keeping 
the reactant feed rate constant) would have no effect on the conversion 
of the reactant (Table 5-1). Accordingly, the speed of the recircula- 
tion pump was varied to change the mecycle i lowarrom 13.2 0 21.7 
liters per minute and finally to 32.8 liters per minute. The con- 
version was found to unaffected by the volumetric flow rate in the 
recycle loop at the reactant partial pressure of 0.479 cm Hg and a 
reaction temperature of 261.7°C. Since these values represented the 
highest reaction temperature and lowest reactant partial pressure em- 
ployed in this study, the bulk diffusion limitation was assumed to be 
negligible in all of the experimental runs. 

(2) Pore Diffusion 
The effects of pore diffusion can be determined by measuring 
reaction rates for different sizes of catalyst particles under other- 
wise identical conditions ‘112) | The catalyst effectiveness factor 
approaches unity (absence of pore diffusion effect), when no increase 
in reaction rate per unit quantity of catalyst occurs on further sub- 


division of the catalyst pellet. Accordingly, three different weights 
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of wafers, and correspondingly, three different catalyst wafer thick- 
nesses were investigated. The results, plotted in Figure 5-10, show 
that the reaction rates (see equation 4.5) were identical for both 
catalysts of 0.1325 and 0.3024 gram weight. It was concluded that no 
mass transfer limitations were present for a catalyst wafer of less 
than 300 mg weight. 

For the series of runs at 252.7°C, if any of the dehydration 
rates were limited by diffusion of the reactant into the wafer, the 


diffusional limitation would be expected to be most severe for large 





RT (-r,) 
values of the Thiele diffusion modulus o = Cor ae gt and 
D(1/L) PA 


hence for experimental runs at the largest ratio of reaction rate 


(-ry) to partial pressure fae 


Calculations, based on equations 
given by Satterfield !08) | showed that the effectiveness factor for 
pore diffusion was close to unity, and thus the intra-particle con- 
centration gradient was negligible, for all the rate data obtained 
in this study. The detailed calculations are given in Appendix D. 

For dehydrogenation reaction, similar results were obtained. 
It is thus concluded that pore diffusion limitations were negligible 


for both catalysts (pure y-alumina and 8% NaOH-doped y-alumina) when 


weights were less than 300 mg. 


33 Stability of the Catalyst Activity During the Decomposition 
of 2-Propanol 


Successful kinetic measurements require reproducible catalytic 
activity during the various experiments. The activity of the 2N NaOH- 


treated alumina catalyst on both dehydration and dehydrogenation at 
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Figure 5-10. 


Dehydration of 2-Propanol 


At 252.7°C over Y-A1,03; 


Catalyst wafer weights: 
4 @ 0.1325 gm. 
b 0 0.3024 gm. 
C @ 0.5257 gm. 


N. flow rate: 9,378 g-mole/hr . 
Recycle rate: 32.8 1/ min. 





20 30 


w/F, g-hr./g-mole 


Evaluation of Pore Diffusion. 
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465°C was previously reported‘ !9) to decrease with process time. This 
was explained by the deposition of elemental carbon on the catalyst 
surface and deactivation by water absorption in the catalyst. How- 
ever, in this study the catalyst activity was found to remain constant 
during 0.5 to 60 hours of use. Since the time of use for any charge 
of catalyst in the various runs never exceeded 48 hours, the deactiva- 
tion of the catalyst during the experimental runs was believed to be 
negligible. For illustration, the dependence of activity, in terms 
of fractional conversion of 2-propanol, on time for both dehydration 
and dehydrogenation is shown in Figure 5-11. Surface spectra were 
also recorded at both one and sixty hours of the reaction time. The 
spectra of 2-propanol adsorbed on pure y-alumina show the presence 
of hydrogen-bonded 2-propanol] as well as the carboxylate surface 
structure. The intensity of the peaks associated with these two 
species remained identical after 60 hours of reaction time. Hence 
at reaction temperatures below 400°C it is safe to assume that the 
catalyst activity will remain unchanged for at least 60 hours. For 
the case of dehydrogenation on NaQH-doped y-alumina, the surface spectra 
show the presence of carboxylate species which again gave identical 
spectra at one and 60 hours of reaction time. 
9.3.4 Consideration of Capillary Condensation 

In essence, capillary condensation occurs when the vapor 
pressure of the reactant and/or products over a liquid-phase concave 
meniscus in a catalyst pore is lower than that of the bulk liquid-phase. 


If this difference is sufficient to lower the vapor pressure to the 
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total pressure of the system, the small pores will be filled at total 
pressures lower than the saturation pressure of the bulk liquid. As 
the pressure increases, larger pores will also be filled. At the 
vapor pressure of the bulk liquid, of course, the pore volume will be 
completely filled. 


Bonneomte) 


» while studying the dehydration of alcohols at 
lower temperatures, listed several points of difficulty which should 
be considered in obtaining kinetic data. He has shown that some of 
his data were complicated by the influences of both pore diffusion 
and capillary condensation effects. However, in this study, the 
lowest temperatures employed for the dehydration and dehydrogenation 
were 220°C and 283°C, respectively. For dehydrogenation, the reaction 
temperatures were above the critical temperatures of 2-propanol (235°C) 
and acetone (235.5°C). For dehydration, the reaction temperatures were 
near the critical temperature of water (374°C) and above that of 
2-propanol. Moreover, infrared spectra of the catalyst surface under 
reaction conditions show no liquified propanol, acetone and water. 
The capillary condensation effects on the reaction rates are thus be- 
lieved to be negligible. 
5.3.5 Side Reactions 

(1) Dehydration 

During the study of catalyst stability (Section bas Steno 

Side reactions were evident from the infrared and gas chromatographic 
analyses. One side reaction which is thermodynamically possible is 


the dehydration step to form ether. The absence of this side reaction 
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is in agreement with the observations of deMourgues et ay (16) for the 
same reaction system. 

(2) Dehydrogenation 

In the case of dehydrogenation, the dehydration reaction 
would be the potential side reaction. No evidence was observed for 
the presence of water and propylene (dehydration products) by both IR 
and GC analyses. It is thus assumed that only dehydrogenation reaction 
occurred on the 8% NaOH-doped y-alumina. 
9.3.6 Kinetics and Mechanism of Dehydration of 2-Propanol 

Infrared spectra of the adsorbed species were obtained under 
reaction conditions. In light of the infrared observations, various 
chemically and mechanistically plausible kinetic models were proposed 
and the kinetic data were used to establish by statistical discrimina- 
tion methods which kinetic model was the most satisfactory one. In principle, 
this would appear to be the usual approach of modern kineticists, however, 
the mechanistic data largely eliminated the majority of the kinetic 
models usually considered on a priori basis. The following sections 
enumerate the details of this approach. 

(1) Spectra of Adsorbed Species 

Infrared spectra of the species adsorbed on y-alumina surface 
under the reaction conditions were recorded for each kinetic run. 
Figure 5-12 shows the typical spectra of the adsorbed species after 
eliminating the base-line spectrum due to Al,03. The close similarity 
of each set of spectra, as shown in Figure 5-12, suggest that the re- 


action mechanisms are identical for all runs. 
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Figure 5-12. 





IR Spectra of Adsorbed Species on 
Y-Alumina at Steady-state Reaction 
Conditions. 
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In the study of dehydrogenation of ethanol on copper catalyst, 


Bradshaw and Davidson reported !!0) 


that the rate-controlling step of 
the reaction changes from surface reaction at lower temperatures to 
adsorption of ethanol at higher temperatures. Thaller and Thodos (111) 
obtained the rate data for the catalytic dehydrogenation of 2-butanol 
to 2-butanone and concluded that the rate-controlling step of the re- 
action at 600°F shifts between surface reaction-controlling and desorp- 
tion of hydrogen-controlling steps, depending upon the experimental 
pressures. In this study, if the reaction mechanism changes as a 
function of either reactant pressure or the reaction temperature the 
infrared spectra of the adsorbed species should also change accordingly. 
However, no evidence for such shift in mechanism was observed for the 
dehydration of 2-propanol on y-alumina. 

Spectra obtained at the reaction temperature of 252.7°C show 
the presence of hydrogen-bonded 2-propanol at 1465, 1373 and 1170 cm” | 
and carboxylate structure at 1575 and 1465 cm |, At this temperature, 
the infrared absorption bands due to the hydrogen-bonded 2-propano] 
increase in intensity with increasing partial pressure of 2-propanol] 
(Figure 5-12). However, the intensity of the bands due to carboxylate 
species remains constant throughout the period of investigation, j.e. 
between one and sixty hours of reaction time. This is also an indica- 
tion that the carboxylate species, although very stable on the alumina 
surface even at high temperatures, do not poison the catalyst surface. 


No bands at 1635 cm™! (due to water) and between 3000 and 3200 cm! 
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(due to olefins) are observed in the surface spectra. This Suggests 
that under the reaction conditions, water and Propylene were not ad- 
sorbed on the catalyst surface to any significant degree. 
(2) Kinetic Models 

For any given reaction system, it is generally assumed that 
there exists a rate model by which one can describe adequately al] 
of the phenomena actually occurring in the system. In the case of 
heterogeneous catalysis, it is a common practice to postulate not 
one, but a number of plausible models which are in agreement with 
the overall] chemistry and, then, to discriminate between these rival 
models with the aid of experimental data. Thus, the success in rate 
modeling a given system is vitally affected by the degree of postulation 
involved in the models to be screened and, of course, is directly af- 
fected by the effectiveness of the model discrimination technique. If 
an adequate model is not included among those models proposed for the 
System under study, any sophisticated model discrimination method will 
generate a meaningless correlation. In fact, all statistical discrimina- 
tion techniques are developed on the basic assumption that one adequate 
model must be contained among the rival models to be discriminated. 

Basically, two model-postulation methods are most commonly 
employed in the correlation of kinetic rate data for solid-catalyzed 
gaseous reactions: 

(a) The over-all reaction rate is correlated empirically as a 


function of the operating conditions (e.g. power-function models). 


i abies Seon s16 2borltom nofdstut2oq=febom ows wvtleotase tes 
deena ne ate ee v0 fi 
iT Mic 





.ersyeq2 ststive off nf bevisede: ee 6) : 


ab oft yd betostts ylletiv. ef mateye nevtp 6 ‘en rab | 


- af j ~ / . 4 
tenimivaetb [sbom oA to 22enevitostts ont " ot3 


aT 
: 


hatsnimiroetb sd ot efsbom (avit odd proms’ benfatnos od “eum Tab 


dinoesb bo “Sne 5 Wot yd fabom: sien avetns ove 


is : i: 
: : me eral t: 
1 



















7) 






i 


fyqor bis istsw , onots totes noftoasey ent sebm” ‘8 ie 


pent? 


tnsottinpte yne oF 906% we Jeufetss ond no’ be dyoe 
2febot attonty™ (s) at 


tt .ms eve notsose7' Nn nevtp yrs 107 


maseya oft af enty 119990 vl fsusos snemonsitg of 
»fiDevq NoMMOD 5 et +t . 2teylstso evosnaporst: of 

. {w Bfsbom efdteustq Yo +edmwn 6 tud 1 Fam 
orscimris2tb ot nend bis vrtetmeno Itsyevo"4 9 


Be 


AT =. 6d8b TsFrvamt ex to bts arid rit Pw at ab on 


s 


7 


/ 


S200 To , bMS Ldililanove ad oF efsebom sit at bev ove 


aS 


32073 proms babufont gon ef Fabon’ os sup nate 

= * ng 7 
s2tb fabom betsotiafriqo2 ae «but2 abn ITeYe 
‘le ,dost nl .nobtefsvi03 aeatgninean 6 ott ne 


tomuezs ofesd sid no bisdefavel 31h zsuptiitaat” nord 







c 


wcmapieiiry sit # beyol en 
is EftO yube) ree 


a 


as 


134 
(b) Rate equations are developed which describe a rational 
mechanism by which the reactions May occur, and the best-fitting re- 
lationship is assumed to describe the reaction mechanism. 
The first method has a limited scope of application since 
it is not safe to extrapolate beyond the range of the experimental 


data( 112). 


The usefulness of this type of analysis is also reduced 
by the fact that little insight into the steps occurring during the 
reaction process is obtained. As a result, this approach will not be 
considered in the study of dehydration of 2-propanol. 

In the second method, the rate data are fitted by a model 
describing a plausible reaction mechanism. The most widely employed 
models to describe the mechanism of gaseous reactions on solid surfaces 
are Hougen-Watson type (113) models. This class of models assumes 
that the reactants are adsorbed on one type of active site on the 
catalyst surface where they react, either with themselves, with other 
adsorbed species on neighboring sites, or with components in the gas 
phase adjacent to the catalyst surface. Further assumptions in formu- 
lating Hougen-Watson rate models include the existence of a single 
rate-controlling step, such as adsorption, surface reaction, or de- 
sorption, within a mechanism involving a series of reaction steps, 
and the absence of interactions between adsorbed molecules. 

Although the assumption of a single rate-controlling step 
greatly simplifies the derivation of the rate expression, rate equations 
that involve more than one rate-controlling step should be considered 


for some reaction systems to obtain a more adequate representation of 
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(110,114) | For example, Bischoff and Froment ‘!!5) 


the experimental data 
reanalyzed the rate data obtained by Thaller and Thodos (111) and pro- 
posed a rate model of Langmuir-Hinshelwood type in which two consecu- 
tive steps, surface reaction and desorption of ketone, are both rate- 
controlling. Despite the lack of experimental justification to 
indicate that these two steps are comparatively slow, the proposed 
model appeared to be adequate. 

The assumptions applied in the formulation of Hougen-Watson 


(113) can be modified in various ways. An obvious ex- 


rate equations 
tension of the assumptions is a heterogeneous surface with different 
types of adsorption sites which produce non-interacting chemisorbed 
gases. For the reaction system studied by Thaller and Thodos ‘111) , 
Kolboe (116.117) proposed a rate model] that involves two chemically 
different chemisorption sites and showed that the goodness-of-fit of 
rate data is improved noticeably compared to models based upon a non- 
interacting, uniform surface. 

Although the approach followed in developing rate equations 
of the Hougen-Watson type is straight forward, the task involved in 
determining the best rate expression for a complex reaction system 
is immense. If all possible reaction models are listed within the 
assumptions of the Hougen-Watson scheme, a relatively large number 
of models is obtained. For example, 84 models have been postulated! 18) 
for methane oxidation and 47 models for ethanol dehydration! 19), 

The presently available model discrimination techniques can not be 


easily applied to a group of rival models, as large as 47 or more, 
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due to the prohibitive amount of computer time that would be required. 
Consequently, it is desirable if not necessary, to eliminate obviously 
inadequate models before utilizing the better model discrimination 
procedures. Not infrequently, the results of the discrimination pro- 
cedure turn out to be ambiguous in that no one model is shown to be 
clearly superior. However, a generally effective screening technique 
is not yet available in the literature. 

The infrared spectroscopic technique has provided consider- 
able insight into the identity of chemisorbed species which may act 
as intermediates in catalytic reactions, and also the nature of the 
interaction between an active site and an adsorbed molecule. This 
information, obtained at the molecular rather than the macroscopic 
level , should serve as a good Screening method because the mechanisms 
proposed according to infrared data would incorporate more detailed 
knowledge of the mechanism taking place on the catalyst surface. 

The infrared study of adsorption of 2-propanol on y-alumina 
Surface at temperatures varying from room temperature to 300°C, has 
been discussed in Section 5.1 and a reaction mechanism has been pro- 
posed based on the observed adsorbed species at the different tempera- 
tures. According to the proposed mechanism, two types of active sites 
are involved in the reaction sequence, surface hydroxy] groups and 
Surface oxide ions. This is also in good agreement with other evidence 
from the literature (see Section 2.6). Both single-site and dual-site 


mechanisms of the Hougen-Watson type rate equation have been proposed 
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in different studies (see Section 2.7) for the dehydration of alcohols, 
however, in formulating the dual-site Hougen-Watson model there seems 
to be no justification for treating the two types of active sites as 
being identical (reaction 5.1). A dual dissimilar-site mechanism was 
therefore used in extending the Langmuir-Hinshelwood assumptions. 

The assumptions are generally valid if the reactant is adsorbed on 
highly selective sites so that the surface coverage is rather low. 
The model thus derived may provide an appropriate approximation of the 
true mechanism. The low surface coverage of 2-propanol under reaction 
conditions is supported by the infrared spectra of the adsorbed species. 
The absorbance of the IR band at 1373 cm”! due to adsorbed 2-propano] 
increased from 0.057 to 0.65 when the temperature of the cell was 
lowered from 252° to 150°C for a 2-propanol partial pressure of 1.9 
cm Hg. The monolayer coverage would be expected to give an absor- 
bance for the band at 1373 cm much higher than 0.65. Since the maxi- 
mum absorbance observed in this Study was 0.33 at 1373 cm | the 
Surface coverage for all the kinetic runs was believed to be quite 
low. 

The plausibility of involving only surface hydroxy] groups 
as active sites is also proposed according to the following representa- 


tion of chemical steps, 




















vel | | | | 
z2ionoots to morsetbyieb ott vo? (1,8 noftoge 992) estbude Ingrstttb, at 7 
amase ayant Tebom nozteW-noguol Stizq|eub ond paitefumot ah, .vevewod i 


2s aoti2 svitos to eaqyt owt odd potdeowh 10d noltepttiseut, onned oF) 
25w petnsrionm st fe~vel fnteetb Ceub A ‘Cae aptapess) feottnebi pated 
.enotsqmees boow/l eran tH-niaipaal os. entbnadxa nt bozu svotevedd : 
no badvoebs: 2t tnbtasey odie yl fsv9nep. 916 2notiqmuers, ; 4 
wor. vsdter 21 sp619VO9 Soest we ghd tend oz 2edte evttoatae vided 
oft to nofdsmixorggs stsiyqorages ne obhverq yam beviveb. 2udd t=bom ¢ 
notsoce7 yabou fonsqoyq-S to spsigyeo sastw2ziwof eT ie nish em 
sjvag2 badio2bs sd3 to sqtosge DeveNtat geld ud: badvoqque et cosile 
fansgova-S badve2bs of ‘ub e mo €{éf ts bred #1 ond to sansdro2ds,.9A » 
agw I'fao ant to ovutiracniad ett noriw 2.0 oF $20.0 mort b 
e.1 to sruezexg Tetdwq TonegongeS s tot °08E of, °S8S. moxk, boxe 
(o2zds M6 eis ot betasqxe ‘sd bi vow lige 19.5 [ onom sdT « + Qh 7 _ 
-ixem sit sonte .23.0 mene yerokd tou! “a AXEL $s. bred alt yor 
ait ,~ mo EXEL, $6 S60 250 ybui 2 tds of bovesedo sonscrozds & 
atiupiad 6 bavatted 26W env sitontd odd ffs NOt SP6BT9VOD 9: 
: | hat Gnaeplanberl 

 aqvorp Tyxovby" sostawe yi no Fanon to sane ek 
saatoaptan eniwor Tot gilt ot entbroass, ReEaReRs: eels. at.gotle avid i . 


tt ADs naa moors 
; ihe 


138 


H H 
| | 
H ei aii — CH, HAC = CH - CH, 
H----- 0 > H (5 7) 
H 
ais SY pee 
H 
NN ' 
H 
0 
| a; 
aN | 
Al 
es 


The generalized mechanisms for single-site and dual dissimilar-sites 
mechanisms are given below: 


(a) Single-site mechanism 


At ) AL, (5.8) 
At; =s P+ We, (5.9) 
WE = Wt hs (5.10) 


where A, P, W and dy represent 2-propanol, propylene, water and 
surface hydroxyl groups, respectively, and propylene is assumed not 
to be adsorbed. The latter assumption is based upon the failure to 
detect adsorbed propylene even at room temperature by the infrared 
technique. Equations (5.8) to (5.10) imply a reversible reaction 


mechanism but a feed consisting of a mixture of propylene and water 
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does not yield any detectable alcohol. Consequently, one of the above 
three steps of the single-site mechanism must be irreversible. Since 
adsorption and desorption of gaseous reactants or products are in- 
herently reversible steps, the surface reaction (step (5.9)) is there- 
fore the most likely step to be irreversible. The irreversible form 


of step (5.9) of the generalized single-site mechanism becomes , 
AL, > P + We, (5.11) 


The rate expressions derived for the single-site mechanism are listed 
in Table 5-3 for all of the possible combinations of the rate-con- 
trolling steps. 

It can be seen from Table 5-3 that the assumption of a 
single rate-controlling step, as proposed by Hougen and Watson, is 
not necessary for this dehydration reaction if one of the reaction 
step is irreversible. This may be the reason why Hougen-Watson 
models generally express the experimental data well, even if the 
reaction may involve more than one step rate-controlling. 

(b) Dual dissimilar-sites mechanism 
From the representation of the chemical reaction steps, 


reaction (5.1), the generalized mechanism can be written in the form, 
Ath, = Ae, ee 


AL, + Ry > PL + We, (5.13) 
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Table 5-3 


Possible Kinetic Models for Dehydration of 2-Propanol, A +P + W 


A = 2-Propanol; P = Propylene; W = Water 


(a) Single-site mechanism 


] Ly 


(2) AR, ae Tew. 


] 


] 
(3) Ww = Wt h, 


ee 


ae Rate Expression Assigned 
Controlling Step(s) (-ry) NodeTeNoe 
ee ee ee 
biP, 
(1) (THBP) yF 
b.P 
LA 3 
(2) cae m2 
I+b5P,+b5P\, 
(3) Same form* as M2 
(1) & (2) Same form as M2 
(2) & (3) Same form as M2 
(1) & (2) & (3) Same form as M2 


a 


*"Same form" implies different constant terms but an identical 
functional form 
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Table 5-3 (continued) 


Dual dissimilar-sites mechanism 
(4) A+ 4 = Ag, 
(5) AL, + 25 > PL + We, 


(6) PL — P+ Ly 


= 
(7) we, S Wee 


2 2 


eee 


; Rate Expression Assigned 
Controlling Step(s) Ty MadeliNa. 
ae Sues ee ei no feney lene into 

| (4) Same form as MI] 
Die 
LA 
(5) (1#55P,) (145?) HS 
b,P 
LA 
(6) (67, )(Hb SP )eb,P, MM 
(7) Same form as M4 
(4) & (5) ene eAeora etn ot 
4) & (5 M5 
I+b5Py I+b,P,, +b yP, 
(5) & (6) Same form as M4 
(5) & (7) Same form as M4 
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(5.14) 


mo] 
20 
— 
=) 
ss 
20 


We, => WHE 


where &.'s represent surface-active sites with Ly = surface hydroxy] 
group, Lo = surface oxide. The step (5.13) is again assumed to be 
irreversible as per the arguments employed for the single-site mechanism. 
The desorption constant of Propylene for step (5.14) is assumed to be 
very large since no detectable amount of Propylene was adsorbed on 
the catalyst surface and, also, introducing propylene into the feed 
had no effect on the observed reaction rate. The possible rate ex- 
pressions for this mechanism are also given in Table 5-3. The com- 
bination of three or more rate-controlling steps is rejected for the 
following reasons: 

(a) The rate expressions become exceptionally complex and the 
available model discrimination methods cannot be employed effectively; 

(b) Water and Propylene were not observed in the surface spectra 
under reaction conditions and thus desorption of water or propylene 
is unlikely to be the rate-controlling step. 

The derivation of rate equations for the single-site mechanism 

is similar to that given by Hougen and Watson $113) | In the case of 
the dual dissimilar-sites mechanism, a sample derivation of the rate 
equation, assuming step (5.13) to be rate-controlling, is given in 


Appendix E. 
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(3) Model Discrimination 

The literature available on model discrimination procedures 
is extensive. Bard and Lapidus ‘9°) give an excellent review of the 
methods available at the present time. In this study, a Bayesian 
approach suggested by Quon and developed by Singh (28) at the University 
of Alberta has been adopted to discriminate among the proposed rival 
kinetic models. 

A set of ten dehydration experiments was conducted in order 
to obtain the best model from the proposed five rate equations, MI] 
to M5. The kinetic data are presented in Table 9-4. From the IR 
calibration curves (Tables 4-2 and 4-4), it was found that the measure- 
ments of the propylene partial pressure were less accurate than those 
of 2-propanol. Since no side reactions were detected, the partial 
pressures of water and propylene in the reactor were computed from 
the difference of the partial Pressures of 2-propanol in the feed and 
in the exit stream. This resulted in a substantially better fit 
than if the partial pressures of water were obtained from the propylene 
calibration. 

It was assumed, initially, that each of the five models was 
equally likely. A 20% prior probability was thus assigned to each 
model. The following criteria were used for model rejection: 

(a) lack of fit of the model to the data, as indicated, by an 
extremely small posterior probability of the examined model, and 


(b) unacceptable characteristics of the estimated parameters. 
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Table 5-4 


Experimental Data for Dehydration Reaction 


Catalyst weight = 0.1383 gm 
Feed composition = 2 propanol with No at .378 g-mole/hr 


Reactor temperature = 252.7°C 


Run No. peaeon fate PAO PA Py P Material* 
hr-g-catalyst cm Hg cm Hg cm Hg balance, % 
eg eee Ee Pe Phe QO 
] 0.03947 1.044 0.479 0.570 100.5 
2 0.06355 1677 0.812 0.859 99.6 
3 0.08569 2.31] 1.162 1.158 100.4 
4 0.12663 3.578 1.907 1.759 102.5 
5 0.16365 4.845 2.700 2.193 101.0 
6 0.19739 6.113 34533 2.49] 98.5 
7 0.27938 9.597 5.958 3. SiO 98.7 
8 0.40000 15.934 TOs 721 D3! 99.8 
9 0.8851 60.110 52.016 11.434 95.0 
10 1.0105 51.090 40.367 10.416 100.5 
11** 0.2901 47.951 40.124 


*Material balance a\P + Pi )/Pao x 100 
*kwith 45% of H,0 in the feed 
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In case (b), for instance, estimates of adsorption/desorption con- 
Stants are generally required to be positive and to possess negative 
temperature coefficients. If several models have acceptable para- 
meter estimates, the model with the highest posterior probability, 
Say larger than 0.95, was selected. 

Incorporating all of the ten data points in Table 5-4 by 
the Bayesian expected likelihood method, the following posterior 


Probabilities were obtained for the five different models: 


Po,10 = 9-476 5 Pz 49 = 0.524 ; P 


1510" 481098 


where the first Subscript refers to the model] number asSigned in 

Table 5-3 and the second Subscript to the number of data points in- 
corporated. Thus, model M3 is regarded to be the most likely candi- 
date but by no means conclusively; actually, both M2 and M3 can be con- 
Sidered equally favored by the ten data points employed. The need 

for further experiments in order to obtain a better discrimination 
between M2 and M3 is apparent. 

It was decided to design the future experiments on a rational 
basis. The criterion of Singh?2) as given in Section 2.10.2 was emp loyed 
and the region yielding the maximum value for the expected change in 
entropy was chosen in obtaining additional data points. Table 5-5 
gives some of the results of this exploration by assuming the variance 
of experimental measurements to be 0.001. The maximum entropy change 


occurs in the region of Py = 50 cm Hg and Py = 30 to 50 cm Hg. Since 
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Table 5-5 
Expected Entropy Change at Various Partial Pressures 


of 2-Propanol (A) and Water (W) for Dehydration Reaction 


ne nr see ee 


Expected Entropy 


Pa » cm Hg Pw » cm Hg Change - R 
er 
10 5 5.8 x 107° 
20 5 0.264 
30 5 0.43 
40 5 0.40 
50 5 Vireo Vf 
10 20 0.094 
20 20 0.154 
30 20 0.209 
40 20 0.285 
50 20 0.396 
10 30 0.09 
20 30 0.18 
30 30 0.26 
40 30 0.35 
50 30 0.46 
10 40 0.08 
20 40 0.18 
30 40 O27 
40 40 0237 
50 40 0.476 
50 50 0.476 
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the partial pressure of water with a feed of Pa = 50 units would be 
expected to be somewhere between 10 and 15 cm Hg, it was decided to 
inject additional water into the feed. Data point 11 in Table 5.4 was 
thus obtained at Pa = 47.951, Pi = 40.124 and (-ry) = 0.2901, ata 
temperature of 252.7°C, with 45% water in the feed. The results ob- 
tained after including this additional data point strongly favored the 
dual dissimilar-site mechanism, M3. The results are given below: 

(a) Single-site mechanism, M2 
Parameters by =$0%124 bo = 0.0234 by = 0.4292 
95% limits + 0.025 + 0.0130 + 0.1360 
Variance (0°) = 2.59 x 1077 
Model likelihood a 1.42 x 10/3 

(b) Dual dissimilar-site mechanism, M3 
Parameters b, = 0.09544 b, = 0.007856 by =e 0825747 
95% limits + 0.00800 + 0.0020 + 0.04830 
Variance (0°) = 4.126 x 107° 
Model likelihood a 3.436 x 10/7 


Variance - covariance matrix: 


1.595 x 10° © -2.858x 10° 9.13 x 197 
-2.86 x 1076 9.81 x 1077 te Ae 


9.13:x-1057%- 7 22.15 x 1072 5.83-x 107° 


The posterior probabilities for M2 and M3 are respectively equa] to 0.00004 
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and 0.99996. Since the posterior probability of M2 was extremely low, 
the single-site model was rejected. The rest of the models were also 
examined with this additional data point. The results showed that 

the posterior probability of each model is much less than 0.001% and 
therefore, they were also rejected. 

From Table 5-3, it can be seen that M2 also represents the 
combined steps, surface reaction and desorption of water, to be rate- 
controlling. However, it was found that water was not present on 
the catalyst surface at any detectable level under the reaction con- 
ditions. Thus, it is clearly shown that the dehydration of 2-propanol 
proceeds via a dual dissimilar-site mechanism in which the surface 
reaction is the rate-controlling step. 

It is interesting to make a comparison between the conven- 
tional dual similar-site Hougen-Watson model and the dual dissimilar- 
site rate equation developed in this thesis. The results are given 
as follows: 


(a) Dual-site Hougen-Watson model, surface reaction-controlling 


b,P 


; PA o* = 1.387 x 1074 
(1+b,P,+bP,,)“ Posterior probability < 0.0001% 





(b) Dual dissimilar-site mechanism, surface reaction-controlling 


byPa of = 4.126 x 10°" 


= ee 


(1+b,Pq) (1+b5P,)) Posterior probability > 99% 


148 


a eid ih : 
- r : ' 
Wy ey i s* 

: Jt } 


wot yfamsrdxo esw SM to yttl idsdog -rofit20q ent -eonte 2200.0 bs 3 
o2is svew 2fobom oft to 204 set .bagosper 26w febam atte-stente: odd 
tert bawode 2tluzey ait? .sefoq edb fenott tbbs ats dd iw bontmsxe c 
bis 8100.0 nent ezaf doum ef Tebow i358 ro yi tdsdorw, notiaszoq-oltt ; 
-betsatey o2zfs oven yond «store ’, 


BPI 

























os sdnsesngen oats SM tent nese od nap st eal: stdsT mort. Sole 

-9t51 od ot ,vatew to noftqveesh bis nohioset sostwe 2qate beri dr . 
fo ¢nseevq ton eew vetew Jadt bnuo? zew tt .18v9woH pnb Forse 

| -f09 noftossy sft tebay Fave ol detostab yas Js sostwe 4eyiatso & 
lonsgovd-S to nottavbydeb sav Jedd wore ulseolo af tr , 2udT a 
aoetrue oft dotdw nt me toedosm ‘otte~veftmtzetb Isub a stv ab 


agate enti fovsnoa-ss4 ‘odd et nota 


mac 


~ - 


at 
. 


-nsvnoo ent nsouted noe tysqueD 6 stem of eniteswint aeal Mi 
-“elimbe2ib Psu oft bre fabom no23 Bh ~ ~neguol at hawus! tke Taub f ) 
nevip eis etfueey AT 2hesitd otay nt beqof evab nots supe otst she 
es sewoltot @ 

ontfforvinoo~notsoss sostwwe ,fabom innit nsquol stte~faud- (sp 


Lit. 
{ Ai) 


Of x TWE.fow Ap ad he 
BH000.0. > yttlidédowg yoinereot “een \ = ¥ 


fh. 


q 7» ro 7A uv? 


7 ej vais ee ‘i 7 
Ln i Tap vx a 


149 

Thus, it may be concluded that assuming a dual-site mechanism without 
a priori knowledge about the nature of active sites could lead to 
poor formulation of kinetic models. 
9.3.7 Activation Energy for the Dehydration of 2-Propano] 

The kinetic data obtained at different reaction temperatures, 
Runs 12 to 18, are presented in Table 5-6. After having finalized 
the reaction mechanism for the dehydration of 2-propanol from isothermal 
rate data at 252.7°C, an estimation of the activation energy was at- 
temped by applying the same reaction mechanism to other temperatures. 
This assumption was based upon the observation that infrared spectra 
of the adsorbed species exhibited the same spectral band characteristics 
as those obtained for the 252.7°C runs. Therefore, the mechanism 
was believed to be applicable in the range of temperatures, 120°C to 
171°C, studied. By the usual Arrhenius rate expression, the following 
kinetic equation can then be used( 112) for the temperature range studied: 


b, 0 exp (- E,/RT) Pa 


(-ry) = (5.16) 


H 
[I+b, 4 exp(- Rr PqlLitb, 9 exp(- Py] 


The statistical analysis can not be expected to give a good estimation 
for the six parameters in equation (5.16) without a very large number 
of data points over a wide range of reaction temperatures. Since it 
was found that the dehydration reaction was extremely slow at temper- 
atures below 110°C and that the side-reactions were Significant at 


temperatures above 350°C (94) | kinetic runs could be performed 
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Table 5-6 
Kinetic Data over Range of Temperatures , 


for Dehydration of 2-Propanol 


Catalyst weight = 0.1489 g 


Feed rate = 0.02386 g-mole of 2-propanol per hour* 


ae Se ape a 


Reactor Reaction Rate (x 10°) 


temp. °K -mole ) 
hr-g-catalyst 
ee ee eae tet ee 


Run No. Pas cm Hg Pie cm Hg 


12 393 1.904 2.286 0.025 
13 409 3.917 2.256 0.055 
14 413 5.091 2.240 0.071 
15 418 5.830 2.230 0.081 
16 427 8.27] 2.196 0.115 
17 439 13.091 2.129 0.182 
18 444 15.204 2.100 0.211 


Ln i. 


*With No flow rate = 0.378 g-mole/hr 
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15] 
only in the range of 110°C - 300°C. This resulted in a restriction 
for obtaining a large number of data points for parameter estimation. 
It was decided to keep Pa and Py deliberately low to ensure the de- 
nominator terms, boP a and b3P\s in model M3 were small compared to 
unity. Therefore, (1+b5Pq) and (1+b,P\) in model M3 would be expected 
to be much less sensitive to temperature than the specific rate con- 
stant, b,. The low values of Pa and Py will also give a uniform 
temperature over the catalyst surface. As an approximation, bo and 
b3 were treated as constants, bo = 0.007856 and b3 = 0.25747. Equation 
(5.16) is now rewritten as equation (5.17) by assuming that the 
Arrhenius form of temperature-dependency for the rate constant was 
applicable to this reaction system, 


6) 9 exp (- iP, 


- Yr (5317) 
At” HBP I+b, W 


The parameters , by 0° the so-called frequency factor and E.. the 


yt 
apparent activation energy, were estimated by the same methods used 


in the previous section. The results with 95% confidence limits are 


= (1.276 + 0.001) x 10° ; E, = (1.46 + 0.05) x 104 —cal 


: g-mole ° 





0 


For illustration, the log of the specific reaction rate constant 
(logy, b, b,) is plotted against the reciprocal of the absolute temperature 
(1/T) as shown in Figure 5- 13(112) and a fairly linear plot is obtained. 


From this plot, the apparent activation energy is calculated to be 


. 
>) Pa 
Aowi 


far yey 

















a 


notiotitaer 5 nf betfuesy atiT .9°00E ~ J°Ofl Yo sonst ott ait Mi 

Nottsmisz2s votemersq 10? etntog sis to vednyn’ ‘spiel 6 onintstdo ‘ ‘ 

-9b ofJ siwens ot wot vistetedhisb’ yi bes. ‘) qaoa ot bsbiosb 26w 1 

ot beisqmoa [Tem2 svsw SM febou di ied bas. ais! « zartsd “10: 
betosaxe sd bluow & Isbom at (yiedr) bas (,%ga+l} ,s10t STENT. 

-09 93617 If thosqe shit nerd s"1vss19qmst od ourereaba 2esf Adum so 

tro tau: 5 vig ocis [i tw y9 bas 49° %0 eaufsy wot oft «yd etnst: 

bas .d ,coftemixorggs nis 2A eoatiue teyfstso et ysv0 onud549 0 

ottsupl .{MT2S.0 = -¢ bas B8BT00.0 = 4 {2#niatenog! és badaerd. sie : 

od tend pnimveeg yd (Xt.2) moftsupe 28 node hye won at (af.a) 

26W tnssenon S357 oft- — yonebrsqeb-swisisqnes to aot 2ut 

smasey2 nottossy etait of orate 


~ 


3 Mf | 
alte | 
at he. 7. . 
inao TEE Narada 4 a 
wee’ SANS <5 
+s iS OC). 3 


ia . os 
21F «3 bas tot2st Yoneupeit belfso-o2 ont ‘9 po «ers3 
. 


bs2u zbodtem smsz ont yd batenttes etsw Npvons. nottevitos. $1 
o16 etimt> sonsbitnes xée dgtw ettveer ait .nottse2 euotverq 


_ 






Ve * Oe Te gf asi 
2 tame OF x (20.0 + apr) =, om x (100.0 + ats. tts “yd 
: : A bs y | er 


ae a 


¥ - ~ yale : ® in a . = 
7 : = ? : 4 ¥ q . 
> : * 


: { } i a ¢ ms ai 
: iP» ae oan ~ > Sa = 7 ae 8) 
t es io PY F We - - 
is . fl rt, - - 


152 





*jouedoirg-z jo uotzeazpAyeqd 
Jo o3ey oyQ Jo Aduepusedseq seanjerzedusay 


t-7e ‘OT X L/T 


"€1-G eanbtd 


qsAjTe jes 6-wo-1zy/aetou-5 "(e°ezx¥a)/ (“4-) 








a 
24 





Ly 


Sa - 
} 


sae 


— 


‘ 


7 = * 
ae ee s 8) 
F 5 
a ‘ 
. 
7 
6 jt . 
cm cole ig - 
; 
: 
7 





to sts eft to yonsbneqsed emssteqmeT .£1-2 exupit 
-ionsgoxrd-S to noitsibyded = 


Sei 


| 

15 K cal/g-mole which is very close to the value of 14.6 obtained by 
Singh's parameter estimation method '98) . 

It can be seen from the derivation Of model M3 in Appendix E 
that by = LyLoKak. or, in terms of Arrhenius relationship over a 
range of temperature, b, 5 LL OKy ok3.0 exp[- (Hy + E_)/RT]. There- 
fore the true activation energy, Er for the dehydration should be 
the sum of apparent activation energy, Ey» and the enthalpy change 
for adsorption of A, - Hy. Since - Hy is approximately equal to 
10 K cal/g-mole( 16) | ES is thus estimated to be 14.6 + 10 = 24.6 which 
is in good agreement with the value of 25 reported in the literature (16). 
5.3.8 Kinetics and Mechanism of the Dehydrogenation of 2-Propanol 

Suitable kinetic and mechanistic studies with which to com- 
pare quantitatively the rate expression for the dehydrogenation of 2- 
propanol on 8% NaOH-doped y-alumina do not exist. Nevertheless, the 
methods for treating the kinetic data are similar to those discussed 
in the case of dehydration. 

(1) Spectra of Adsorbed Species 

Infrared spectra of the adsorbed Species on 8% sodium hydroxide- 
treated y-alumina under Steady-state reaction conditions are Shown in 
Figure 5-14. The alkoxide structure and carboxylate species were 
again observed on the catalyst surface. The infrared bands at 1560, 
1445, 1360 and 1155 cm”! were attributed to the carboxylate species 
and the bands at 1545, 1360, 1155 and 1125 cm”! were attributed to the 
alkoxide species, as was found in the room temperature spectra obtained 


in the adsorption Study (Section 5.1). Since the IR band at 1705 cm! 
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CATALYST WT.=0.2747 gm. 
TEMP ERATURE=343 .2°C, 


PARTIAL PRESSURE OF 


2-PROPANOL, cm Hg. 


1.04 





1600 1400 1200 1000 
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Figure 5-14. IR Spectra of Adsorbed Species on 8% 
NaOH-doped +Alumina at Steady-state 
Conditions. 
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which corresponds to the stretching vibration of C = 0, was not ob- 
served in the surface spectra, it was concluded that acetone, one of 
the dehydrogenation products, was not adsorbed on the catalyst under 
the reaction conditions. Hydrogen, in the form of H which splits off 
from 2-propanol was proposed(45) to be adsorbed on the surface oxide 
ion to form a surface OH group. The parent alcohol thus formed an 


adsorbed alkoxide. The reaction can be represented as follows: 


| 
BeS | | 
/ H y. Al 
# NS 
Al 0 
eS. 


The presence of the newly formed OH groups due to the reaction (5.17) 
may be evident from the increase in peak intensity of the surface 
spectra in the region of 3600 to 3800 cm! which are attributed to 


OH stretching vibrations (34) 


However, a strong, broad hydrogen- 
bonding peak centered at 3550 cm”! developed when alkoxide and ee 
late species were present on the catalyst surface; thus, the change 
in peak intensity in the OH stretching region cannot be effectively 
measured. 

Again, as in the case of dehydration, the spectra of the 


alkoxide species showed an increase in intensity, with increasing 


partial pressure of 2-propanol although to a lesser extent than the 
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intensity increase of the hydrogen-bonded alcohol in the case of de- 
hydration. The bands at 1560, 1445, 1360 and 1155 cm”! are attributed 
to the surface carboxylate species, and are observed to remain un- 
changed for reaction times of one to sixty hours. These results 
Suggest that the carboxylate species had very little effect, if any, 
on the catalyst activity. Since the alkoxide structure, and the 
carboxylate species, were the only components observed on the catalyst 
Surface for all of the kinetic runs, it was believed that the reaction 
mechanism remained unchanged for those runs. The mechanism proposed 
for the dehydrogenation from the adsorption study (Section 5.1) seems 
applicable to the high temperature kinetic runs. 

(2) Kinetic Models 
From the study of infrared spectra of the adsorbed species, 
the possible reaction mechanisms for the dehydrogenation reaction can 


be represented as follows: 


Scheme A, Single-site mechanism 


CH CH 


3 | 3 
eee) es He —» H.C —C—CH 
2 | ; 3 | 3 I: 3 + Ho (5.18) 
0 0. 0 
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Scheme B, Dual dissimilar-sites mechanism 
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Scheme C, Alternate dual dissimilar-sites mechanism 
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Scheme D, Dual dissimilar-sites mechanism, reduced from Scheme B, 
by assuming 2-propanol to be dissociatively adsorbed simultaneously on 


the oxide and the aluminum ion. 
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, H 
| 
HC —C— CH. H3C —C— CH. H3C — C — CH, 
| <= | Aa 1 | 
0 0 0 treed 
\ ale | 
H 0 Al 0 Al 
Dar Ale aD J ve ee 
ee oe 


— H3C — C —CH + H (S219 


grene'0 Tene 0 
LE Ta 
The single-site mechanism (Scheme A) is obtained by assuming 
that only the surface aluminum ion participates in the reaction. Al- 
though the infrared spectrum of the surface alkoxide is very similar 
to that of aluminum isopropyl oxide, the presence of the structure [XIII] 
(a loose hydrogen atom bonded to the alcoholic oxygen) in the reaction 
step (5.18) is still possible. A feed of acetone vapor mixed with 
hydrogen does not produce any detectable amounts of alcoho]. It is 
thus assumed that the dehydrogenation reaction is also essentially an 
irreversible reaction. The adsorption of alcohol and the formation 
of alkoxide are assumed to be reversible steps because the alkoxide 
species can be desorbed to give the corresponding alcohol on evacuation! 4) , 
Also, acetone is found to adsorb at lower temperatures. The irreveri- 
bility of the dehydrogenation reaction thus should result from the 
surface reaction step (decomposition of the alkoxide species). 
The generalized mechanism of Scheme A consists of the follow- 


ing steps: 
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A+ = AL, (5.22) 
AL, ~ Ka, - H, (5.23) 
Ke, Sa K+e, (5.24) 


where A, K and a represent alcohol, ketone and surface aluminum 
ion, respectively. 

The participation of surface aluminum and oxide ions in the 
dual-site mechanism resulted in several possible reaction mechanisms 
and this has been already discussed partly in Section 5.1. Schemes B 
to D show these possibilities. In order to develop rate equations 
for all possible combinations of rate-controlling step(s) for each 
mechanism, it is more convenient to rewrite these mechanisms in the 
generalized form with the additional notations, Xo = Surface oxide 
ions, and XL, = surface alkoxide species. 


Scheme B 


A + Ly = AL, (5.22) 
AL, + ky = Xu, + Heo (5.25) 
Xe, + dy * Kg. + He, (5.26) 


2H2, = Hp + 2h (5927) 
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Scheme C 
A+ 2%, =e Ag, (5 
AL, +h) = X21 + H, (5. 
XQ, + ky * Kk. + H2, (5 
2Hk. = Hy + 22, (5. 
Ki, —— Ly (5 
Scheme D 
A+ 2, +h = X21 + H2, (6. 
Xk, + 25 * Key + He, (5. 
2Ht, = Hy + 22, (5 
Ke, = K+ dy (5. 


Since acetone was found not to adsorb on the surface under 


reaction conditions and since the addition of acetone to the liquid 
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feed stream does not inhibit the reaction rate, it is assumed that the 
desorption constant for step (5.24) is very large. Using this assump- 
tion, the possible rate expressions derived for the dehydrogenation 

in a similar way to those of dehydration are summarized in Table 5-7. 
A sample derivation is given in Appendix E. It should be noted that 
in the case of Scheme C, the derived rate expressions are highly non- 
linear and complex, except when the surface reaction (step 5.26) alone 
is assumed to be the rate-controlling step. A sample derivation, 
showing the complexity of this kind of rate equation, is again given 
in Appendix E. It is thus expected to be very difficult to perform 

a reasonable parameter estimation and mode] discrimination with exist- 
ing statistical methods. However, the surface spectrum of 2-propanol 
on 8% NaOQH-doped y-alumina shows the presence of alkoxide species but 
not of hydrogen-bonded 2-propanol. This is supported by the spectra 
of adsorption of 2-propanol at room temperature and then evacuation 

at the same temperature. In both cases the IR spectra showed the 
same characteristic bands. Since hydrogen-bonded alcohol, but not 

the alkoxide-species, was found to be removed upon evacuation at room 
temperature, these results suggest that only the alkoxide species, 

not the hydrogen-bonded 2-propanol, was formed as a reaction inter- 
mediate on this catalyst. This is also in agreement with the results 


(85) who showed that decomposition of hydrogen-bonded 


of Knozinger 
alcohol leads only to the formation of olefin. Therefore no attempt 
was made to analyze the rate equations for the Scheme C mechanism 

which involved hydrogen-bonded alcohol as one of the reaction inter- 


mediates. 
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In the case of Schemes B and D, the rate expressions obtained 
are also very complex whenever the combination of two hydrogen atoms 
(step (5.27)) is assumed to be involved in the rate-controlling step. 
(See Appendix E for a sample derivation). Since some earlier studies 74 >120) 
have indicated that hydrogen is dissociatively adsorbed on alumina 
surface only to an extremely small extent, it is presumed that the 
amount of adsorbed hydrogen in atomic form is so small that it is not a 
rate-controlling consideration. 

(3) Model Discrimination 

Model discrimination for the dehydrogenation reaction was 
also achieved by the Bayesian approach as was used in the case of de- 
hydration. The prior probabilities of the nine proposed kinetic models, 
N] to N9, were taken to be 1/9 for each, showing no initial preference 
for any specific model. Twenty-four kinetic data points were obtained 
at the conditions as shown in Table 5-8. These data were examined to 
determine which of the reaction mechanisms was best fitting. The re- 
sults are listed in Table 5-9. Models N1 and N2 were rejected because 
these two equations indicate that the reaction rate is independent of 
hydrogen pressure and therefore give a very poor fit compared with the 
other models. Among all of the other models examined, N3, N5, and N8& 
were found to represent the experimental data quite well (see Table 5-9 
for model probabilities). 

The parameters estimated for N3 gave a value for the coeffi- 
cient of Pre bos equal to zero. In terms of adsorption-desorption 


coefficients for the mechanism of Scheme B, derived with the adsorption 
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Table 5-8 


Experimental Data for Dehydrogenation Reaction 


Catalyst weight = 0.2747 gm 
Feed composition - 2-propanol except when specified* 


Reactor temperature - 343.2°C 


er 
Run R 


eaction rate (x 10°) p 
A? 


No. -mole 
hr-g-catalyst 





19 3.49 0.362 0.111 
20 4.52 0.529 0.134 
2] 5.42 0.698 0.156 
22 6.22 0.870 0.174 
23 6.98 1.040 0.194 
24 i .67 i215 0.209 
25 8.34 1.386 0.228 
26 8.94 1.564 0.244 
roll 9.58 Poo) O=207 
28 10.65 2.091 0.283 
29 Te 7t 2.440 Oslo 
30 I2e7) Gel og 05336 
3] 13.66 3.156 Oma09 
32 14.56 3.514 0.38] 
33 15.40 0/3 0.402 
34 16.25 4.231 0.424 
2 20.26 6.222 0.524 
36 26.6] 10.052 0.686 
37 BOcoy 20.001 1.000 
38 47.89 29.280 12227 
39 58.28 42 1.492 
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Table 5-8 (continued) 





Run Reaction Rate P p 
P,, cm Hg Ho K 
No. (x 1000) A emHg cm Hg 
40* 34.85 16.306 0.895 13615 
4|* 64.72 51.662 1.656 6.998 
42* 30.22 iZn62e 0-777 Zable 
43** 72.59 19.972 0.397 0.397 


ne ——————— 


*(1) feed composition: 10% acetone in 2-propanol 
(2) Py is computed from Pao - Pa 
2 
** Catalyst weight = 0.0586 gm 


+ No flow rate = 0.378 g-mole/hr 
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Table 5-9 


Model Performance for Dehydrogenation of 2-Propanol 


ose Wo. Yartnee (A) SRST 
N3 1.59 x 1079 81% 
N4 2 iTuxslOae 1% 
NS Reducible to N8 
N6 3.45 x 107? ~ 0 
N7 2,48 x 10°° = 0 
N8 1.75 x 107° 18% 


No 3.36 x 10°” at 
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of alcohol and surface decomposition of alkoxide assumed to be rate- 


controlling, the constant by is equivalent to, 


_ Kogll + Kos) 
1 + aKa. 7Ko7 


BY 
where the subscripts of K's refer to the reaction step number, and 

a is an arbitrary constant. The value, bo = 0, is possible if either 
Koo = Q0or Ko7 =o. For the reaction to proceed in the forward 
direction to a significant extent, Ko9 cannot be equal to zero. If 
Ko9 =, then the equilibrium of step 5.27 is highly in favor of the 
formation of hydrogen gas. The latter possibility appears to be in 


(72) which showed that 


agreement with the previous reported results 
hydrogen only adsorbed on y-alumina in minute amounts. The model, 

N3, could also be obtained from the mechanistic Scheme D by assuming 

the decomposition of alkoxide, step (5.26), and the desorption of 
acetone, step (5.24), to be rate-controlling. Again, expressing bo 

in terms of the adsorption and desorption constants for this mechanism 
gives b, = K39 = 0. This implies that 2-propanol was not adsorbed, 
which is contrary to the experimental evidence. Also, acetone was not 
observed on the catalyst surface under reaction temperatures by 
infrared analysis. Since acetone is formed, it is presumed that the 
surface concentration is so small that desorption of acetone cannot 

be rate-controlling. This ruled out the possibility that mode] N3 

was obtained from the mechanistic Scheme D. Thus N3 should be regarded 
as obtained only from Scheme B with both steps (5.22) and (5.26) to be 
rate-controlling. Comparing model N5 to N8 it can be seen that if b, 


of NS is zero then N5 is reducible to N8. By fitting the existing 24 
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data points, to model N5, it was found that by = - 0.001 + 0.0016 
(95% confidence interval) and this implies that b3 indeed is insigni- 
ficant. Therefore, if N8 represents the data best, the kinetic 
model for dehydrogenation could be derived from Scheme D with the 
Surface reaction step being rate-controlling. 

The posterior probabilities of N3 and N8 were computed and 
were found to be 0.82 and 0.18, respectively. Although N3 showed a 
better fit than N8, its superiority is not conclusive. To discriminate 
further between these two rival models, the region of feasible experi- 
mentation was examined using the expected entropy change criterion (22), 
However, the scope of conditions available for predicting expected 
entropy changes was greatly restricted due to the inability of the 
IR analysis to detect hydrogen partial pressure. Theoretically, 
introducing excess hydrogen into the reaction system was desirable 
for model discrimination, but due to the inability of the IR to 
detect Hos it was decided to adjust the weight of the catalyst pellet 
to give desired dehydrogenation conversions. At such conditions, one 
obtains a hydrogen partial pressure which will give a maximum entropy 
change but without adding hydrogen gas into the feed stream. The 
partial pressure of hydrogen can thus be varied by adjusting the 
fractional conversion of the reactant and the partial pressure of 
hydrogen can still be assumed to be equal to the partial pressure of 
acetone,which can be accurately measured, because no side reactions 
were detected. 


The selected values of entropy change at various partial 
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Pressures of 2-propanol and hydrogen are shown in Table 5-10. Ac- 
cordingly, the data point at Pa = 20 cm Hg and PH ~ 0.4 was chosen 
by reducing the catalyst weight to 0.0586 gm from the 0.2747 gm used 
in the previous experiments. The new data point taken at the maximum 
expected entropy change yielded (-r,) = 0.07259, Pa = 19.972 cm Hg, 
and PH S Py = 0.397 cm Hg. This data point together with the pre- 
vious 24 data points provided enough information to rule out model 
N8. In fact the posterior probability for N8 was practically zero 
and that for N3 was unity. The results of the Calculation are summarized 
as follows: 
N3 
Parameters b, = 0.0045 b, = -0.00012 bs = 0.00734 by = 1.16624 
95% limits + 1.787 x 107° + 0.0043 + 0.000124 + 0.00887 
Variance of fit (0%) = 1.479 x 1079 


Posterior probability = 1] 


N8 
Parameters b, = 0.004 b, = 0.954 by = 0.0038 
95% limits +6x 10° + 0.034 + 0.0003 
Variance of fit (0°) =3x10° 


Posterior probability = 4.44 x 107!7 


The proposed models other than N3 and N8 were again fitted 
with the 25 data points. The results again myenicn that these models 
were inadequate. Thus, the successful model involved the dual dissimilar- 
sites mechanism (Scheme B) with adsorption of 2-propanol and decomposi- 


tion of alkoxide jointly rate-controlling. 
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Table 5-10 
Expected Entropy Change at Various Partial Pressure of 


2-Propanol (A) and Hydrogen (H.) for Dehydrogenation Reaction 


eee 
Expected Entropy 


Pas cm Hg PH? cm Hg Change - R 
Ai | J CE ES PD NMC ene 
0.5 ex 10° 
1.0 0.0005 
5 2.0 0.0006 
10 Gnee 1.02 x 1076 
10 mal 67 1.36 x 1077 
10 1.00 0.001 
20 0.4 0.0196 
20 te 0.0190 
20 0.6 0.0160 
20 1.0 6.0 x 1077 
30 0.3 0.016 
30 0.5 0.02 
30 0.7 0.02 
30 0.8 0.018 
30 1.0 0.0056 
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5.3.9 Activation Energy for the Dehydrogenation of 2-Propanol 


The most successful model, N3, obtained from the isothermal 


data at 343.2°C, is rewritten below according to Arrhenius expression: 


E 
] 
b1 0 exp (- RY? fe 
1+ b3 9 exp (~pF) m, * ba 0 exp(- RT "Py 
2 


In order to estimate b1 0 and Ey» low conversion data at various 
temperatures were obtained and are shown in Table 5-1]. The low value 
of P, together with the small value (0.00734) of bz at 343.2°C en- 
sured that b3 9 exp(-E./RT) Ca de in equation 5.31 is much smaller 


than unity. In this way equation 5.31 may be reduced to 5.32, 


E 
] 
Car (5.32) 


Fitting the experimental data with the above equation, it was found 
that the parameter estimation using Singh's method did not converge. 


It has been reported that similar difficulty may be due to the inherent 


121) 


functional form of the model under investigation| A particular 


model leading to such difficulty was 


1+ by Pa 


(-r,) =: (5.33) 
A b, + boPy 


Bf! 
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Table 5-11 
Data for Estimation of the Activation Energy 


of the Dehydrogenation Reaction 


Catalyst weight = 0.2853 gm 


Feed rate - 0.04088 g-mole of 2-propanol per hour* 


———— pre app re nenenepe nyse ere 


Run No. Reaction Reaction rate Pas cm Hg Py » cm Hg 
Temperature, °K (x 10°) 2 
ag 556 2.961 3.810 0.085 
45 567 4.110 oa7 Ga Othe 
46 577 5.462 3.750 0.147 
47 585 6.780 3.709 0.186 
48 593 8.33] 3.669 0.226 
49 599 9.659 3.633 0.262 
50 607 11.700 aeons, 0.318 


*No flow rate = 0.378 g-mole/hr 
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The form of equation (5.32) is similar to that of (5.33) and this 

is the reason for its lack of convergence. However, the Arrehnius 
plot, by assuming ba 0 exp (-E,/RT) = by = 1.166, shows a linear re- 
lationship with uniform distributed errors (Figure 5-15). Therefore 
Ey would be expected to be a small] value and by to be temperature in- 
sensitive. The apparent activation energy, Ey» is estimated to be 
30.3 k cal/g-mole from the slope in Figure 5-15. Again Singh's 
method for parameter estimation was applied. The constants, b, and 
bys which are functions of the equilibrium constants for each elementary 
step, are now assumed to be temperature insensitive for the range of 
temperatures studied, with the values, by = 0.00734 and by = 1.166, 
as obtained from the isothermal data. The data of Runs 44 to 50 in 
Table 5-11 together with the first 16 isothermal runs (Table 5-8) 
were correlated for equation (5.31). The following results were 
obtained: 

Parameters by = 2.148 x 108 —_E, = 30,090 

95% limits + 4.66 x 10” + 257 

Variance of fit = 2.34 x 1079 

An apparent activation energy of 30 k cal/g-mole is obtained which 

is in agreement with the Arrehnius plot and which is much higher than 
that of dehydration (14.6 k-cal/g-mole). It is thus expected that the 
dehydration reaction will be predominant on the y-alumina surface at 
lower temperatures. Since aluminum ions, which are responsible for 
the dehydrogenation reaction, are believed to be abundant on pure y- 


alumina, the dehydrogenation of 2-propanol on y-alumina is possible 


ot! 
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Figure 5-15. 


Temperature Dependency of the Rate of 


Dehydrogenation of 2 


-Propanol. 
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for temperatures above 500° (94) | However side reactions, such as 
the thermal decomposition of 2-propanol, become significant at high 
temperatures (2) so that little information can be obtained for the 
study of simultaneous dehydration and dehydrogenation of alcohols 


on pure y-alumina. 
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CHAPTER VI 
CONCLUSIONS 


6.1 Adsorption of Secondary Alcohols on y-Alumina and NaOH-doped 
y-Alumina 
Secondary alcohols are primarily dehydrated by hydrogen 


bonding to surface hydroxy] groups on y-alumina at high temperatures. 

At lower temperatures, an alkoxide surface Species also appeared and 

reacted at higher temperatures to form a surface carboxylate structure, 

which was stable up to 300°C. Contrary to other reports (679) the 

alkoxide did not lead to the formation of detectable amounts of ether. 
Doping of y-alumina with sodium hydroxide reduced the 

number of surface hydroxy] groups at 3785 and 3680 cm! , thus leading 

to a reduction in the hydrogen-bonded surface species formed from 

secondary alcohols. The amounts of the surface Lewis sites were found 

to be unaffected by increasing the catalyst NaOH content. Consequently, 

the surface alkoxide species, which were obtained by adsorption of 

the alcohol on Lewis acid sites, were found to be identical for al] 

of the catalyst samples investigated. An adsorbed ketone-like inter- 

mediate was observed analogous to the adsorbed secondary alcohol for 

the catalyst containing more than 1% of NaOH. With further heating 

to 300°C the ketone-like intermediate formed a stable carboxylate 

surface structure accompanied by the release of methane from secondary 


alcohols, 2-propanol and 2-butanol. 
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6.2 Catalytic Properties of y-Alumina Towards Decomposition of 
Secondary Alcohols 


The Bronsted-type acidity exhibited by y-alumina surfaces 
is related to the hydroxy] bands appearing at 3785, 3720, and 3680 cm”! 
Roughly 10 to 15% of the Bronsted-acid sites are stronger acidic sites 
than the remaining 85 to 90% sites. The Lewis-type acidity behavior 
is related to incompletely coordinated aluminum ions. The 3720 cm”! 
hydroxyl groups and the Lewis-acid sites appear to be less accessible 
than the 3785 and 3680 cm”! hydroxy] bands to doping agents such as NaOH 
or to chemical reactants. The 3785 cm”! hydroxyl groups are more 
catalytically active and acidic than the 3680 cm”! hydroxy] groups, con- 
trary to the model of the surface of y-alumina developed by Peri. 

The structures, = Al-0-Na, seem to be formed when alumina 
is doped with NaOH. The number of aluminum ion sites remains un- 
affected by increased NaOH content but their accessability increases, 
hence the catalyst selectivity may be continuously modified from de- 
hydration alone, to dehydration plus dehydrogenation, to dehydrogena- 
tion alone by proper doping procedures. 

The selectivity of various catalysts towards simultaneous 
dehydration and dehydrogenation reactions could be related to the 
number and accessability of Bronsted-acid and Lewis-acid sites on the 
surface. It also appears that the catalytic activity of the weak 
Bronsted-acid sites and the Lewis-acid sites do not differ substantially, 
but the Bronsted- acid sites (hydroxyl groups) are more accessible 


than the Lewis-acid sites (aluminum ions) because of their greater 
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proximity to the exterior surface planes. 

The mechanisms for dehydration and dehydrogenation of 3- 
pentanol on pure y-alumina or NaQH-doped y-aluminas involve hydrogen 
bonding of the alcohol to the Bronsted-acid sites (via the oxygen 
in the alcohol) or formation of a surface alkoxide species with the 
Lewis-acid sites, respectively. NaOH-doping reduces the number of 
Bronsted-acid sites thus suppressing the dehydrating activity of 


alumina. 


6.3 Kinetics of Dehydration and Dehydrogenation of 2-Propanol on 
y-Alumina Catalysts 


The recirculation reaction system for the study of simul- 
taneous kinetics and mechanism of 2-propanol decomposition using 
infrared spectroscopy was found to be versatile, simple and accurate. 
Since mass transfer limitations for the pressed catalysts were shown 
to be negligible, the pressed wafers were suitable for the kinetic 
studies. 

On the basis of the surface spectra obtained under the re- 
action conditions, together with the results from the adsorption studies, 
the dehydration and dehydrogenation of 2-propanol appear to involve 
two types of active sites for each reaction: surface hydroxyl and 
surface oxide for dehydration, and surface aluminum ion and surface 
oxide for dehydrogenation. The mechanism proposed from the adsorp- 
tion studies for the decomposition of 2-propanol on pure y-alumina 
and on 8% NaOH-doped y-alumina are in agreement with the results from 


the kinetic data obtained. The kinetic data were correlated for the 
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proposed mechanisms by a Baysian approach to model discrimination. 
The reaction rate expressions that fit the data best are Gury) = 

by Pa/(1+b5Pq) (1+b5P),) for the dehydration reaction and (eta) = 
Pal (1+b, iat, + aie for the dehydrogenation reaction. 
The apparent activation energies for the dehydration and dehydro- 


genation reactions are 14.6 and 30.1 k cal/g mole, respectively. 
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CHAPTER VII 
RECOMMENDATIONS 


(1) The design of the infrared cells in the present study was 
restricted by the width of the sample compartment of the infrared 
spectrophotometer (11 cm). In order to determine low concentration 
components accurately, a long infrared cell may be required. This 
can be achieved by a combination of three mirrors, similar to those 
used in commercially available multi-path gas cells, as shown in 
Figure 7-1. The cells should be enclosed in an isolated system which 
has been purged with dry nitrogen gas. This will Prevent the spectral 
interference due to atmospheric water vapor and carbon dioxide. 

(2) Figure 7-2 illustrates a proposed scheme whereby the flexi- 
bility and reliability of a high-quality infrared spectrophotometer 
can be used as an “infrared analyzer". The experimental parameters, 
such as temperature or reactant flow-rate, are under direct digital 
control during either a planned sequence or computer-determined se- 
quence of experimental conditions. The IR spectrophotometer monitors 
a particular band continuously via the computer interface or alter- 
natively, it is instructed by the computer to commence a scan over a 
pre-set range. The latter scan may be repeated as often as desired 
under computer program control. Both qualitative and quantitative 


analyses are feasible with such a scheme. 
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Figure 7-1. Schematic Diagram for the Use of 
Long IR Cell in Kinetic Studies. 
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INFRARED "ANALYZER" CONST. TEMP. BATH 
PRODUCT 
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Figure 7-2. Online Monitoring and Control of 
Chemical Reaction/Kinetics Studies. 
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(3) Simultaneous dehydration and dehydrogenation reactions were 
observed for the catalysts containing one to eight per cent NaOH. It 
is further recommended that additional experiments need to be per- 
formed on these catalysts in order to determine the effects of NaOH 
content in y-alumina on the decomposition of alcohols. By assuming 
the same mechanisms obtained from dehydration and dehydrogenation 
Studies, the rate expressions for simultaneous dehydration and de- 
hydrogenation of 2-propanol on NaOH-doped y-alumina can be derived. 
They are given as follows: 


(1) Generalized mechanism for dehydration reaction: 


Ath = AL, (1) 
AL, + hy > PL, + We, (2) 
PR) = P+ a (3) 
Wo = Wt Lo (4) 


(2) Generalized mechanism for dehydrogenation reaction: 


a (5) 


AL, + 2, = XL, + Hi. (6) 
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XL, + hy > K2q + HL, (7) 
Ht, = He + 22, (8) 
Kt, = K+ he (9) 


where has Lo and 2, are surface hydroxyl, surface oxide and surface 
aluminum jon, respectively. 

(a) The rate expression for the dehydration reaction is derived 
by assuming that the surface reaction, step (3), is rate-controlling. 


Then 
Le pM eS ily 0 
rate of dehydration = r, = k, C,. Cg 
eC 2 Ax, 2 


A further assumption is that the total number of sites of type 2, Los 


is constant and L, = constant = Lig f(Na), where Ly 9 is number of OH 


groups on pure y-alumina and f(Na) is some function of NaOH content 


which can be determined experimentally. Thus 


+ 


Ly = L100 f(Na) = C2, + “ag, PL, 


Pp 
= C2,(1 + K)Py + ree 
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Let the coefficient of Pp + 0, then the rate of dehydration of 2- 


propanol on NaOH-doped y-alumina can be expressed as 


be derived in a way similar to that for the dehydration. 
that the adsorption of 2-propanol (step (5)) and the surface reaction 


(step (7)), are both rate-controlling, the rate of dehydrogenation 


: b, f(Na) Pr 


r 
2 
(Tb Pe) (+ bP by Pu) 


(b) The rate expression for the dehydrogenation reaction can 


can be written as: 


beP A 


rate of dehydrogenation = 


Leb aPw thay Pu, eas tPonatUp pai + be 


Thus the rate of 2-propanol] conversion is the sum of the rate of 


dehydration and rate of dehydrogenation: 


By assuming 
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LIST OF NOMENCLATURE 





The nomenclature in use throughout this thesis is given in 


the following list. Symbols and constants are also defined where they 


first occur. Special symbols which are defined and used in the Appendix 


are not included here. 


English Letter Symbols 


Alcohol 

Absorbance at wavelength i 

Arbitrary constant 

Vector of unknown parameters corresponding to the qth model 
Total number of active sites which can be Occupied by jn 
component 

Arbitrary constant of jth 

Concentration (g-mole/2z) 

Variance covariance matrix of the parameters in a kinetic mode] 
Effective diffusivity 

Activation energy 

Entropy as defined in equation (2.32) 

Apparent activation energy (K cal/g-mole) 

True activation energy (K cal/g-mole) 


th 


Error term in the ji kinetic model 


Bel 
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Molal flow rate of the reactant (g-mole/hr) 

Enthalpy change for adsorption of component i 

Intensity of the transmitted beam 

Intensity of the incident beam 

Ketone 

Equilibrium constant for the reaction step i, or adsorption 


constant for jth 


species 

Reaction rate constant for the reaction step i 

Total number of active sites of type i on the catalyst 
surface 

External surface area-to-volume ratio for the catalyst wafer 


Length of the IR cell 


Surface active sites of type i 


Assigned model number i for the dehydration reaction 


Assigned model number i for the dehydrogenation reaction 
Total number of data point obtained 


Propylene 


h 


Probability that the it model form is correct, or partial 


pressure of jth component 


h 


Partial pressure of it component in the product mixture 


Partial pressure of reactant in the feed 


h component measured in the reference 


Partial pressure of iv 
cell 
Ideal gas law constant 


Expected entropy change 
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as 


Rate of reaction; g-mole of reactant converted per hour 


per gm of catalyst 

Rate of reaction for the reaction step i 
Initial rate of reaction 

Minimum sum of errors squared; equals 


min 


n 
oh oe itepsaruedaa J Vcaled, i 


2 


Absolute temperature (°K) 
Water 
Weight of catalyst (gm) 


te neLicunodel 


Design matrix for j 
Surface alkoxide species 

Independent variable 

Fraction of reactant converted in to product 


Dependent variable 
Greek Letter Symbols 


Extinction coefficient; equals 1/C log (Ip/1) 
Error of variance of a single measurement 
The unbiased estimate of error variance 


h 


Proposed i* kinetic model 


Thiele diffusion modulus 
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A. ON-LINE OPERATION OF A PE621 - IBM 1800 SYSTEM 
ee a eet = ADM 16U0 SYSTEM 


Analytical instruments, such as infrared spectrophotometers 
equipped with digital encoders, can generate a great deal more data, 
and of higher precision, than is normally available with conventional 
instruments. This additional data and precision is easily handled 
by a real-time digital computer, such as an IBM 1800, and sophisticated 
data processing methods, including signal filtering and averaging, can 
be used in generating quantitative descriptions of spectral bands 
which characterize molecular species. Moreover the manner of process- 
ing each set of data is left for the user to decide after data acquisi- 
tion; this gives a wide range of flexibility. Several references in 
which such concepts are developed are already available‘ !22>123) 

A block diagram of the infrared Spectrophotometer inter- 
faced to the computer is shown in Figure A-1. The multi-application 
capability of the system is also indicated but emphasis is placed 
upon the participating role of the IR unit. The diagram shows the 
flow of data originating as the signals, Y' and X', through to their 
presentation in several forms of digitized output, Y and X. 

The advantages of this type of computerized application 
arise from the elimination of conventional input-output contro] 
systems, greater quantity, and often superior quality of data taken 
in the same scan time, and benefits occruing from rapid, detailed 


processing of the data. The following section will describe in de- 
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A-3 
tail the three major equipment items comprising the PE621 - IBM 1800 


system, and their methods of communicating information or instructions. 


A-1 System Components 
A-1.1 PE621 Infrared Spectrophotometer 


This instrument amd its capabilities are well known and so, 
detailed comments on its design are unnecessary in this thesis. The 
conventional output signals are percent transmittance (Y' = source 
signal, Y = digitized signal) and wave-length (X' = source signal, 

X = digitized signal). The PE621 was purchased with Perkin-Elmer 
X' and Y' encoders installed. The two encoders are one-brush abso- 
lute position types and their shaft positions are tracked continuous ly 


124) | Freedom from gearing, backlash and 


by flip-flop registers: 
possibly servo errors, particularly those inherent in systems taking 
data from the recorder servo unit, are eliminated through the use 
of encoders. The Y signal in three digits is obtained from the comb 
servosystem. The X signal, resolved to a tenth of a wave-number, 
is obtained from a 100-count encoder mounted to the wave-length shaft 
(monochromator). An up-down counter (three digits) in the interface 
unit counts the number of revolutions of the 100-count encoder (two 
digits), thereby giving a total wave number of five digits. Since 
subsequent transmission of the eeiotiea Signals to the computer is in 
BCD code, this system permits the monitoring of Y-X signals with high 
resolution, the resolution being limited only by that of the encoders. 
A-1.2 Encoder Readout-Computer Interface 

The Perkin-Elmer Model 8RLS encoder readout-computer inter- 


face includes an operating mode for direct transmission of Y-X data 
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A-4 
to a digital computer. The signals are received continuously from 
the IR in encoder switching code and are then translated to BCD and 
stored in an encoder register. The BCD data may be decimalized for 
visual Nixie display or alternatively, transmitted to a data-logging 
device such as an external paper tape punch coupler. The data in the 
encoder register is strobed into a buffer storage register which 
freezes the information. A translator converts the buffer register 
Output into 8-4-2-1 BCD output in positive logic (logical one = +8 
volts, zero = - 12 volts) and transmits it to a rear panel connector 
for direct hookup to the IBM 1800 computer. 

The encoder readout has a maximum capacity of 32 bits 

(eight decimal digits) and can accept up to 300 bits/sec (roughly, 

ten Y-X points per second). The data encoding rate is operated asyn- 
chronously with the computer since synchronized real-time monitoring 
of data points is not desirable with the hardware-control led operation 
of the spectrophotometer. The wave number interval which determines 
the data sampling rate is specified by operator selection of multiples 
of the minimum increment (0.1 cm!) available for the X-axis. The 
sampling rate is generated by a data trigger which strobes the Signals 
at the specified intervals of X from the encoder register into the 
buffer register for parallel transmission of all bits in the 8-digit 
data word to the computer. Each time that a Y-X point has been trans- 
mitted from the encoder register to the buffer register, and thus to 
the computer, the encoder register is left in a "reset" position. The 


computer acknowledges receipt of the data point by sending a voltage 
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A-5 
pulse to the interface unit (one pair of wire needed), which returns 


the encoder register to "set" position. Ifa second point should be 
Strobed into the buffer register without the preceeding computer 
acknowledgement, the "reset" position indicates loss of data and 
initiates an electronic contact operation alarm. 

In the kinetics studies, the IR spectrophotometer monitors 
a particular frequency as a function of time via the encoder read- 
out - computer interface. Since the data trigger is not activated 
without varying the wave length, X, a pulse generator was installed 
to tramsit data from the encoder register to the buffer register. 
The sampling rate is controlled by the rate of pulse generated. 
The other functions are identical to the normal operation mode. 

The computer and IR-interface are separated by about 
125 ft. of the communication cable (two 18-pair 22-gauge individually 
shielded twisted-pair cables). The cables are totally shielded as 
well and are run through conduits for mechanical Protection, pro- 
viding rather excessive shielding protection. The BCD signal, the 
Start-up communication, and the computer acknowledgement pulse require 
34 pairs leaving two pairs in reserve. It is interesting to note 
that the essentially noise-free 32 pairs of cable could be replaced 
by at most two pairs if the analog Signals, Y-X, were to be trans- 
mitted. Noise pick-up in transmitted Y-xX Signals is more likely and 
could reduce the data accuracy considerably. 

The direct generation of digital values and the transmission 


of the digital BCD signals should be essentially noise-free. The 
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quality of data input to the computer therefore is directly dependent 
on the accuracy of the encoder readout and the basic IR components. 
A-1.3 IBM 1800 System 

The hardware configuration most pertinent to the on-line 
operation of the IR unit is shown in Figure A-2. The equipment de- 
signated by the boxes with slightly thicker margins has been involved 
in the current studies. Access to such an extensive facility with its 
large core and disk storage facilities, full range of conventional 
input-output peripheral hardware, substantial computing power, and 
supporting system analysts greatly simplified the task of developing 
a powerful and convenient-to-use system. 

The system was initially operating with an IR data acquisition 
program permanently resident in core when the IBM 1800 was operated 
under Time-Sharing Executive (TSX) System. After the computer was 
converted to operate under Multi-Programming Executive (MPX) System 
in 1970, the IR data acquisition program was stored on disk and alter- 
nated with the Department's Gas Chromatograph Monitoring Program in 
core load area three. The program responds to interrupts generated 
by the interface unit and reads in 32 bits of data. The 32 BCD input 
data bits are converted to binary integer numbers and stored into 
two 16-bit words of a core buffer. When the 320-word core buffer is 
filled, it is transferred to one sector of disk storage. In view that 
this is a time-shared operation, the disk could be busy completing a 
disk operation for another program. Thus provision had to be made 


for queuing time and disk-seek operations as well as for the actual 
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10 ms. of write time. Four possible methods (125) 


of insuring a disk 
write time compatible with the IR data acquisition program limitations 
are dedication of a disk, dual internal buffering, look ahead method, 
and secondary overflow buffering. Only the last method makes efficient 
use of both core and disk. Since minimized core requirements were 
a design specification for the IR data acquisition program, it was 
decided to use the secondary overflow buffering technique. 

The provision of a secondary overflow buffer would allow the 
IR routine to accept a limited number of data samples while the main 
buffer is being written to disk. When the main buffer is freed, data 
stored temporarily in the secondary buffer would be transferred into 
the main buffer. For the maximum data rate (100 ms. sampling interval), 
a 40-word secondary buffer would allow 2 seconds for disk write 
operations. The size of the auxiliary buffer can be changed by an 
equate card at system generation time. 

The disk data files were stored in a form accessible by 
FORTRAN coded programs so that the users could write their own 
data processing programs. Since the five-digit wave number could be 
as large as 40,000 and since the largest number that can be represented 
by the 16-bit data word of the IBM 1800 is 32,767 it was necessary to 
use either a two word integer or a special data code. To conserve core, 
it was initially decided to store the four least Significant digits 
of the wave number in a one word integer format with the associated 
transmittance in the adjacent word. Each time that the fifth (highest 


order) digit changes, it is recorded by putting a special entry in the 
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file with a transmittance of "-1" (a physical impossibility) and a 
wave number equal to the fifth digit. The users' program can then 
access the file and handle it as one-word integers, or convert it to 
two-word integers or floating-point numbers. Negative values for 
percent transmittance are also used to "flag" special codes to indi- 
cate run termination or alarm conditions. They are shown in Table 


A-1. 


A-2 Operation of the System 


The flow diagrams in Figures A-3 to A-5 describe in simpli- 
fied terms the role of the computer during the startup, operation and 
shutdown of program-controlled data acquisition. The IR operator in 
the remote laboratory actuates a "request to computer" electronic 
Signal when an IR scan is to be performed. 

In the "start" procedure, upon assuming the "ready" state, 
the computer informs its operator via print message that it has 
initialized data acquisition and also informs the IR operator of 
its "ready" state. The latter action is accomplished via a light on 
the panel of the interface unit. It then returns to its time-sharing 
service role until the IR operator, at his convenience, starts the 
hardware-controlled IR scan. 

When reading data, the computer files all of the raw data 
before any data processing is commenced. When the IR operator signals 
by turning the on-line switch to the "off" position that the scan has 


been completed, the computer again advises its operator that the job 
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Table A-1 


The Usage of the Special Indicators in Data File 


me eg 


Transmittance Wave Number Explanation 


-2 Undefined This record indicates the end 
Of valid data for the most 
recent run and that the run 
was terminated normally by 
either the operator or the 
preset range switch of the IR 
spectrophotometer. 


-3 Undefined This record indicates the end 
of valid data for the most 
recent run and that the run 
was terminated due to either 
an external or internal data 
over flow. 
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Figure A-3,. Start Procedure. 


Figure A-5. Stop Procedure. 
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A-12 
has been completed and it then summarizes the system statistics. 
In the above procedures, all logging and error messages 
issued by the IR routine are typed on the log typewriter in a standard 
format for reasons of simplicity and core economy. The message format 


is: 


IR CLOCK CLOCK XXXXX ‘els X 


ioave File used 
Number of data points 
(message 02 only) 
elapsed time (message 02 only) 
time of message 


message identifier 
IR routine identifier (message origin) 
The messages are further identified by color, with normal logging 


messages being typed in black and error messages typed in red. The 


present system messages are defined as follows: 
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Message Identifier 
01 


02 


03 


05 
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Explanation 
Log-on message (indicates that the IR 
system was initialized for a data 
acquisition run) 
Log-off message (indicates that data 
acquisition was terminated normally) 
Internal data overflow (indicates that 
data acquisition was terminated be- 
cause of the failure of a disk write 
operation to be completed before the 
maximum allowable number of overflow 
samples were accumulated) 
External data overflow (indicates that 
data acquisition was terminated due to 
a signal from the IR spectrophotometer 
Signifying that a data sample was not 
acknowledged by the system before the 
next sample was ready) 
Data file overflow (indicates that data 
acquisition was terminated due to the 
fact that the maximum allowable number 
of samples, as specified by the size 


of the defined file, were accumulated) 
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On the initiation of a run, the IR program reads the last 


four positions of the "Special Digital Input Switches" to determine 


at what wave number to start collecting data and what file to put the 


data into. 


Digit Switch 


3 
OFF 


The significance of these four switches is as follows: 


14 
OFF 
OFF 
ON 


ON 


Data Acquisition to Start 


Le 

OFF 09999 or 
ON 19999 or 
OFF 29999 or 
ON 39999 or 
OFF 49999 or 


Data stored in file "IRDTO" 
Data stored in file "IRDTI" 


at Wave Number 


Since October, 1970, the IBM 1800 computer was converted 


from the TSX operating system to the MPX operating system. The IR 


data acquisition program was thus revised according to the require- 


ments of the MPX System. 


However, the system role remains much the 


Same. Figure A-6 shows the details of the operational procedures 


for the present data acquisition system. 


A-3 Data Processing 


The treatment of infrared spectroscopic data which has been 
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FINISH 











CORRECT 
FAULT IF 
POSSIBLE 
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I.R. RUN 













DETERMINE 
FAULT IF 
POSSIBLE 


*STATUS BOARD FLAG 


Figure A-6. IR On-line Data Acquisition Procedures 
(MPX System) . 
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A-16 
obtained with the system described herein involves processing up to 
38,000 Y-X points if the full range of 3800 cm”! is scanned. These 
points would occupy 76,000 words of memory, or more if several com- 
plete spectral scans were to be retained. Only with access to a time 
Sharing computer system, with relatively large disk storage and rapid 
computing capability, can the experimenter obtain reduction of such 
large amounts of data in a speedy manner. 

Data processing routines such as data filtering, peak 
picking and other numerical methods can be used on the data file, and 
where advantageous, the benefits from varying the computational para- 
meters can be exploited. The raw data and/or the resulting smoothed 
data can be stored on disk cartridges or punched cards, printed on the 
line printer or typewriters, or displayed using the digital plotter or 
the oscilloscope. 

For the sake of efficiency, for both the user and the com- 
puter, the development of adequate data reduction programing must 
involve a careful evaluation of the user's needs. Smoothing of the 
data is advantageous in research applications such as this Study if 
spectral characteristics may be ambiguously interpreted because of 
high random-noise levels in the percent transmittance or Y data, j.e., 
random error superimposed upon the significant output signal. If the 
Signal-to-noise ratio is very high, routine smoothing of data may be 
wasteful through unnecessary use of computational time in a time- 
sharing system. If selection of peaks, i.e., the Significant spectral 


bands which characterize the sample molecule, is required, then some 
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A-17 
level of data smoothing may be advantageous so that those changes in 
slope which define significant peaks can be more readily detected. 
A considerable amount of work on such problems has already been re- 


ported by Savitzky (127) , 


Much of this thinking and programming! !26) 
has been incorporated into the data reduction methods employed in 
this work. The methods involve least Squares fitting and assume that 
the error is associated with the Y Signal (three significant digits) 
and not with the X signal (five significant digits). 

Figure A-7 presents a simplified data-reduction flow dia- 
gram which relates the major steps involved in the data smoothing and 
peak-picking routines. The data acquisition involves multiple-user, 
on-line operation and so, the data reduction and Output functions 
are not commenced until the scan is completed. Although not shown 
in Figure A-7, the data points are read from the disk file, sector by 
sector. They are subjected to least-squares smoothing and a slope 
calculation except where the transmittance reads 3% or less. The 
Slope is calculated at each point using a quadratic function fitted 
to the smoothed data. When the sign of the Slope inverts, defining 
Passage across an absorption band, a 15-point cubic equation is used 
to determine the ordinate and abscissa at the peak minimum. The raw 
data is completely treated, the smoothed data replaces it in the disk 
file. At present, if a plot of the raw data is required, this out- 
put must proceed the smoothing since the latter is destructive to the 
original memory file. Upon completion of these steps, the Y-X plotter 
programs calculate appropriate scales and automatically plot the 


function. The evaluation of the effect of computational parameters 
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READ X, Y 
DISK 


CALL SMOOTH SUBR, 


9 - pt. Smoothing, 
least squares 


CALC. (Y, X) AT dY/dx = 0 
15 - pt. cubic equation 






PLOT Y-X 
PRINT PEAK 


Figure A-7. Data Reduction Flow Diagram. 
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A-19 
upon smoothing and peak-picking, e.g., different levels of filtering, 
was usually done by off-line processing. 

Figure A-8 illustrates a comparison on an expanded scale 
between some raw and smoothed data for a polystyrene film in an 
ATR cell. The figure was actually plotted by programmed control, in 
the manner shown, to conform to a preset 18-inch ordinate scale. This 
range of the X-axis was chosen because of the relatively high noise 
level visible on the PE621 drum recorder chart, a line of thickness 
roughly 1% transmittance. The four curves on Figure A-8 show the raw 
data, taken at highest sampling rate of every 0.1 cm”! interval, and 
the influence of smoothing these raw data. Some differences in data 
filtering are indicated from using all of the data, half or one-quarter 
of the data. 

In the smoothing procedure, it would be expected that the 
nine-point filter band (or “smoothing gate") is effective in Suppress- 
ing random noise if the width of the gate is considerably larger than 
the period of the noise frequency to be removed. The application of 
smoothing to the complete file of data sample involved a 0.9 cm”! gate 
width (as shown by the second plot of Figure A-8); this did not prove 
to be effective. Since the wider spacing of the data points produced 
by using every Nth point in the file generates a wider effective band 
with the nine-point filter, the smoothing became more noticeable as 
can be seen on plots 3 and 4. Because less data is employed when the 
gate is increased, the quality or significance of the resulting plot 


could also be reduced. Replication of a run, and averaging to eliminate 
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random noise, is essential to establish the best smoothed function 
which could be used as a criterion of smoothing effectiveness. 

An alternative explanation, originating with Savitzky(!22) 
relates the lack of smoothing which resulted from use of all of the 
data to the difference in magnitudes of time constants between that 
of the signal-generating instrument and that associated with the span 
of smoothing-gate. Until the time constant of the latter exceeds that 
of the instrument, it would be unrealistic to expect beneficial smoothing. 
This may be the case for the ATR signal until the time constant associ- 


ated with a 1.8 cm” | gate span was attained. This span was obtained 


with the nine-point gate using each second Y-X point. 


A.4 Research Applications 


The research programs, which utilized on-line capabilities, 
have presented two types of data interpretation problems. The first 
involves a comparison of infrared spectra of catalyst surface "before 
and after" the adsorption experiments to determine the spectra of ad- 
sorbed species. The second problem involves obtaining kinetic data 
from monitoring the characteristic peak intensities of the reactant 
and product. Both of these studies have benefitted from computerized 
treatment of data: the former, by more sensitive diagnostic analysis 
resulting from increased numbers of data points and by elimination 
of tedious point-by-point visual measurements, and the latter by 
replacement of manual calculations with computer operations thus per- 
mitting a direct output of kinetic data when the calibration functions 
were stored in the disk files. This feature enables the operator 
in attendance to change experiments or conditions quickly when non- 


routine unanticipated results are encountered. 
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A-22 
A-4.1 Measurement of Differential Spectra on a Single Catalyst 
Figure A-9 illustrates some of the spectral characteristics 
which are encountered during infrared studies of adsorbed molecules 
and surface reactions on solid catalytic substances using the Eischens 


technique! !27) . 


By plotting the difference in percent transmittance 
between the two spectra at the same frequency, taken before and after 
the experiment and noting also the mid-points for all peaks, the ex- 
tensive data are reduced to that which appears to be Significant. 
Often, the full spectral range spanning 3800 cm”! must be 
Studied for the "before" as well as the "after" experiments so that 
the significant differences between the two sets of spectra can be 
plotted. The differences between the two sets of spectra may be 
plotted in the form of the "delta-plot" shown in Figure A-9. This 
being a tedious but simple task, it is ideally suited for computer 
Processing, especially when smoothing is advantageous. The delta- 
plot shows only the differences in transmittance at each frequency 
and thus classifies the changes according to whether they are 
positive or negative peaks. Positive peaks imply an increase in an 
existing band or the appearance of a new band, except when a frequency 
shift occurs. The negative peaks indicate either disappearance or a 
decrease in transmittance for an existing band. The use of the delta- 
plot together with a list of mid-peak frequencies for both before and 
after spectral scans, all provided by computer output, would present 


savings in calculation time and an improved reliability in the results. 
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Figure A-9. Comparison of IR Spectra "Before and 
After" the Experiment. 
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A-4.2 Infrared Monitoring of a Reaction system 
In this case the IR spectrophotometer monitors a particular 
band continuously via the encoder readout-computer interface. Since 
large amounts of data were generated, it was able to determine the 
reaction steady state by fitting the data with a least squares equation 
over several periods of time and observed the change of slope of the 
concentration curve. Both qualitative and quantitative analyses are 
feasible with such a scheme. The computer programs used for the data 
processing are listed as follows: 
1. IRIST: Remove special indicators in data file and smooth 
the raw data 
2. PCNTR: Determine peak centre by 15-point cubic equation 
3. PKFND: Calculate slope by 7-point quadratic equation 
4. DETRN: Determine if transmittance < 3%, define range for 
peak picking routines 
5. BKGRD: Remove special indicators in data file, store data 
in data file of either IRDT1 or IRDTO for plotting 
programs . 
6. IRPLT: Plot data from file IRDTO or IRDT1 
7. DELTA: Compute the difference in transmittance for the 
| data between IRDTO and IRDT1, plot the resulted data 
8. KINET: Obtain kinetic data by following the procedures 


described in Section 4.2.3. 
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A= 26 
C MAINLINE IR1ST 
MAINLINE IR1ST 
MAIN LINE OF IR PROGRAM 


SHORTEN RECORD WHEN WAVENUMBER PASSES MULTIPLE OF 
100060 CM$1 . 
SETS CONDITIONS FOR FIRST AND LAST RECORDS 


DIMENSION DAT(80) »sWAVE(80) »¥KDATA(210) sIWAV(210) » 
INDATA (210) sLDATA(210) sLWAV(210) sPP(10) »83AT(80) 
2sMWAV(25) sMDATA(25) 

COMMON ITIL sJUWAVs CONST 

DEFINE FILE 1015 93209UsKI) 92(3093209UsIK) 93 
1(3093209UsKO) 

KN=0 

CALL TSTOP 

IY=1 

MEND=0 

IN=1 

IP=1] 

IQ@=1] 

DO 695 JK=1915 

READ(1'JK) (IWAV(LQ) sKDATA(LQ) sLQ=15160) 

DO 667 JJ=19160 

LL=160 

ITER=0 


THIS INDICATES WITH END OF DATA OR A MULTIPLE OF 


100060 CM$1 


IF(KDATA(JJ)+1) 66896699667 
668 IF(4J=-1)745597455670 
670 LL=JU=1 

ITER=1 

GO TO 798 

669 CONST=FLOAT(IWAV(UU))*1000e0 
WRITE(69671) IWAV(Uu) 

671 FORMAT(1H1s21HWAVENUMBER HAS PASSED®I591l0HO00e0 CM=1) 
DO 672 LJU=JUs160 
KDATA(LU)=KDATA(LU+1) 
IWAV( LU) =IWAV(LU+1) 

672 CONTINUE 
LL=LL=1 
GO TO 798 

667 CONTINUE 
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THIS IS A 


MAINLINE IR1ST eee (CONT'D) 


N INSERTION 


A= 26 


STORE UNSMOOTHED DATA IN FILE FOR PLOTTING PROGRAM 


DO 799 LM=lsLL 


DAT(IN) =F 
WAVECIN)= 
IF (IN=80) 
WRITE(2'*I 
IP=IP+1 
IN=0 

IF (ITER) 
IF (LM=LL) 
DAT(IN+1) 
WAVE(IN+1 
KN=IN+1 
GO TO 796 
IN=IN+1 
CONTINUE 
GO TO 123 
WRITE(2"1 
IF (JK=1) 
LL=139 


SMOOTHING 
DEFINE CO 


IQP=LL+20 
IQQ=LL+16 
LF (JK=1) 
IBA==-1 
IBB=5 
1BD=1 
IBE=1 

DO 12 Iwe 
NDATA( IW) 
GO TO 677 
IBA=11 
IBB=17 
LA=LL 

IF (JK=1)6 
IBG=20 

GO TO 627 
IBG=I1BD=1 
DO 624 K= 
LDATA(K)= 


LOAT(KDATA(LM) ) 

(FLOAT CIWAV(LM))) +CONST *10¢0 
78897919791 

P) (DAT(IU) sWAVE(TIU) sIU=1580) 


77597759789 
77596019601 

=-1.0 
Y=EFLOAT(IP#804+IN) 


P)(DAT(UN) sWAVE( UN) sUN=19KN) 
67436749673 


LOOP IS INTRODUCED 


NTROL CONSTANTS 


67596759676 


194 
=KDATA( IW) 


2596259626 


1sIBG 
MDATA(K) 


LWAV(K)=MWAV(K) 


CONTINUE 
IF(UK=1) 
LA=LL+20 


68196819679 




























aS =A ; 
(G'THOD) «se T2Lnl BWLINIAM De "as 

WOITHIeMT NA Sl SrHT > 

MASOORA DAITTOIA HOA 3419 WI ATA pent sey aAnoTe _ * 


JaeLeMs eer 00 gen 

{ (MJVATAGX) TAOUT= (METAR >% 

Oe0L* T2vVOD+ CC (MJ)VAWT )TACIF)2 (KT PSVAW 3 
FOF eLevaeser(Osev1Al | 2. 

10601 A Ee (GI'S)atiaw rey 
5. appr ti: , 


- OSKL. 
eats ayreatt (satel 6e8t i 
{0defOdeeTTCLJaMsr AT eer 
Osle=(lL+n1)TAd ne 
(AT +O8RGT I TAO IA RC L+AT VAN I 


ps tpn — 
oe eer : 
eSi of OD 


. (vive Pee (AL) BWAWe (WL TAQ) EST! S) STEAM eet <5 
ETSe A TaePTO( LeALy AT est 
sin eget oT 

(2 23..° 1 7RAS 

a5UQORTAT 21 400) aATHTOOMe- 

2THAT2NOD ETHER pt a 
OS+iJ=901 eta 
olejJ2001 

5 GCaR TQ SGe: Oe Stone 
fosagl era 





eae a3 SRP 
en aicatinee Pa 

- eit. bad hs ane eT 
free AEBS 4 8 


eet slide. Papin 


4 2% re i535 ar 


. % 


ll 


200 


741 


740 


G MAINLINE IR1ST 


DO 678 IXZ=1lsLA 
KK=IBD+1XZ=1 
LDATA(KK) =KDATA(IXZ) 
LWAV(KK)=IWAV(IXZ) 


INITIALIAE SMOOTHING 


DO 10 KZ=299 

JZ=KZ+IBA 
PP(KZ)=FLOAT(LDATA(UZ)) 
IF(JK"1) 66496649665 
IBC=LL+16 

GO TO 867 

IBC=LL+IBD=5 

CONTINUE 


MAIN SOOTHING LOOP 


DO 200 IA=IBBsIBC 
JA=IA+4 

DO 11 KY=198 

KA=KY+1 

PP(KY)=PP(KA) 
PP(9)=FLOAT(LDATA(UA)) 
SUM1=59 e0*PP(5)/231.0 


eee (CONT'D) 


SUM2=540e0#(PP(4)+PP(6))/23160 


SUM3=396%(PP(3)+PP(7))/231e 


SUM4=1460*(PP(2)+PP(8))/23140 
SUM5=214e0*(PP(1)+PP(9))/231¢40 


SUM=SUM1+SUM2+SUM3+SUM4=SUM5 


NDATACIA)=IFIX( SUM) 
CONTINUE 


CALL DETRN(NDATAsLWAVsLL 9UK9IsITsMEND) 


DO 202 IL=IBEsIBC 
BAT(IY)=FLOAT(NDATA(IL)) 
BAVE(TY)=FLOAT(LWAV(IL)) 
IF(IY=80) 74097419741 


+CONST *10¢0 


WRITE(3'IQ) (BAT(1IZ) sBAVE(IZ) s1Z=1980) 


IQ=IQ+1 

TY=0 

IFCITER=1) 74397423742 
BAT(1Y+1)==1.0 


BAVE(IY+1)=FLOAT(IQ*80+IY) 


IXX=1Y+1 


WRITE(3'IQ) (BAT(K) sBAVE(K) sK=19IXxX) 


GO TO 745 
TY=IY+1 
CONTINUE 


Arm ZF 


PREVIOUS RECORD IS TACKED TO BEGINNING OF NEXT RECORD 


TS =A 
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C MAINLINE IR1ST 


IF(IBD=20) 68436859686 
TAQ=LL+1 

IBE=17 

GO TO 687 

TAQ=LL 

IBE=IBD=5 

GO TO 687 

TAQ=LL+1l 

IBE=I1BD=4 

IBC=IBC+4 

DO 680 IAA=IAQsI8C 
TIW=TAA@=LL 
NDATACIW)=NDATA(IAA) 
MDATA(IW)=LDATA(IAA) 
MWAV( IW) =LWAV(IAA) 
18D=21 

CONTINUE 

CALL EXIT 

END 


eee (CONT'D) 
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C SUBROUTINE PCNTR 


SUBROUTINE PCNTR 
SUBROUTINE PCNTR(KIDAT sMWAVsICLoICK) 


COMPUTATION OF PRECISE SPECTROSCOPIC PEAK POSITION 


AND INTENSITY USING A CUBIC FIT TO THE CENTRAL 15 PTS 
OF THE "ARRAY! KIDAT HE WAVE NUMBER POSITION OF THE 
PEAK MUST BE WITHIN THE GATE AND INTENSITY IS THE 
LOWER ROOT 


EQUATION IS Y=A30+A32*SQ+A33*CUB 
COEFICIENTS ARE DETERMINED BY ORTHOGONAL POLYNOMIALS 


411 


412 


530 


520 


INPUT 

KIDAT@=THE ARRAY OF 15 PTS IN WHICH SLOPE CHANGED 
FROM NEGATIVE TO POSITIVE 

MWAV==wAVE NUMBER POSITION OF THESE PTS 

ICL-= BEGININGS OF ARRAY 

ICK== END OF ARRAY 


DIMENSION KIDAT(25) sMWAV(25) 9DNS(25) 9P (7) 9PM(7) 9 


LDWAY(25) 


COMMON IIL sJUWAV» CONST 

DO 411 NI=1915 

ND=NI 

DNS(NI)=FLOAT(KIDAT(ND) ) 

DWAV(NI)=FLOAT (MWAVI(ND) ) 

DO 412 NJ=1l97 

NIP=16=Nu 

P(NJ)=DNS(NU)+DNS(NIP) 

PM(NJ)=DNS(NJ)I=DNS(NIP) 

AZ0=—78 e#P(1)-136#P(2)4+420#P(3)4+876#P(4)41220%P (5) 


1+147e#P(6) +1620%*P(7)+167e*DNS(8) 


A30=A30/1105.6 
A312129220#PM(1)=41216%PM(2)=14150.%PM(3)-18334e%PM(4) 


%=178426% 
1 PM(5)-138436*PM(6)=-75060#PM(7) 


A31=A31/3341526 
AZ2=91e*P(1)4520*P(2)4+190#P(3)=8e%P(4)—2964P(5) 4d, 


L¥P(6) =530*P(7)—56e#DNS(8) 


A32=A32/123766 
A3ZZ=—91e#PM(1)—130*PM(2)4+350#PM(3)+586%PM(4)4+610*PM(5) 


1+490*PM(6)+27e#PM(7) 


A33=A33/477366 

IF (ABS(A33)-00000001)53095305420 

X="A31/(20e0%*A32) 

Y=A30+X*(A21+X*A32) 

SWAV=(DWAV(8)4+X)/100e+ CONST 

Y=Y/106 

WRITE (69520) YsSwAv 

FORMAT (1HOs5Xs14HQUADRATIC PEAKs3Xs14HTRANSMITANCE =y5 


1F10¢393Xs SHPRCNTs5Xs13HWAVE NUMBER =9F100393X94HCM=1) 


GO TO 499 
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C SUBROUTINE PCNTR eee(CONT'D) 


420 DENOM = 3e *A33 
DISC=A32*A32=A31*DENOM 
IF (DISC)4219422 9423 

421 RAD=SQRT(=DISC) 
REAL=-A32/DENOM 
ROOT=RAD/DENOM 
DWAV(8)=DWAV(8)/10e+CONST 
WRITE (69502) REAL »ROOT sDWAV(8) 

302 FORMAT (1HOs5X*19HCOMPLEX SOLUTION ’ 
L22HWAVENUMBER=REAL PART 4 
2F100393Xs16HIMAGINARY PART =9F10e3910HPOSITION =y 
3F100393X s4HCM=1 ) 

GO TO 499 

422 X==A32/DENOM 
YEA3B04+X*(AZ1+X#*(A32+X#A33) ) 
SWAV= (DWAV(8) +X) /10e0+CONST 
Y=Y/106¢ 
WRITE (69521) YosSwAv 


5921 FORMAT (1HO*s16HONE ROOT PEAK 9l14HTRANSMITANCE = 


1F100393X+» 
LSHPRCNT s95Xs12HWAVENUMBER =9F100393X94HCM=1) 
GO TO 499 
423 RAD=SQRT(DISC) 
XM1=(-A32+RAD) /DENOM 
XM2=(=$A32=RAD) /DENOM 
ORD1L=A30+XM1*(A31+XM1*(A32+XM1#A33)) 
ORD2=A30+XM2* (AB1+XM2*(AZB2+XM2%A33) ) 
IF (ABS(XM1)=—700)42494259425 
425 IF(ABS(XM2) =720)42694279427 
427 SWAV=DWAV(8)/10+60+CONST 
Y=DNS(8)/10.0 
WRITE (69428) YsSWAV 
428 FORMAT (1HO95Xs36HSOLUTION IS OUTSIDE THE ARRAY 
ls1l2HR OF ARRAY=s14HTRANSMITANCE =9F100393X95HPRON 
Z295Xel2HWAVENUMBER =9F10e39 
23X94HCM=1) 
GO TO 499 
426 SWAV=(XM2+DWAV(8))/1060+CONST 
Y=ORD2/106 
432 WRITE (69429) YsSWAV 
429 FORMAT(1LHOs5Xs17HBEST SOLUTION = 4 
LIGHTRANSMITANCE =9F10¢393X» 
LS5HPRCNT s5Xs1l2HWAVENUMBER =9F10e6393X94HCM=1) 
GO TO 499 
424 IF(ABS(XM2)—720)43094319431 
431 SWAV=(XM1+DWAV(8))/10060+CONST 
Y=ORD1/106 
GO TO 432 
430 IF (ORD1"ORD2) 426942659431 
499 CONTINUE 
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SUBROUTINE PCNTR eee (CONT'D) 
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a SUBROUTINE PKFND 


SUBROUTINE PKFND 

SUBROUTINE PKFND(NDATAsLWAV»LL »JKsI9IT»LEND) 
DIMENSION NDATA(210) sLWAV(210) sKIDAT(25) sMWAV(25) 
COMMON IIL sJUWAV s CONST 


DEFINE LIMITS OF PK FIND ROUTINE 


ILK=0 

JU=0 

IF(LEND) 33393335334 
333 INF=I+6 

GO TO 335 
334 INF=I<9 
235 Lhlert<=7 


CHECK WHETHER TRANSMITANCE IS LESS THAN 360( 


DO 260 IB=INFolIFI 
IF(UK=6) 19293 
3 IF(UK=7) 19294 
G IF(UK=8) lse2sl 
2 JU=1 
IF (LWAV(IB)=7057)59591 
5 GALLCOTSTRT 
GO TO 6 
R? CALE! TSTOP 
6 CONTINUE 
IF (NDATA(IB)=30)2069207%9205 
206 IF(NDATA(1IB=-1)—30) 26092609208 
208 TWAV=FLOATI(LWAV(I8))/100e+CONST 
WRITE (69211) TWAV 
211 FORMAT (1HOs5Xs36HTRANSMITTANCE IS LESS THAN 3 PERCENT>s 
116H AT WAVENUMBER= 9Fl12e3) 
GO TO 260 
207 IF(NDATA(1B8=1)=30)26092609208 
205 IF(NDATA(IB=1)=30)20992099210 
209 TWAV=FLOAT(LWAV(IB))/100e+CONST 
WRITE (69212) TWAV 
212 FORMAT(1HOs5XsS52HTRANSMITTANCE IS MORE THAN 3 PERCENT 
* AT WAVENUMBE 
lR= »Fl2e3) 
GO TO 260 
210 IF(NDATA( IB) =800)440 9440 9442 
442 ISUMT=0 
DO 443 J=196 
IBJ=IBtU 
IBK=IB=J 
ISUM=U*(NDATA(IBK)=NDATA(IBJ)) 
ISUMT=ISUMT+ISUM 
443 CONTINUE 
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444 
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216 
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224 


221 


220 
222 


223 


225 FORMAT(1HOs5Xs32HFLAT PEAK=SLOPE CONSTANT LONGER 
127HTHAN1*5CMS1 AVER WAVE NOe =9 
2F100393Xs14HTRANSMITANCE =916911H/10 PERCENT} 


250 
260 


C SUBROUTINE PKFND eeel(CONT'D) 


GO TO 444 

ISUMT=0 

DO 441 J=194 

IBJ=IBt+u 

IBK=IB<uU 
ISUM=J*(NDATA(IBK)=NDATA(IBYU) ) 
ISUMT=ISUMT+ISUM 

CONTINUE 

IF(UK"1) 44554459214 
IF(CIB@INF) 25092509214 
CONTINUE 

IF (ISUMT) 2169218250 

IF (1TASU) 25092209219 
JIWAV=LWAV(IB) 

GO TO 250 

ICL=IB=-7 

ICK=IB+7 

DO 221 IKL=ICL»sICK 

ILK=ILK+1 
MWAVCILK)=LWAV(IKL) 
KIDATCILK)=NDATACIKL) 

ILK=0 

CALL PCNTR(KIDAT sMWAVsICLoICk) 
GO TO 250 

IF ( JUWAV=LWAV(1B)-15)22292235223 
ICL=1B=12 

ICK=18+2 

GO TO 224 

XWAV= (FLOAT ( JUWAV) )/2000 
YWAV=(FLOAT(LWAV(IB)))/2020 
ZWAV=XWAV+YWAV+CONST 

WRITE (69225) ZWAVsNDATA(IB) 


TASU=ISUMT 

CONTINUE 
CALL TSTOP 

RETURN 

END 
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A= 34 
< SUBROUTINE DETRN 


SUBROUTINE DETRN 
SUBROUTINE DETRN(NDATA sLWAV sLL sUKoI9IT »MEND) 


THIS SUBROUTINE DETERMINES WHETHER PEAKING FINDING 
ROUT SHOULD BE 

CALLEDe THE CONDITION TRIGGERING THE CALLING OF PKFIND 
IS THAT THE 

DATA POINTS MUST BE O«1 WAVENUMBERS APART. 


DIMENSION NDATA(210) sLWAV(210) 

COMMON IIL sJJWAV sCONST 

IF(MEND) 50595059506 

LEND=0 

GO TO 507 

LEND=1 

MEND=0 

T=0 

IT=0 

ITLO=LL+16 

IF(JK"1) 68196815682 
IF(LWAV(2)"LWAV(1)-1) 68396839684 

T=1 

DO 685 IB=2sILO 
IF(LWAV(IB+1)“LWAV(IB)=1) 686596863687 
IF(LWAV(IBI-LWAV(IB=1)-1) 68556855689 
1=IB 

GO TO 685 

IF(LWAV(IB)=LWAV(IB=1)=1) 69096909685 
IT=IB=1 

IF(IT=1-20) 945994559946 

CALL PKFND(NDATASLWAVsLLoJKsIesITsLEND) 
GO TO 685 

1=0 

IT=0 

CONTINUE 

GO TO 688 

IF (LWAV(17)-LWAV(16)=1)69156915692 
l=17 

DO 693 IB=17sILO 
IF(LWAV(IB+1)—-LWAV(IB)=1) 69496949940 


IF(LWAV(IB)-LWAV(IB=1)-1) 69396939941 
I=1B 

GO TO 693 

IF(LWAV(IB)@“LWAV(IB=1l)—1) 94299429693 
IT=IB=1 


IF (IT#1=20) 94399435944 

CALL PKFND(NDATAsLWAVsLLoJUKoI eo ITsLEND) 
GO TO 693 

T=0 

IT=0 
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C SUBROUTINE DETRN- eee (CONT'D) 


CONTINUE 

IF( I) 94799473948 

IF (1T)94999499947 

IT=ILO 

MEND = 1 

CALL PKFND(NDATAsLWAV»LL sUKoIsITsLEND) 
RETURN 

END 
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A= 36 
C MAINLINE BKGRD 
MAINLINE BKGRD 


SHORTEN RECORD WHEN WAVENUMBER PASSES MULTIPLE 
OF 10000e«0 CM$1 
SETS CONDITIONS FOR FIRST AND LAST RECORDS 


DIMENSION DAT(80) sWAVE(80) sKDATA(210) sIWAV(210) 
DEFINE FILE 1(459320sUsKI) 93(2693209UsIK) 
KN=0 

IY=1l 

MEND=0 

IN=1 

IP=1 

IQ=1 

DO 695 JK=1913 

READ(1L' UK) (IWAV(LQ) sKDATA(LQ) sLQ=15160) 
IF(UK=1)10910911 

CONST=FLOAT(IWAV(1)) *1000¢ 

CONTINUE 

DO 667 JJ=19160 

LL=160 

ITER=0 


THIS INDICATES WITH END OF DATA OR A MULTIPLE 
OF 1000+e0 CM$1 


IF(KDATA( JJ) +1) 66896699667 

IF (JJ@1)74557459670 

LL=uJ=1 

ITER=1 

GO TO 798 
CONST=FLOAT (IWAV( JJ) )*1000¢e0 
WRITE(69671) IWAV(UJ) 

FORMAT (1H1s21HWAVENUMBER HAS PASSEDs15910HO0060 CM=1) 
DO 672 Lyu=JU9159 
KDATA( LJ) =KDATA(LU+1) 

IWAV( LU) = IWAV(LU+1)+41L0O*ITFIX( CONST) 
CONTINUE 

LL=LL=1 

GO TO 798 

IWAV( Jd) = IWAV( JJ) +1l0*ITFIX( CONST) 


THIS IS AN INSERTION 
STORE UNSMOOTHED DATA IN FILE FOR PLOTTING PROGRAM 


DO 799 LM=leLL 
DAT( IN) =FLOAT(KDATA(LM) ) 
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* MAINLINE BKGRD eee (CONT'D) 


WAVE( IN) =FLOAT(IWAV(LM) ) 

IF ( IN=80)78897919791 

WRITE(3'IP) (DAT(IU) sWAVE(IU) » IU=1980) 
IP=IP+1 

IN=0 

IFCITER) 775977558789 

IF (LM=LL) 77596019601 

KN=IN+1 

DAT(KN) ==200 

WAVE (KN) =FLOAT(IP*804+IN) 

GO.TO 796 

IN=IN+1 

CONTINUE 

CONTINUE 
WRITE(3'IP) (DATCUN) sWAVE (UN) »UN=1 9KN) 
WRITE(69655) IP 

FORMAT(IXs'END OF DATAcceIP=' 514) 

CALL EXIT 

END 
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. MAINLINE IRPLT 
MAINLINE IRPLT 
PLOTTING IR SPECTRUM FROM FILES NOw 2 AND 3 


DEFINE FILE 2(3093209U9M) 9 3(3093209U9M) 

DIMENSION WMNMX(2) 5s TMNMX(2) 

DIMENSION DAT(80)»s WAV(83) 

DIMENSION BX(9)s BY(9) 

DATA BX / 4#0e09102900e2093*000 / 

DATA BY / 2%0009105708 s0eO 900 FD s0elOs0e0 9m1le5708 9000 / 


DETERMINE MINIMUM AND MAXIMUM VALUE FOR WAVENUMBER AND 
TRANSMITTANCE 


WMNMX(1) = 9999960 
WMNMX(2) = =9999960 
TMNMX(1) = 99990 

TMNMX(2) = =-999940 


DO 20 1=293 

DO 20 J=1930 

READ (I'J) (DAT(K)»s WAV(K)»s K=1980) 
DO 20 K = 1980 

IF (DAT(K)) 219595 


IF (WAV(K) = WMNMX(1)) 69797 
WMNMX(1) = WAV(K) 

IF ( WAV(K) = WMNMX(2)) 99998 
WMNMX(2) = WAV(K) 

IF ( DAT(K) = TMNMX(1)) 10911911 


TMNMX(1) = DAT(K) 
IF ( DAT(K) = TMNMX(2)) 20920912 


TMNMX(2) = DAT(K) 
CONTINUE 

DO 25 K=192 

WMNMX(K) = WMNMX(K)/1000 
TMNMX(K) = TMNMX(K)/10¢0 


PLOT GRAPHS FROM FILES 2 AND 3 


DO 200 1=293 

SELECT ORIGIN AND SCALE FOR THE GRAPH FOR BOTH AXES 
CALL COSGR (WMNMX929XMINsXMAX9159XSCFA) 

CALL COSGR (TMNMX929YMINsYMAX910sYSCFA) 


SET THE SCALE AND DRAW THE GRIDS 
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As. 29 
¢ MAINLINE IRPLT eee (CONT'D) 


CALL SCALF ( leOsleOs=Oe59—-165) 
CALL FGRID (0:0e03:00e09066330) 
CALL FGRID (1900090e090045320) 


LABEL BOTH AXES 


CALL FCHAR (=00709203590017900253105708) 
WRITE (7930) 

FORMAT ('TRANSMITTANCE ( PERCENT )!') 
CALL FCHAR (7039007290017 90025 5000) 
WRITE (7940) 

FORMAT ('WAVENUMBER ( CM=1 )#) 

CALL FPLOT (090003060) 


ANNOTATE BOTH AXES 


CALL SCALF (1e2/XSCFA9009/YSCFAsXMINsYMIN) 


BX(1) = XMIN 

Bx(2) = YMIN 

BX(4) = XMAX = XMIN 
BX(5) = XSCFA 

Bx(9) = XMIN 

CALL ANOTS (BXs0608900179le2/XSCFAs009/YSCFA) 
BY(1) = XMIN 

BY(2) = YMIN 

BY(4) = YMAX = YMIN 
BY(5) = YSCFA 

BY(9) = YMIN 


CALL ANOTS (BYs0e08900e17s1le2/XSCFA9069/YSCFA) 


PLOT THE GRAPH 


ICNTL = =2 

DO 100 Jv=1930 

READ (I'JU) (DAT(K) sWAV(K)»s K=1980) 
DO 100 K = 1980 

IF (DAT(K) )150s50950 

X = {(WAV(K)/10¢0) 

Y = (DAT(K)/10¢0) 

CALL FPLOT (ICNTL9XsY) 

ICNTL = 0 

CALL FPLOT (19XMAX#20e5/(le2/XSCFA) sYMIN=220/ 
1(O069/YSCFA)) 

CONTINUE 

CALL EXIT 

END 
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24 
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A= 40 
a MAINLINE DELTA 


MAINLINE DELTA 

DIMENSION DAT2(80) sWAV2( 80) sDAT3 (80) sWAV3(80) 
DIMENSION DEL23(80) 

DIMENSION WMNMX(2) 9s OMNMX(2) 

DIMENSION BX(9)s BY(9) 

DEFINE FILE 2(4593200UsIK) 93(3093209UsIL) 
194(3093209UsI0) 

EQUIVALENCE (DEL23(1) »DAT3(1)) 

DATA BX / 4*¥OeOs1629002093%*000 / 

DATA BY / 280009105708 0 e0 Oe sVelOsDeOem1le5708 9000 / 


CALCULATE THE DIFFERENCE B/W FILES 2 AND 3 AND STORE 
IN FILE 4 

CALCULATE MINIMUM AND MAXIMUM FOR BOTH VARIABLES 

OF FILE 4 SIMULTANEOUSLY 


WMNMX(1)=9999960 

WMNMX (2) ==9999940 

DMNMX (1)=999960 

DMNMX (2) ="999960 

DO 20 I=1970 
READ(2'1) (DAT3 (4) sWAV3 (4) 9U=1980) 
READ(3'1) (DAT2(U) sWAV2(4) 9J=1980) 
DO 15 y=1s80 

IF (DAT2(4)) 219595 

IF (WAV2(4) -WAV3(U) 19996999 
DEL23(U)=DAT2Z2(U)=DAT3 (4) 

IF (WAV2 (4) —-WMNMX(1))79898 
WMNMX (1) =WAV2( J) 

IF (WAV2(0)=WMNMX(2))1091099 
WMNMX (2) =WAV2 (J) 

IF (DEL23(0)=DMNMX(1))11912912 
DMNMX(1)=DEL23(u) 

IF (DEL23(U)=DMNMX(2))15913913 
DMNMX(2)=DEL23(J) 

CONTINUE 
WRITE(4'1) (DEL23 (4) sWAV2(U) »y=1980) 
CONTINUE 

GO TO 24 

DEL23(u) = 99999640 
WRITE(4'I1) (DEL23(K) sWAV2Z(K) 9K=1 9¥U) 
DO 25 I=le2 

WMNMX (1) =WMNMX(1) 71060 

OMNMX (1) =OMNMX(1)/1000 

CALL COSGR (WMNMX929XMIN*sXMAX915sXSCFA) 
CALL COSGR (DMNMX929YMINsYMAX910sYSCFA) 
CALL SCALF ( leOsleOs=Oe59—-105) 
CALL FGRID (09000904e09006930) 

CALL FGRID (19000900090045:20) 

CALL FCHAR (=02709105090017900259165708) 
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c MAINLINE DELTA eee (CONT!D) 


WRITE (7930) 
30 FORMAT('TRANSMITTANCE DIFFERENCE ( PERCENT J!) 
CALL FCHAR (7039=0072900179002590s0) 
WRITE (7940) 
40 FORMAT ('WAVENUMBER ( CM=1 )*) 
CALL FPLOT (090609000) 
CALL SCALF (102/XSCFA9009/YSCFA sXMIN»YMIN) 


BX(1) = XMIN 

BX(2) = YMIN 

BX(4) = XMAX = XMIN 

BX(5) = XSCFA 

BX(9) = XMIN 

CALL ANOTS (BXs0eO08s00el7sle2/XSCFAs0e9/YSCFA) 
BY(1) = XMIN 

BY(2) = YMIN 

BY(4) = YMAX = YMIN 

BY(5) = YSCFA 

BY(9) = YMIN 

CALL ANOTS (BY 90 e08s0el7sle2/XSCFAs00e9/YSCFA) 
ICNTL = =2 


DO 100 J=1s70 

READ (4'J5) (DEL23(K) sWAV2(K) sK=1980) 

DO 100 K = 1,580 

IF (DEL23(K) = 9999.0) 5091509150 
50 X = (WAV2(K)/1000) 

Y=(DEL23(K)/10.0) 

CALL FPLOT (ICNTL9XsyY) 
100 ICNTL = 0 
150 CALL FPLOT (19XMAX+205/(102/XSCFA) sYMIN=260/ 

1(0e9/YSCFA) ) 

CALL EXIT 
99 WRITE (69999) WAV2(U)»s wWAv3(u) 


Am SL 


999 FORMAT (10X»'WAVENUMBER SYNCe ERROR = '32(2X9F90l)) 


GO TO 21 
END 
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Ar. “2 
G MAINLINE KINET 
MAINLINE KINET 


THIS PROGRAM COMPUT PERCENT CONVERSION OF ALCOHOL 
ALSO THE PARTIAL PRESSURE OF EACH COMPONENT IN THE 
PRODUCT MIXTURE AND THE REACTION RATEs 


INPUT DATA REQUIRED==FEED RATE AND CATALYST WEIGHT 


We==WEIGHT OF THE CATALYST» GM 
FAO@===FEED RATE sG=MOLE/HRe 


DIMENSION IWAVE(160) sIDATA(160) 
DEFINE FILE 1(4553209Usk0O) 
I1T=0 
J=0 
RTAL=O06 
400 J=Utl 
READ(1'J) (IWAVE(K) sIDATA(K) 9K=19160) 
IF(J-1)39391 


3 N=2 
GO TO 4 
1 N=1 


4 DO 2 IR=N9160 
IFC IDATACIR)=1)75735 
5 RTAL=RTAL+FLOAT(IDATA(IR)) 
Il=I1+1 
GO TO 2 
7 WRITE(69102) IWAVE(IR) sIDATACIR) sII 
L102 FORMAT(1Xs!WAVE NOe=!sT4e'TRANSe=' sI4e'AT POINT' s I6 
*»' TNCORRECT') 
IF (II=2000)4013:40196 
401 I[1=I1I+1 
2 CONTINUE 
GO TO 400 
6 RI=RTAL/II 
RZ2=(RI/100e0)+6e0*(RI=“300e)/706 
R =(RI/1000)+6e0*(R2—300)/706 
ABS=(AL0G(0e81) =ALOG(R/10060) )/263025 
WRITE(69100)RoII 
100 FORMAT(1Xs'AVGse PERCNT TRANSe=! 9Fl00293Xs'!TOTAL 
* POINT='o9159/) 
WRITE(69101)ABS 
101 FORMAT (1X»s'ABSORBANCE='! 9F15e6) 
READ( 5910) FAQ9wW 
10 FORMAT(2F 1566) 
PO=02609518+92e694480*FAO 


COMPUT PA FROM CALIBRATION FUNCTION 


PCAL=(-02e0027+ABS) /0015966 
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c MAINLINE KINET eee (CONT!D) 


CORRECT FOR TEMP=170 Ce 


PARS=PCAL#(273e+1700)/(2730+29¢8) 

CONV= (Lem(PAR/PO) )/(Let(PAR#0 eO951E*FAC) /PO**2 ) 
PCONV=1006*CONV 

WRITE (69103) PCONV 

FORMAT(1Xs'PERCNT CONVERSION=! 9F 1064) 
RATE=CONV*FAO/W 

WRITE (69104)RATE 

FORMAT(/91Xs'*REACTION RATE='sF1506) 
PA=PAR*®(PO+(OeO9518*FAO*#CONV) ) /POX*2 
PW=PO=PA 

WRITE(69105)PAsPW 

FORMAT(/91X9 'PA='sF100395Xs'!PWsPP=' 9F1003) 
CALL EXIT 

END 
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B. COMPUTER PROGRAMS FOR MODEL DISCRIMINATION 


The appendix contains the programs used in this thesis. 


The notation used in the following programs is described here. 


AO 


IT 


LAM 


Dispersion matrix xx 

Prior estimate of the dispersion matrix 

Least squares estimate of the parameters (vector) 
Variance-covariance matrix 

Determinant, normalizing matrix 

A function subprogram to set up the design matrix 
Errors vector 

A fun-tion subprogram to compute the model predictions 
Mesh size in the design of experiments 

Reference index for the competing models 

Expected likelihoods 

A parameter in Marquardt's algorithm 

Number of parameters in the given model 

Prior estimate of the parameters (vector) 


Number of data points 


-Total number of competing models 


Independent variables 
Negative of the expected entropy change 


Minimum sum of errors squared 


ren bedit32eb 


if beeu 
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ef 
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SETUP 


B-2 
A function subprogram to specify the prior parameter 
distribution for the given model 
Error variance 
Vector of observations 


Expected value of the prediction from the entire methamodel 
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MGT 
0798712325 On ONG 3 She? S06 O, 4358432714 PARAMEVER VALUES AT 
0.08485636293 0.01001486825 0.1888224163 VARIOUS ETERALTOUS 


0.09233897556 0.008185821263 0.2382591273 OF HMARQUARDT'S ALG- 


B-3 
SAMPLE OUTPUT USING THE PROGRAMS OF APPENDIX B. 
DEHYDRATION MODEL M3,SPECIFIED HERE BY LM=2. 
Lile«2 
SPECIFY THE PARAMETER PRIOR VARIANCE-COVARIANCE 
MATRIX AND THE STARTING VALUE FOR THE PARAMETERS. 
AO0+3 390 
Mize. €; DN3 ek 1 
MARQURDT'S(1963) ALGORITHM IS USED TO COMPUTE THE 
LEAST SQUARE ESTIMATE OF PARAMETERS AND VARIAINCE- 
COVARIANCE MATRIX IS COMPUTED AS DESCRIBED IW 
CHAPTER TWO. PROGRAM ‘IIST' IS USED TO OBTAIN 


LHE MODEL LIKELIHOOD.CALL '4GT' TO PERFORM THE 


ABOVE CALCULATIONS. A SAMPLE OUTPUT IS; 


0.09519688775 0.007873153225 0,.2560221165 ORITHH. 
Pe O2Sh 3h 7 2478" 0 007855672 14 © On 25747044119 


PARAMETERS 


0.09540525955 -907007855869105 0.2574720037 


95 PCWT LMT: 0.007987637281 0.001980287403  0.04827574289 


VARIANCE 


VAR-COVAR: 


4.125640078F 5 STD DEV: 0.00642311457 


1. 3950587736 5 ~2,.857511393F 6 9.-180G706332 5 


2eaeWsi1s9In 6 -9.80agu5500R 7 °2.1596601978 75 
9.1306706337 5 2.1526691977 5 5.826368378F 4 


MODEL LIKLIUD: 3.436255498217 
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PROGRAiIL =‘ INy! 


LUIS PROGRAG COMPUTES THE IuVERSE OF A NOH- 


OINGULAR UATRIX BY GAUSL-JORDOI PIVOT, 


VINVLIIIV 

V 2eTAY* ths Ted sK3P 
IR Epo] A= / pil 
>O,pli<'VO INVERSI POUND" 
1#Q(1 0 +o!) p(,Q/1) Jel <Pei1Te1tpi! 
MUR, 23 eC KEK PRA EM T2145 | 
P<10d/? ,OpPLK,1)]+P[1,x] 
ML 31t+pod Je~d 
-2x 1h 90> /)0231 340 / 
't#1910[1)] M-( 7x0 31] 
+x 1O047¢I-1 


feat 
joes Pe let(h 1s Jl 134) 
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PROGRAGE Lp ay a 


TOES) OBTAINS, TUR -VRTERUIUART.-OF 4 


MQUARE HATRIX BY CHIO'S NETHUOD., 
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PROGRAM ‘ysy' 


LHIS PROGRAM COMPUTES THE EXPECTED LIKELI- 
HOODS FROM THE SUFFICIENT STATISTICS OF THE 


DATA WHEN IMPROPER UNIFORM PRIORS ARE USED. 


VASTCUOIV 
V HST3M3N 
4+ oB 
N<«oX0 31) 
L«( (2x9 VAR) «(-M42))*x(0.5%"(i52))x(e(-SSt(2xVAR))) 


PROGRAH = "SS" 


WET OCOURUTES THE EXPRCTEDSLTIKELINOOD OF a GIVEN 
MOBHIDARROM MDI QUE BICTION TONGA REDTICS-OF DHE DKA 


We ee LOPED Til CHaArien 2 Of YHIS Tileoia. £1 


Lae 
Uv 


POREERIOR DISTRIBUTION AND £2 USES HE PRIOR 


Prete pUrloy OF Tie PA VIEVERS . 


JNISSCTIIV 
V HOSA ZAA SAB SAC 3511340 30D 54 
W+o£l 31 J 
DVD*«DETCA0Ot+A ) 
L11<€C2*xOVAR )*(-142))xC*(-SS4(2xVAR))) 
AA+(A0+.x/1U0 )+A+.xB 
AR<«(A0+. X/1U/0 )+2xAt,. x3 
AD<( (QB )+.*%A+.*B)+(QHAA 4.x CLTUV Auta )+. xia 
L1e(*-0.5*x(AD-(Qi3 Dt. x (INV AU4+2xA )+.X%a55))xL11 
hi<(L1i (OEY (AQ+AxK2))*0.5)*xDD*0.5 
AR<((Q3)+4.*A+t.*2 )-(Q4A )4+.* (TV AUtA)+.%*dA 
AC#(AL+( (QitU0 )+.*A0+.*.1U0) ) 
L2*LiixC( (DET Ad)*0.5)x(* O.5xAC)E(DV*0.5) 
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C. CALIBRATION OF LIQUID SYRINGE FEEDER 


The Sage syringe pump was used to inject liquid feed into 
the reaction system. The syringe feeder was calibrated by feeding 
freshly distilled water at room temperature of 24°C. The water was 
collected in a 25 ml pycnometer after it has passed through a syringe 
needle which is just big enough to fit in the capillary plug. To 
prevent the loss of water due to evaporation, the pycnometer was cooled 
in an ice-water bath. The duration of each run was one hour. The 
amount of water was obtained by weighing the water at 24°C. Each 
piston speed was calibrated four times in such a way that the whole 
range of the piston travel (about 12 cm) was covered. The results 
showed that the speed was reproducible within + 0.1% with respect to 
the position of the syringe piston. The average values of the four 
runs for each speed is used to obtain calibration curves. For the 
speed range of 1/1000 X, six data points (average values) were collected 
and of 1/100 X five were obtained. These two sets of data were fitted 
by the least squares method and the results are shown in Tables C-1 
and C-2 for each speed range. The computer programs, LEAST, of the 


least squares method are also listed below. 
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TABLE Cle CALIBRATION OF THE SYRINGE FEEDER BY 


LEAST SQUARES METHOD 


X= % OF (1/100 


Y= FEED RATE= CC PER HOUR AT 24 DEGe 


THE COEFFICIENTS OF THE POLYNOMIAL» 


AQ = 0 200690 
Al = 0205274 
A2 = -0.00000 
X MEASURED Y OBSERVED 
100000 005211 
2020000 120433 
3020000 15666 
4020000 200935 
600000 301300 
8020000 41640 
VAR TANCE = 
STANDARD DEVIATION = 
MAXIMUM PCT ERROR 


OX) ON THE CONTROL DIAL OF THE FEEDER 


REGENERATED DATA 


= 00255461 


06000003 


06001803 


Y CALCULATED 


005197 
160449 
165685 
200907 
301304 


401641 


(LOW SPEED RANGE) 


YEAOtTAL*X+AZ#X#HH2 5 


PCT ERROR 


022554 
Oe1l564 
0«1264 
001329 
020120 


020029 


C~2 
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TABLE C™2e CALIBRATION OF THE SYRINGE FEEDER BY 


LEAST SQUARES METHOD (HIGH SPEED RANGE) 


X= % OF (1/100X) ON THE CONTROL DIAL OF THE FEEDER 


Y= FEED RATE= CC PER HOUR AT 24 DEGe 


THE COEFFICIENTS OF THE POLYNOMIAL» Y = AO + Al * Xo 


AQ = =0 203847 
Al = 0652202 
REGENERATED DATA 
X MEASURED Y OBSERVED Y CALCULATED PCT ERROR 
102000 Se211 50181 02561 
202000 10439 106402 0«360 
302000 152575 156622 00304 
40000 204765 204842 O«372 
606000 370341 312283 02184 
VARIANCE = 006003467 
STANDARD DEVIATION = 02058884 
MAXIMUM PCT ERROR = 04561762 
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C MAINLINE LEAST 


MAINLINE LEAST 


ee ee ee ee ee ee ee ee ee ee ee ee ee ee 


ANY ORDER UP 


X(T) on 
Y(T) - 


ke OK OK Ok OK ok OK Ok KK ke ok OK KK OK OK Kk KK OK OK OK OK 


MAINLINE LEAST 


THIS PROGRAM WAS WRITTEN FOR FITTING A MAXIMUM OF 
50 DATA POINTS TO POWER SERIES TYPE POLYNOMIALS OF 
TO A MAXIMUM OF FOURTH DEGREEes 


INPUT DATA 

NCASE = NUMBER OF SETS OF DATA 

NCOPY = NUMBER OF COPIES OF OUTPUT DESIRED 

N - NUMBER OF DATA POINTS 

M - DEGREE OF POLYNOMIAL 

NTL = NUMBER OF CARDS FOR TITLE 

NPAGE = PAGE NUMBER OF OUTPUT 

NPLT = DATA REGENERATION FLAG 
ee eO@=REGENERATE GIVEN DATA ONLY 
eee lREGENERATE GIVEN DATA PLUS 20 

INTERMEDIATE POINTS 

DES(K) = ALPHANUMERIC DESCRIPTION OF THE TITLE 

XNAME = ALPHANUMERIC DESCRIPTION OFX 

YNAME = ALPHANUMERIC DESCRIPTION OF Y 


INDEPENDENT VARIABLE 
DEPENDENT VARIABLE 


HER HH HIE HH HT HHH HHH HHH HHH HHH EH KEK EH EK EE IE SE 4G HE SE HE SE SE ee EH 


KK RK KK OK KK OK KK OK KK OK KK KK KK KK OK OK 


DIMENSION X(50) 9¥(50) 8A(5095) »P (20920) 9V(20) »Z(20)»5 
1DES(10%15) »sSNAM(5) »XNAME(100) »YNAME(100) 


DATA SNAM/'AO ='s'Al ='9!A2 


READ(59s1) NCASE sNCOPY 
DO 9 NC#=1lsNCASE 
READ( 591) 


FORMAT (515) 

DO 11 NT#=1leNTL 
READ( 5912) 
FORMAT (15A4) 
READ (5920) 
FORMAT (15A4) 
READ (5921) 
FORMAT (15A4) 
FORMAT (10X915A4/) 
MM=M+1 

DO 2 I=leN 
READ(593) X(1T)sY(1) 
FORMAT (2F1005) 

DO 4 I=ls9N 

DO 4 J=l»sMM 
A(IsJ)=X(1)**(U=1) 


='sgta3s alytag =t/ 


NoMsNTL eNPAGE sNPLT 


(DES(NT 9K) sK=1915) 
(XNAME(1) 9 I=1915) 


(YNAME (J) 9v=1915) 
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C MAINLINE LEAST eee (CONT'D) 


DO 5 T=l»sMM 

DO 5 J=lsMM 

P(IsU)=O06 

DO 5 K=loeN 
P(ITsJJHFP(ITsU)+tA( Ko I) *#A( Kou) 
DO 6 I=l»MM 

V(IT)=O6 

DO 6 J=leN 

VIT HVT) FY 00) #AlCUeT) 

CALL GAUSS(P9VsMMoZ) 

DO 16 ICOP=1lsNCOPY 
WRITE(6910) NPAGE 
FORMAT('1's///966Xs'A='sI29/) 
DO 17 IT=leNTL 

WRITE (6913) (DES( 19K) sK=1915) 
WRITE (6930) (XNAME( I) sI=1915) 
FORMAT (///7912X915A4) 

WRITE (6931) (YNAME(JU) »¥=1915) 
FORMAT (12X915A4) 

WRITE(698) 

FORMAT( ///910X9'THE COEFFICIENTS OF THE POLYNOMIAL 
ls'AREA'/) 

DO 15 IT=1l9MM 

WRITE(697) SNAM( I) 9Z(1) 
FORMAT (15XsA49F11¢e5/) 

CALL REGEN(XsY9ZoMMoN) 
IF(NPLT) 999914 

CALL POLYT (X9Z9N9MM) 
CONTINUE 

CALL EXIT 

END 
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¢ SUBROUTINE POLY 

SUBROUTINE POLYT 
HHH KH HEH HHH HK HE IE HEH HE HK HE IE HE IE SE IE IE IE SEE FE SE HE SE FE IE SE IE JE 9 96 E96 36 26 56 96 26 
SUBROUTINE POLYT 


* 
* 

* POLYT SUPPLIES REGENERATED DATA AT POINTS INTER= 
* MEDIATE TO THE GIVEN DATAs 
* 
% 


HK OK kK Ok Ok 


HEHE TE HE HE KE HE HE HIE HE HE HE HE HE HE FE HE HEHE IE SE IE HEE FE IE HE FEE HE IE SE IE IE IE SE IE XE SE IE ESE SE BEE 5 5B 5K 56 EK 9K 


SUBROUTINE POLYT(X9Z9NoMM) 

DIMENSION X(50)92(20) 

WRITE (691) 

FORMAT(///932Xs'PLOT TEST DATA'//25X!X CALCULATED! 54x 


ls'Y CALCULATED!/) 


XMAX=06 

XMIN=999996 

DO 2 I=loeN 

IF (XMAX=X(1)) 39394 
XMAX=X(1) 

IF (XC I)—XMIN) 595592 
XMIN=X(1) 

CONTINUE 
DELX=(XMAX=XMIN)/206 
XY=XMIN 

DO 6 I=1s20 

CALHA06 

DO15 J=1lsMM 
CAL=CAL4+Z (4) #XY**(U=1) 
WRITE(697) XYsCAL 
FORMAT (24X%92(F100395x)) 
XY=XY+DELX 

RETURN 

END 
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C SUBROUTINE REGE 
SUBROUTINE REGEN 


HEHEHE HEHE HE HEHEHE IE HEHEHE HEHE HE HEH HEHEHE HEHE EK HEHE IE IEE HE HE IE HE EEE EE IE SE HEHE HE 96 98 96 96 HK 


SUBROUTINE REGEN 


CALCULATES THE VARIANCE AND STANDARD DEVIATION OF 
THE FlTs 


* OK K kK OK K 


* 

* 

* 

* THIS SUBROUTINE REGENERATES THE GIVEN DATA AND 
* 

* 

* 

* 
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SUBROUTINE REGEN(X9Y9Z9MMoN) 
DIMENSION X(50)9Y(50)9Z(20) 
WRITE(691) 
1 FORMAT (///929X9 "REGENERATED DATA'//10X9'X MEASURED! 55x 
ls'Y OBSERVED' s5xs'y CALCULATED! 93Xs'PCT ERROR',g/) 
VAR=0+6 
HI =06 
DO 2 T=leN 
CAL=06 
DO 32 J=l»MM 
3 CALSCAL+Z(4U) *X( 1) **(U=1) 
CAT=ABS(Y(1I)=CAL) 
PCE=CAT/Y(I)#¥100. 
VAR=VAR+CAT#¥2 
IF{HI=PCE)49452 
HI=PCE 
WRITE(695) X( 1) 9Y(1I) sCALsPCE 
FORMAT( 9X94(F106395X)/) 
VAR=VAR/(N=1) 
DEV=VAR*#005 
WRITE(696) VARsDEVsHI 


un + 


& FORMAT(//910X'VARIANCE ="9F10066//10X> 
1'STANDARD DEVIATION ='59F1006//10X¢9 

Z2'MAXIMUM PCT ERROR =!5F1006) 
RETURN 


END 
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e SUBROUTINE GAUS 
SUBROUTINE GAUSS 
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SUBROUTINE GAUSS 


* * 
* + 
%* THE FUNCTION OF THIS SUBROUTINE IS TO SOLVE THE * 
* SET OF EQUATIONS A*X=B USING GAUSSIAN ELIMINATION * 
* AND BACK SUBSTITUTION ROTATING ASOUT THE ELEMENT # 
* OF MAXIMUM MODULUS. * 
* * 
% x 
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SUBROUTINE GAUSS (AsRoNox) 
DIMENSION A(20920) sR(20) sX(20) 
M=N=-1 

DO 11 JY=leM 

S=0.6 

DO 12 I=JU9N 

U= ABS(A(I9uU)) 

IF(U=S) 1291298112 

$=U 

L=] 

CONTINUE 

IF(L=J) 1199199119 

DO 14 I=JUeN 

S=A(LosI) 

Al(LseI)=A(UsI) 

AlJsI)=S 

S=R(L) 

R(LI=HER(4U) 

R(J)=S 

IF( ABS(A(J9J))—10eE=30) 1159115915 
WRITE (693) 

GO TO 500 

MM=J+1 

DO 11 I=MMsN 

IF( ABS(A(I9J))—-1eE=30) llslllsill 
S=A(Js4)/A(I9u) 

A(IsJ)=040 

DO 16 K=MMeN 

A(T sK)=AlCU9K)=S#A(I 5K) 
R(II=R(4J)—S¥R(T) 

CONTINUE 

DO 17 K#1l9N 

I=N+1=“K 

$=020 

IFC I-N) 1179179117 

MM=I4+1 

DO 18 J=MMsN 
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SUBROUTINE GAUSS eee(CONT'D) 


S=StA( 194) *X(U) 
X(TI=(R(1)=S)/A(I91) 


RETURN 
FORMAT 
END 


(1H 
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D. PORE DIFFUSION LIMITATIONS 


The method of estimating the effectiveness factor, n (the 


ratio of the actual rate to the rate which would be observed in 


the absence of diffusion limitations), for kinetic expressions such 


as Hougen- 


Watson type of rate equations, as reported by Satterfield! !08) | 


was employed to determine whether concentration gradients existed 


within the alumina catalyst wafer. The following parameter values 


are required: 


and 


The terms 


follows: 


Da 


Yeat 


Dy = 9700 + y, + (T/M)!/2 (D.1) 
_Lé 
Kee ewe Dy Z(K.v./D.) 1/w (D.3) 
1 
WaT A EK, (Py + (PyvjDQ/D;)1 (D.4) 


within the expression for Das dy» K and w are defined as 


effective diffusivity of reactant based on Knudsen diffu- 
sion within the pores 


pore radius, cm 
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7, .42 
M = 


re 


[> = 


K,w = 


i= 
Vv = 


Ca 


cussed in 


using the 


absolute temperature, °K 
molecular weight 
Thiele diffusion modulus 

the observed reaction rate, g-mole/sec. cc of catalyst 
ratio of the catalyst volume to outside surface area 
through which reactant has access, cm 

parameters used to evaluate n 

subscript for species other than reactant A 
stoichiometric coefficient 


concentration of the reactant, g-mole/cc 


The rate of 2-propanol dehydration is determined as dis- 
Section 5.3.6. The effectiveness factor was calculated 


highest -1y/Cy value obtained in this work at 252.7°C 


LPs = 0.03947 g-mole/hr-g-catalyst, Pa = 0.479 cm Hg, and Py, = 


0.570 cm Hg; see Table 5-4) which would give the most severe, if any, 


diffusion 


Table 5-4 


limitation. 
fe) 
From Section 3.2, the Ye was measured to be 84A and from 


Bi is taken to be 


0.03947 x 0.138/(0.012 x ¥ x (2.54)°x 3600) Sere 


where the 


sec-cc catalyst 


thickness of the wafer is about 0.012 cm, weight of the 


catalyst is 0.138 gm and diameter of the wafer is 2.54 cm. Using 


ideal gas 


D-2 


law, the external concentration of 2-propanol] at 252.7°C is | 
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D-3 
Ca = 0.479/(76 x 82.06 x 523.7) g-mole/cc. The ratio of catalyst 
volume to outside surface area, L, was about 0.006 cm. Thus, from 


equations (D.1) and (D.2), we get 
Dy = 7.1596 x 107° and 4, = 0.21 (D.5) 


The best form of rate equation for the dehydration was found to be 


kPa 


(Tq) = (ERP, 


(D.6) 
A wt Kaky 


Pay) 


The diffusivity is approximately inversely proportional to the square 
root of molecule weight by equation (D.1). Thus D,/D,, = /18/60. 

From Equations (D.3) and (D.4), by assuming Buen rption constants of 
2-propanol, Ky, to be 0.007856 (b, in M3) and of water, K,, to be 
0.25747 (b, in M3), we have 


w = 1.03 and K = - 3.6 


108) | the above 


Hence KP) = -0.0227. Using Figure 4-2 of Satterfield! 
parameters dy and KPa give an effectiveness factor for this reaction 
system of virtually unity. Since the “1y/Cy ratios employed in 
other experimental runs were always lower than the value used in 
the calculation, pore diffusion can thus be neglected. 

In case of dehydrogenation, the effectiveness factor calcu- 


lated using the same criterion as that of dehydration was also found 


to be unity. Therefore, only chemical reactions are considered to be the 


rate-controlling steps. 


#7 iT it 


e-9 | 
$eyistso to ofssy NT .29\stomp (S/ES2 x 30.89. agyversio “pp 
movt ,2uiT .mo 300.0 bods! dow ,J ,sore sostwe ‘sb teduo od smuTov rf 
t9p aw ,(S.0) bas (F.0d) enotisups 
(2.0) [8-0-5 bos “OPin O8at.t = 40 


2 = 
: @ 


od 02 bivot 2aw neidavbyieb. afd yor notteups ater to mot teed onT 
av : 


(9.0) eae ; + = (,1-). 
Cia ytat tye ee OR 


gisup2 sit oF fenotsiroqorsq y[oersvit \issemtxovggs at V3ivieuvtttbh At 


Gaver’, yar 


to 2tnedenos nottaroebs oes Fran yd .(b. a) bie (€.G) snot dupa no 


\y0 aunT . (T. q) no tt eups vd tip tow sfuosfom to 100 


sd oj cyt «Stew Yo bas (EM nt of) aastod. 0 sd of Wi « Tonsaor-§ 
oved ow c(EM nt od). rO88.0 


, "9 96 a 
8.6 - = } bas EO.T @w = : 
ba 


(20f) | | | (oe 
svods ort . biaftisdtee to S-) sy uptt ante .VSS0.0- = av sonsH — 


notionay ated 40? 10F967 ezonev ti29tte ae Bid. es bas ail 
| q | . wenger! na a" arts eohta atau uttaurtty “Yo msde (2 a 
+ | ee ely an 
“< + Beau eu s¥ ‘ei ) bien inal sins Hy A ca o 
: a r 104 \ notsal ~~ > f . oft 
weptaed a =e eis $i 


| i 5 { 
L 









E-] 


E. SAMPLE DERIVATIONS OF RATE EQUATIONS 
E.1 Dehydration, (dual dissimilar-site mechanism) 


A+ 2, = At, (Eel) 
AL, + &o + PR, + Wen CE) 
PRy =P P +R (E.3) 
We, = WH hy (E.4) 


where notations are identical to that given in Section 5.3.6. Assume 


surface reaction (E.2) to be rate-controlling, then - ry = V9 
2% as, C25 where the subscripts refer to the reaction step numbers 
in the above reaction scheme. The other reaction steps are assumed 
to be in equilibrium thus Mya" els Sela 0 

or expressing in equilibrium constants 


K, - ae /Pate. 4 K3 = peta pg. and Ka = Pe Cue, 


Now, total number of active sites, hy and hos are assumed to be constant 


at Ly and Los respectively, then 
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+ 
Cy 
! 


= ce, (1 + KPa + Pp/K) 


L, = C2, +C : = Ch, (1 + Py/Kq) 


The rate expression can now be written as 


KoKy LiLo Pa 


Ty = RkyPace Ch, = 


Pp PY 


By assuming K3 + , as discussed in Section 5.3.6, the above equation 


reduces to 
Saad ve KK LL oP A/(1 + Py/Ky) (1 4 KyPa) 


which is the form of M3 in Table 5-3. 


Ewe Dehydrogenation, dual dissimilar-sites mechanism (Scheme B) 


A+ g, — Ag, (ey) 
AL, tho = XR. + Heo (E.6) 
X24 7 Xo ao Key a Ho (E.7) 


2Heo = Ho + 2hn (E.8) 
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E-3 
Kose ==) Ket hy (E.9) 


again the notations are identical to that given in Section 5.3.8. 
(i) Assume both adsorption of A, step (E.5), and surface reaction 
step (E.7) to be rate-controlling, then 


ey ili coal") |) a eomee gee keP aCe, - K sone, = 7X, CL, (E.10) 


For steps other than E.5 and E.7, they are assumed to be in equilibrium, 


i.e. Cee Gr gg 0 or by definition of the equilibrium constant 





2 
Cxe Cue Eh eee P, CL 
Kote + Kg = and k, = 4 
AX 2 C Kg 
1 H25 1 
On rearrangement, we get 
Ka, = Py Co 4/Kg (ETA 
“He, ~ vPa ~ CLoN Ke (Ee le) 
Cyo ~ Key Ke Cae IVP. (E213) 
| ] 2 
From equation (E.10), 
ke ko 


Cae, “ee. VAIS tae ee “Xn, C2, (E.14) 
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Combining equations (E.11) to (E.14), we have 


on wade: 
AL, = ke PCs /k_. (1 12 isp Cho TWP, ond 
3 2 
and 
k, KevKe 
2 


Again, total number of active sites for each type, hy and Los is as- 


sumed to be constant, thus 


Substitute expressions of fie ne Eke, and “He, into the above 


two equations, then 


Ph 
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k. PAC] eee 
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Substitute L, and L, into the rate expression (E.10), one gets 


k7KevKe ke pee 12 


: : Ace; C2. 
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“541+ RD (Iti + pp) Py (Tg) (kePayPy, +KeVKak2P 4) 
9 8 Ho 2 8 2 
Since acetone has no effect on the reaction rate and since it was not 
observed on catalyst surface, it is assumed that Ky +o, By incorporating 


this assumption into the above rate expression, we get 


1: Ky Kevkg ke Pa Lyle 
A Py K-K-vK, P k5LK-VK 
2 56a A dee 06 
Combining all the constants the rate equation is in the form of equation 
(E.15)~which is N3 in Table 5-7. 


Bian aA 
a rn = Pre @| m2 boPa + b3 7P i + by An (E.15) 
2 2 


(ii) Again assume both reaction steps E.7 and E.8 to be rate con- 


trolling, then 
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= = = = = 2 
(-r,) ails peopl ko “xe, CL, ke “He, ks Po CL. 


and, “AL, = KEP aCe, and “Xa, = Ke “An Pee, 


Substitute expressions of Cry and Cye into the rate expression, we 
1 1 
have 


k 3 


2 s 2 
From the above equation, it can be seen that the expression for Cue 
2 


is not a simple function. 


E.3 Dehydrogenation, dual dissimilar-sites mechanism (Scheme C 





(E.16) 


AL, + Ly — XL, + He, CESa7) 
XR, + Lo > Ki, + He, (E.7) 
2H2. —— Hn + Zhe (E.8) 


Ka, = 1K A Ly (E.9) 
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the notations are also identical to that given in Section E.2. Assume 
both adsorption of A (step E.16) and surface reaction (step E.7) to be 


rate-controlling, then (-r,) tan 6 D5 


or (-r,) = ky Bis CL, = kKigPa Ch. - k_16 “AL, (E.18) 


In terms of equilibrium constants for steps E.17, E.8, and E.9, one 





gets 
PCR Pu * cg Ky ng Ch, 
Cy = RS Cie +1 a ange eye = 6 - 
1 \ 2 8 He, 
From equation E.18, 
ky Paleo - K-16°AL, Lae 
Cy. = oe E.20 


Solve for Cy and Cy from equations (E.19) and (E.20), we have 
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Assume Ly and Lo are constants as in the previous section, 


Ly = Ca, + ae a TH 


Ky ok, Pace PCE 
= Ch + TB A) eg (E.21) 


-16 
Rants + Ky 7k7Ca, 


ret hd 0S ek se & 
2 2 AL, HR, 


kKi¢ Pa ¥PHo Cho me Ch, 


tee 
2 VK 


It can be seen from equations E.21 and E.22 that C2, and Ch. can not be 


solved to give simple expressions. 
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F. THERMODYNAMIC ANALYSIS OF THE DEHYDRATION 
AND DEHYDROGENATION OF 2-PROPANOL 


In planning chemical research, it is customary first to make 
a thermodynamic analysis of the system in order to determine the 
feasibility of the reaction. However, once the Gibbs energy change, 
AG, is established for a reaction, its feasibility may readily be calcu- 


lated in terms of the equilibrium constant, K, by the following equation: 
AG=- RT InK. (ea) 


The AG value for the reaction can be computed from equation (F.2). 


AG = I[G, (products) ] - 2[G, (reactants) ] (Eee) 


where Ge is the standard Gibbs energy of formation. 

For the dehydration of 2-propanol at 252.7°C, G. (2-propanol), 
G_(propy lene) and G, (water) are - 22.55, - 4.33 and - 52.04 K cal/g-mole, 
respectively. These values were obtained by interpolating the data 
given in Reference (129). The AG value for the dehydration reaction 
is thus calculated to be - 33.82 K cal/g-mole, according to equation 
(F.2). The equilibrium constant, computed from equation (F.1), is 


SS or, in terms of partial pressure: 
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It is clear from the above expression that the equilibrium conversion, 
X, 1S essentially unity when Pro is less than one atm. 

Again in the case of the dehydrogenation reaction at 343.2°C, 
Ge (2-propanol) and Ge (acetone) are - 14.40 and - 18.72 K cal/g-mole, 
respectively. At the highest Pag» 52 cm Hg, used in this investigation, 
an equilibrium conversion of 0.98 is obtained by solving equation (F.4) 


fOr X. 


BRE sont2tioexs:. (F.4) 
Pro (52776) ~=—-‘1-x : 


From these results it appears that the assumption that both dehydration 


and dehydrogenation reaction are irreversible reactions is reasonable. 
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Isopropoxide 


Pareyrylateo” 


Hydrogen-bonded 
2-propanol: 


*Abbreviations: i? 


S 
m 
W 
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LIST OF IR BANDS ON THE ADSORPTION 


OF 2-PROPANOL ON y-ALUMINA 


Frequency, cm 


strong 
medium 
weak 
broad 


1130 
1170 
1340 
1380 
1460 
2870 
2960 


1380 
1470 
1575 


2960: 


1170 
1373 
1465 
2870 
2930 
2960 
3550 
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Intensity* 
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